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•Experiments  conducted  with  RP-1  in  a  Convective  Combustion  Environment 
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Oscillating  Liquid  Spray  Combustion  (continued) 
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Intrinsic  Instabilities  Identified  for  Diffusion  Flames 

•  Relevant  to  Nonlinear  Phenomena  in  Liquid-Propellant  Rocket  Instability 

•  Activation-Energy  Asymptotics  Applied  to  Derive  Dispersion  Relations 

•  Lewis-Number  Effects  are  Important;  Related  to  Cellular  Flames 

•  Rate-Ratio  Asymptotic  Analyses  are  Needed  in  Future  Work 
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TEN  YEARS  OF  RESEARCH  ON  H,-0,  INJECTION 
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YEARS  OF  RESEARCH  ON  H-Oo  INJECTION 
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ABSTRACT 

In  coimection  with  cost  reduction  and 
improvement  of  combustion  efficiency  of 
LOX/Hydrogen  or  LOX/Methane  engines, 
employment  of  swirl-coaxial  injectors  seem  to  be 
wise,  because  with  swirl-coaxial  injectors,  even 
employing  large  injection  elements,  which  are 
easy  to  manufacture  and  thus  can  decrease 
production  cost,  high  combustion  efficiency  may 
be  expected.  To  establish  the  design  methodology 
of  swirl-coaxial  injectors,  a  comparative  study  of 
swirl-coaxial  injectors  and  shear-coaxial  injectors 
has  been  conducted.*’^’ Steady  state  combustion 
characteristics  were  examined  with  multi-element 
injectors  with  a  nominal  thrust  of  10  kN  and  a  14 
kN,  and  a  chamber  pressure  of  3.5  MPa. 
Combustion  characteristics  were  evaluated  in 
relation  to  chamber  pressure  distributions,  wall 
temperature  distributions,  and  ■nc*s. 

The  main  findings  on  the  steady  state 
characteristics  were  as  follows: 

1)  Shear-coaxial  injectors  require  a  longer 
chamber  to  complete  vaporization  as  compared 
with  those  of  swirl  injectors;  therefore,  the  qc*  of 
shear-coaxial  injectors  is  apt  to  be  limited  by 
vaporization  efficiency.  2)  Combustion  with  swirl 
injectors  starts  near  the  injector  face  and  ends 
within  a  short  chamber  length.  The  qc*  of  swirl 
injectors  is,  therefore,  apt  to  be  limited  by  mixing 
efficiency.  3)  Injection  velocity  of  fuel  affects 
secondary  atomization  of  LOX  spray.  That  is,  the 
higher  the  fuel  injection  velocity,  the  finer  the 
size  of  LOX  droplets.  Fine  atomization  promotes 
combustion  near  the  injector  face. 

Because  of  the  remarkable  differences  of 
steady  state  combustion  characteristics, 
combustion  stability  characteristics  are  also 
expected  to  be  significantly  different  from  each 
other.  Stability  characteristics  were  examined  by 
introducing  a  bomb  disturbance  during  steady 
state  combustion.’’  LOX/hydrogen  and 
LOX/metihane  were  used  to  see  the  effect  of  fuel 
density  on  the  stability  characteristics.  Although 


all  tests  were  judged  to  be  stable,  the  differences 
of  stability  characteristics  of  shear-coaxial  and 
swirl-coaxial  injectors  were  evaluated  based  on 
amplitudes  and  frequencies  of  the  induced 
pressure  disturbances. 

The  main  findings  on  the  stability 
characteristics  were  as  follows. 

1 )  The  magnitudes  of  the  induced  maximum 
pressure  disturbances,  P'max/Pc,  decrease  with  an 
increase  of  injection  velocity  ratio,  VfTVo,  for 
both  injector  types.  2)  The  magnitudes  of 
P'max/Pc  with  swirl-coaxial  injectors  were  shown 
to  be  smaller  than  those  with  coaxial  injectors  at 
the  same  VfiVo.  3 )  The  correlations  between 
P'max/Pc  and  VWo  were  grouped  in  a  single 
curve  regardless  of  the  difference  of  propellant 
combinations  for  coaxial  injectors  and  for 
swirl-coaxial  injectors,  respectively.  4)  A 
plausible  reason  for  the  observed  trends  is  the 
interaction  between  the  shock  wave  induced  by  a 
bomb  and  LOX  droplets. 

A  cooperative  study  between  NAL  and  DLR 
is  now  on  going  to  see  the  combusting  flame 
during  unstable  combustion.  The  rectangular 
chamber  whose  acoustic  frequency  is  the  same 
order  with  a  typical  large  rocket  engine  will  be 
designed.  The  chamber  has  12  injection-elements 
with  glass  windows  which  allows  optical 
observation.  Bomb  disturbance  is  planed  to  be 
introduced. 
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Purpose  of  the  Study  and  Method  of  Approach 

Establishment  of  design  methodology  for  swirl-coaxial  injectors 
with  large  injection  elements. 

Steady  state  characteristics  were  evaluated  in  relation  to 
chamber  pressure  distributions,  wall  temperature  distributions,  and 
Tic*s  with  multi-element  injectors  with  nominal  thrusts  of  14  kN  and 
10  kN  and  a  chamber  pressure  of  3.5  MPa. 

Stability  characteristics  were  studied  by  introducing  a  bomb 
disturbance  during  steady  state  combustion.  LOX/hydrogen  and 
LOX/methane  were  used  to  see  the  effect  of  fuel  density.  Stability 
characteristics  were  evaluated  based  on  amplitudes  and  frequencies 
of  induced  pressure  disturbances. 
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a:  Non-recessed  swirl-coaxial 


b;  Recessed  swirl-coaxial 


Fig.  2  Comparison  of  recessed  and  non-recessed 
swirl-coaxial  injectors 

Swirl-coaxial  injectors  with  a  recessed  LOX  post  showed  a 
self-pulsation  of  the  spray.  To  avoid  instability  related  to  this 
self-pulsation,  recessed  injectors  with  swirl-coaxial  was  not  used. 


PARIS-France  May  27-28,  1999 


a:  Details  of  14  kN  injector  b:  Details  of  10  kN  injector 

Fig.  4  Injector  configuration 
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a:  Shear-coaxial  injectors  b:  Swirl-coaxial  injectors 

(  10  kN  injectors  ) 

Fig.  5  Comparison  of  axial  pressure  distribution 

A  pressure  drop  in  a  cylindrical  section  of  the  combustor  is 
caused  by  the  momentum  pressure  loss,  thus  pressure  distribution 
indicates  the  progress  of  combustion. 
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Fig.  6  Comparison  of  L50s  as  a  function  of  mixture  ratio 

The  length  required  for  chamber  pressure  drop  to  the  mean  value 
of  the  chamber  pressure,  L50,  is  a  measure  for  the  required  chamber 
length  for  complete  combustion. 
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Comparison  of  swirl-coaxial  and  shear-coaxial  injectors 

Shear-coaxial  injectors 

Swirl-coaxial  injectors 

Coarse  LOX  spray 

Vaporization  limited 

Large  VfA^o  improves  rjc* 

Combustion  starts  far  injector  face 

Low  heat  load  near  injector  face 

Require  longer  chamber 

Fine  LOX  spray 

Mixing  limited 

Small  YffVo  improves  qc* 

Combustion  starts  near  injector  face 

High  heat  load  near  injector  face 

Shorter  chamber  is  possible 
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20kHz 


Dynamic  pressure  trace  Pulse  spectrum 

a:  Coaxial  injector:  HI,  Pc=3.44  MPa,  0/F=6.16,  VfA^o=29.1,  RDX=  2  grain 


Dynamic  pressure  trace  Pulse  spectrum 

b:  Swirl  injector:  H2,  Pc=3.47  MPa,  O/F=6.03,  Vfi^o=4.11,  RDX=  2  grain 

Fig.  9  Dynamic  pressure  trace  and  its  pulse  spectrum 

(LOX/H2,  10  kN  injectors) 
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a;  Coaxial  injector;  Cl,  Pc=3.59  MPa,  O/F-3.69,  VfA^o-12.7,  RDX  2  grain 


Dynamic  pressure  trace 


Pulse  spectrum 


b:  Swirl  injector:  C3,  Pc=3.66  MPa,  0/F=3.57,  VfA^o=1.25,  RDX=  2  grain 

Fig.  10  Dynamic  pressure  trace  and  its  puise  spectrum 

(LOX/CH4,  10  kN  injectors) 
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a:  14  kN  chamber  b:  10  kN  chamber 

Fig.  11  Induced  maximum  pressure,  P’max/Pc, 
vs  Injection  velocity  ratio,  VfA^o 
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--mode  dB/ms 


a:  Dynamic  pressure  trace 


b:  Eliminate  low  frequency  components 


— - - - - - - ^ - - -  0  10  20 

c:  Hilbert  transformation  d:  Obtain  decay  rates  using  log-scale 

Fig.  12  Decay  rate  of  induced  pressure  disturbance 
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a:  Decay  rates  of  IL-mode  b:  Decay  rates  of  IT  mode 

Fig.  13  Correlation  between  decay  rates  of  induced 
pressure  disturbance  and  injection  velocity  ratio 

(  LOX/CH4,  coaxial  injectors,  14  kN  chamber ) 
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Semmary  of  stability  rating  test  results 

1 )  Induced  pressure  disturbances,  P'max/Pc,  decreases  with  an  increase 
of  injection  velocity  ratio,  Vf/Vo. 

2 )  The  P'max/Pc  with  swirl-coaxial  injectors  are  smaller  than  those 
with  coaxial  injectors  at  the  same  Vf/Vo. 

3  )  The  correlation  between  P'max/Pc  and  Vf/Vo  are  able  to  be 
grouped  in  a  single  curve  regardless  of  the  difference  of  propellant 
combinations  for  coaxial  and  for  swirl-coaxial  injectors,  respectively. 

4 )  Observed  correlations  indicate  existence  of  the  interaction  between 
the  shock  wave  induced  by  a  bomb  and  LOX  droplets.  That  is,  the 
larger  the  LOX  droplets,  the  more  the  sensitivity  to  the  shock 
disturbance. 
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Current  program 


A  cooperative  study  between  NAL  and  DLR  is  now  on  going.  The 
purpose  of  the  study  is  visualization  of  combusting  flame  during 
unstable  combustion.  The  optically  accessible  rectangular  chamber,  its 
acoustic  frequency  is  the  same  order  to  a  typical  large  rocket  engine,  is 
designed.  The  chamber  has  12-swirl  and  non- swirl  coaxial  injection 
elements. 


W.  Mayer  and  H.  Tamura 
"Flow  Visualization  of  Super¬ 
critical  Propellant  Injection  in 
a  Firing  LOX/OH2  Rocket 
Engine,"  The  Journal  of 
Propulsion  and  Power,  vol.  12 
No.  5,  1996. 


Our  previous  study  on  flow  visualization  in  a  LOX/BQ  rocket 
chamber  (LOX/H2,  Pc=4.5  MPa,  Vo=30  m/s,  Vf=300  m/s,  do=  1mm) 
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ABSTRACT : 

Three  experimental  investigations  have  been  done  on  the  cryogenic  test  bench  MASCOTTE  at  ONERA 
Palaiseau  center  with  the  objective  to  study  the  LOX  spray  characteristics  in  combustion. 

Informations  concerning  the  breakup  process  at  0.1  and  1  MPa  were  obtained  by  the  mean  of  a  visualization 
technique.  Quantitative  and  qualitative  description  of  the  spray  were  then  obtained  by  the  use  of  a  Phase  Doppler 
system.A  map  of  droplets  size  and  velocity  has  been  realized.  The  influence  of  the  momentum  ratio  J  on  the  first 
atomization  process  has  been  also  shown. 

Recently,  with  the  objective  to  analyse  the  local  Weber  number  by  the  mean  of  discrimination  liquid  and  gas 
velocity,  a  last  experiment  has  been  performed  with  a  2D  LDV  system. 

All  those  results  participate  for  a  better  understanding  of  atomization  in  cryogenic  sprays  in  combustion. 


NOTATION 

J 

momentum  ratio  (gas/liquid) 

Vm 

axial  mean  velocity 

Do 

initial  droplet  diameter 

me 

gas  mass  flowrate 

te 

evaporation  time 

niL 

liquid  mass  flowrate 

ts 

residence  time 

MPa 

pressure  (in  MPascal) 

evaporation  constant 

Vl 

liquid  velocity 

L 

combustor  length 

Vg 

gas  velocity 

Vg 

droplet  mean  velocity 

We 

Weber  number 

INTRODUCTION 

Cryogenic  propellant  combustion  systems  have  been  now  applied  successfully  for  many  years  in  rocket  engines 
(Ariane  1  to  5  ,  SSME,  ..,).  However,  to  improve  the  performance  of  this  type  of  combustion  systems  and  get  a 
good  modelling  of  the  combustion,  we  need  a  better  understanding  of  the  elementary  mechanisms,  in  particular 
atomization  and  vaporization  processes. 

A  lot  of  works  has  been  performed  last  10  years  to  study  the  behaviour  of  many  sprays.  All  those  studies  were 
more  often  conducted  with  simulant  fluids  (water  for  the  main  off  them).  It  is  of  course  because  it’s  more  difficult 
to  performed  those  studies  on  real  fluids  with  real  hot  fire  conditions. 

To  contribute  answering  those  questions  a  cryogenic  test  bench,  called  MASCOTTE,  was  built  and  a  test 
program  was  conducted  under  a  consortium  of  laboratories  and  manufacturers  associating  ONERA,  CNRS, 
CNES  and  SEP.  The  main  objective  was  to  provide  a  good  description  of  atomization  and  vaporization  processes 
under  hot  fire  conditions  and,  by  this  way,  to  have  a  better  understanding  of  this  kind  of  combustion. 

This  article  describes  experiments  performed  during  combustion  tests  mainly  focused  on  the  LOX  droplets  size 
analysis. 

First,  with  the  goal  to  provide  preliminary  informations  on  the  fluid  dynamics  of  the  LOX  core  breaking  process, 
a  visualization  technique,  using  a  stroboscopic  laser  sheet  associated  with  a  high  speeed  camera  (2000  ffames/s), 
was  used. 

After  those  initial  visualizations  the  main  objective  of  this  work  consisted  in  measuring  LOX  droplets  size  and 
velocity  on  the  axial  and  radial  flow  to  get  a  good  mapping  of  droplets  size  and  velocity  in  the  combustion 
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chamber.  Two  pressures  have  been  perform,  0.1  and  1  MPa.  A  PDPA  system  with  one  component  of  velocity 
was  used  for  those  measurements. 

The  two  experimental  investigations  (at  O.l  MPa  and  1  MPa)  give  a  good  qualitative  (by  visualization)  and 
quantitative  (by  Phase  Doppler  technique)  description  of  the  LOX  spray  in  term  of  atomization  and  and 
vaporization  process  in  hot  fire  conditions.  At  0.1  MPa  three  different  zones  has  been  detected:  near  the  post  tip  a 
first  atomization  zone,  after  a  relative  long  ligaments  zone  and  finally  a  well  defined  droplets  spray. 

A  comparison  between  LOX  droplets  size  and  velocity  at  0.1  and  1  MPa  give  us  informations  on  the  influence  of 
the  pressure  on  the  atomization  process. 

The  influence  of  the  J  number  (momentum  ratio)  on  the  atomization  has  also  been  performed  at  1  MPa. 

A  last  characterization  has  been  perform  recently,  for  instance  in  cold  flow  conditions,  with  the  objective  to 
analyse  the  local  Weber  number  by  the  mean  of  discrimination  of  liquid  and  gas  velocity. 


OPTICAL  SETUP 

Three  optical  configurations  have  been  used  during  those  experiments.  On  the  picture  of  the  figure  1  the  PDPA 
setup  is  shown  on  the  test  bench.  Like  is  well  know  today  the  PDPA  is  a  size  and  velocity  particles  optical 
measurement  system.  The  PDPA  complete  theory  has  been  described  and  explained  in  details  in  reference  [1]. 
The  schematic  drawing  of  figure  2  shows  the  optical  configuration  used  for  the  visualization  tests.  This 
visualization  has  been  performed  by  the  mean  of  a  stroboscopic  laser  sheet  associated  with  a  high  speed  camera. 
The  laser  sheet  was  formed  by  a  stroboscopic  copper  laser  beam  (at  A.  =  510.6  nm)  focused  through  a  cylindrical 
and  spherical  lens.  This  sheet  was  about  70  mm  wide  and  1  nun  thick  near  the  chamber  axis.  The  strobe 
frequency  was  2  kHz  and  the  pulse  time  duration  about  40  ns.  A  high  speed  16  mm  movie  camera  (NAC  E  -  10) 
was  used  in  association  with  black  and  white  high  speed  sensibility  movie  film  (500  ISO).  This  camera  was  used 
at  20(X)  frames/s  and  synchronised  via  an  optical  pick-up  to  the  strobe  frequency. 

On  the  picture  of  figure  3  the  optical  2D  LDV  system  used  recently  for  a  2D  velocity  analysis  of  both  LOX  and 
gas  hydrogen  is  shown. 

DROPLET  SIZE  AND  VELOCITY  RESULTS 

The  operating  conditions  are  listed  in  the  table  1.  Those  conditions  were  used  for  the  two  first  campaings  in 
combustion  (drop  size  and  visualisation). 


Case 

Pressure 

MPa 

Vl 

m/s 

Vg 

m/s 

j 

M 

We 

10® 

Re 

10® 

A 

0.1 

2.21 

15.00 

893 

13.4 

3.30 

0.598 

0.67 

A-10 

1 

2.23 

23.70 

308 

14.5 

2.10 

0.20 

0.67 

C 

0.1 

2.21 

10.00 

628 

6.3 

5.00 

0.304 

0.67 

C-10 

1 

2.23 

15.80 

207 

6.5 

3.20 

0.09 

0.67 

Table  1 :  Operating  conditions 


For  all  cases  the  LOX  mass  flowrate  was  maintained  constant  at  50  g/s. 

Results  at  0.1  MPa  12. 3, 41 

The  LOX  spray  evolution  in  the  combustion  chamber  has  been  describe  quantitatively  and  qualitatively.  A  map 
of  the  LOX  droplet  sizes  and  velocity  has  been  done  on  all  the  length  of  the  combustion  chamber.  One 
measurement  point  (on  the  axe)  has  been  done  also  outside  of  the  confining  chamber  and  shows  that  droplets 
exist  yet  in  this  zone  far  from  the  front  flame  and  cross  it.  Comparison  between  residence  time  (defined  as  :  4  =  L 
/  Vg),  and  evaporation  time  (defined  as  te  =  (Dq)^  /  by  A.  H.  Lefebvre  [5])  shows  that  at  0.1  MPa  droplet  size 
up  to  200  \im  was  partially  evaporated. 

A  comparison  with  computed  empirical  correlations  [6,  7,  8]  shows  that  the  droplet  size  measurements  are  in  the 
same  order  that  results  give  by  the  empirical  correlations.  For  more  informations  see  the  reference  [4]. 

An  «  artistic  drawing  »  of  the  qualitative  description  of  the  spray  is  shown  on  figure  4. 

Results  at  1  MiPa 

Like  above,  for  the  0.1  MPa  campaing,  a  qualitative  and  quantitative  description  of  the  spray  has  been  done. 
Those  complementary  results  have  allowed  to  do  a  comparison,  at  the  same  cross  section  of  the  chamber, 
between  0.1  MPa  and  1  MPa.  The  results  at  1  MPa  show  a  lower  number  density  but  droplet  size  (Dio,  D32)  are 
bigger  than  at  O.lMPa  (figure  5).  The  measurement  validation  rates  are  also  better  at  1  MPa.  This  shows  that  at 
high  pressure  all  the  small  droplets  are  evaporated  very  quickly  near  the  post  tip  of  the  injector  and  after  only  big 
structures  and  droplets  are  present  in  the  spray. 
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A  droplet  size  evolution  versus  J  number  has  also  demonstrated  the  importance  of  the  momentum  ratio  for  the 
first  atomization  process  (figure  6). 

VISUALIZATION  RESULTS 

Results  at  0.1  MPa 

During  this  first  campaing  only  the  zone  near  the  injector  has  been  analysed  (from  the  injector  post  tip  to  70  mm 
backward).  A  macro  description  of  the  breakup  process  has  been  done  (figure  7)  and  has  allowed  to  give  an 
evolution,  on  4  ms,  of  an  attach  ligament  to  simple  droplets  (spherical  or  not). 

Results  at  1  MPa 

The  aim  of  this  test  campaing  at  1  MPa  was  to  follow  the  LOX  structure  (if  she  exists)  along  the  combustion 
chamber.  Like  a  first  result,  a  comparison  with  results  at  atmospheric  pressure  shows  that  the  flame  emission  was 
clearly  detectable  at  high  pressure  and  not  at  0.1  MPa.  On  figure  8  two  pictures  recorded  during  the  test 
campaing  are  presented.  First  the  phenomena  were  recorded  without  any  filter  to  reject  the  luminous  flame  light 
(a).  Here  it  is  easily  to  see  the  LOX  structure  inside  the  flame.  In  a  second  configuration,  a  bandpass  filter, 
centered  at  the  wavelength  of  the  laser  beam,  around  510  nm,  was  fitted  in  front  of  the  camera  lens.  In  this 
configuration  only  the  scattered  light  of  the  liquid  structure  was  visualized  (b). 

The  analysis  of  the  pictures  at  same  axial  positions  confirms  the  PDPA  results.  The  spray  is  denser  at 
atmospheric  pressure.  Pictures  far  from  the  injector  (X/d  >  30)  shows  only  big  structures. 

LDV  TWO  PHASES  FLOW  RESULTS 

A  last  campaing,  perform  in  1999,  has  allowed  to  get  a  better  description,  in  term  of  velocity,  of  a  coaxial  LOX 
stucture  with  an  anulus  gas  flowrate.  The  2D  LDV  measurements  has  been  done,  for  instant,  in  none  reactif  flow 
(LOX  and  gas  N2  and  He).The  main  objective  was  to  get  the  local  Weber  number  by  the  mean  of  velocity 
discrimination  between  LOX  droplets  and  gas.  For  gas  velocity  analysis  a  seeding  system  was  added  to  to  the  gas 
distribution  panel.  The  particles  were  chosen  as  zirconium  oxide,  with  the  objective  to  be  used  later  in 
combustion  tests. 

A  first  theoretical  analysis  has  been  done  by  Onofri[9]  to  find  the  best  optical  configuration  of  the  2D  LDV 
system.  This  analysis  has  given  a  numerical  optimisation  based  on  the  diffusion  of  the  two  kind  of  particles 
(droplets  and  solid  particles).  For  solid  particles  less  than  1  pm,  used  for  seeding  the  gas  mass  flowrate,  and 
angular  zone  between  64°  and  1 14°  seem  to  be  the  best  optical  configuration. 

For  instance  the  first  results  show  an  important  recirculation  zone  near  the  injector.  An  other  important  result 
shows  that  the  seeding  solid  particles  stay  in  the  gas  zone  and  never  reach  the  two  phases  flow  in  the  LOX  spray. 
At  the  present  time  a  software  is  on  development  with  the  objective  to  be  able  to  analysed  the  raw  signals  and 
extract  from  the  intensity  levels  of  each  particle  the  velocity  of  LOX  droplets  and  seeding  particles  of  the  gas. 

CONCLUSION 

This  paper  has  presented  the  main  results  got  during  a  combustion  test  program  performed  on  the  ONERA 
cryogenic  test  bench  MASCOTTE.  Three  experimental  investigations  were  performed  until  today.  Those 
investigations  were  mainly  focused  on  the  LOX  droplet  size  and  velocity  analysis.The  visualization  technique 
give  us  first  informations  on  the  LOX  structure  and  particularly  the  difference  between  tests  at  atmospheric 
pressure  and  1  MPa.  By  the  use  of  a  PDPA  measurement  technique,  a  good  mapping  of  droplet  size  and  velocity 
around  the  LOX  core  has  been  obtained. 

As  main  conclusions  today  we  can  say  that  we  have  got : 

a  good  qualitative  (by  visualization)  and  quantitative  (by  PDPA  technique)  description  of  the  LOX 
sprays ; 

those  results  provide  us  informations  on  the  atomization  process  under  hot  fire  conditions,  with  a 
comparison  between  atmospheric  pressure  and  1  MPa ; 

the  last  step  is  yet  under  analysis,  but  would  be  able  to  help  us  to  understand  the  local  spray 
evolution  by  the  knowledge  of  the  local  Weber  number. 


ACKNOWLEDGMENT 

This  work  was  supported  by  CNES  and  SEP  under  contract  N°  875366  and  876691  ;  the  authors  would  like  to 
express  their  thanks  for  this  financial  support. 

REFERENCES: 


57 


[1]  W.  D.  Bachalo,  M.  J.  Houser,  Phase  Doppler  spray  analyser  for  simultaneous  measurements  of  drop 
size  and  velocity.  Optical  Engineering,  volume  23,  N°  5,  1984. 

[2]  P.  Gicquel,  L.  Vingert  D.  Lourme^  L.  Menoret,  Coaxial  Injectors  Atomization,  Progress  in 
Astronautics  and  Aeronautics,  vol.  169,  Chap.6,  pp.  145-189,  1995. 

[3]  P.  Gicquel,  L.  Vingert  et  E.  Brisson,  Caracterisation  experimental  d'un  brouillard  LOX-GH2  en 
combustion  issu  d*un  injecteur  coaxial.  Actes  du  5  ieme  «  Colloque  du  GDR  Combustion  dans  les 
moteurs  fusees  » 

[4]  P.Gicquel,  E,  Brisson,  L.  Vingert,  Experimental  investigation  of  a  LOX  spray  under  hot  fire 
conditions.  ICLASS  97,  Seoul  (Korea),  August  18-22,  1997. 

[5]  A,  H.  Lefebvre,  Atomization  and  Sprays,  Hemispher  Press,  NY,  1989 

[6]  S.  Nukiyama  and  Y.  Tanasawa,  Experiments  on  the  atomization  of  liquids  in  an  air  stream,  Trans. 
Soc.  Mech.  Eng.  Jpn.,  vol.  5,  N°  18,  pp.  62-67,  1939. 

[7]  E.G.  Lorenzetto  and  A.H.  Lefebvre,  Measurements  of  drops  size  on  a  plain  jet  airblast  atomizer, 
AIAA  Journal,  vol  2,  N®  3,  pp.  583-585,  1964. 

[8]  K.  Y.  Kim  and  W.  R.  Marchal  Drop  size  distributions  from  pneumatic  atomizers,  Journal  of  Am. 
Chem.  Eng.,  vol  17,  N°  3,  pp.  575-584, 1971. 

[9]  F.  Onofri,  Optimisation  numerique  de  la  configuration  experimentale  a  adopter  lors  de  la  campagne 
de  mesure  des  champs  de  vitesse  particulaire/fluide  du  foyer  LOX/GH2  «  MASCOTTE  »  Rapport 
d’etape  CNRS  du  3  sept.  1998 


58 


Windows 


Fig  2.  Optical  visualization  setup 


Fig  3.  2D  LDV  system  optical  arrangement 


Fig  4.  Qualitative  description  of  the  spray  at  0.1  MPa 
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Fig  7.  Evolution  of  the  breakup  process  during  4  ms 


(a) 


(b) 


Fig  8.  Visualization  of  flame  and  spray  (a)  and  spray  only  with  an  additional  green  filter  (b) 
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Quantitative  Density  Measurements  in  Supercriticai 
LN2-Sprays  Using  Raman  Scattering 
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conditions  in  a  cryogenic  high  pressure  rocket  combustor: 
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Deutsches  Zentrum  fur  Luft-  und  Raumfahrt  e.V.  German  Aerospace  Center 
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workshop  May  27-28  1999.DOC 


Studies  on  the  Break-up  Regimes 

of  LOk  Droplets 
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37‘^  AIAA  Aerospace  Sciences  Meeting  and  Exhibit  January  11-14, 1990  /  Reno,  NvI 


Introduction 
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Conclusions 


Vttteatn  OMjtiat  injector 
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Gaseous  hydrogen 


^^ctivites  du  L.C.S.R. 


Experimentale  :  en  configuration  monogoutte  : 


goutte 

initiation  de 

gouttes 

fin  de 

primaire 

la  cassure 

secondaires 

vaporisation 

QJ! 

c  » 

o 

•  • 

*  ■  • 

Temps 

d’initiation 


regime  de 
cassure 


duree  de  |  temps  (taux)  de 
cassure  vaporisation  des 

gouttes  secondaires 


temps  total  de  cassure  nombre  et  taille  dcs 

gouttes  secondaires 


Experimentale :  Spray  inerte 


Numcrique  :  vaporisation  /  combustion  goutte  Isolde 
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PARAMETRES  INrLUANj:ANTS 


Temps  d^initiation  etrcgime^e  cassure  : 

-  Reynolds 
-Weber 

-  Ohnesorgc 

-  pression  (via  rapport  de  densite  gaz/goutte) 


Duree  de  la  cassure  (temps  total  -  temps  d^lnitiation) : 

-  Weber 

-  report  de  densites 

independant  de  la  pression 


Nombre  et  taille  des  gouttes  secondaires  : 

-  Reynolds 
-Weber 

-  pression 

-  Do  (diametre  initial  goutte  primaire) 


Temps  de  vaporisation ; 


-  Revnolds  turbulent 


k'“d 


-  Schmidt 


-  ^midt 


turbulent 


'12 


-  Lewis 

-  Pression 
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Droptet  -  jet  interaction 
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Diagnostics  and  intake  processing 


particle  image  velocimetry  for  the  characterisation 
of  the  aerodynamic  field. 


Exploration  area  bt^  f/ideo  camera 

Pjl  Cryogenic 
_ i  injector 
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Conditions  explored  for  the  aerodynamic  fie 

characterisation 
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Conditions  exptored 


o 

Z>  b 

D) 

Cl 

li 
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UJ 


■quid  oxygen/  33.3500  ^.^30  320. 3300 


Superppsitim  of  the  aerodynamic  field  and 

droplet  hreak*up  imayes 
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Droptet  deformation  and  break- up  regimes 
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increasing  I  2  :  Bag  break-up  regime 
Weber  number  I  3  :  Transitionai  break-up  regime 

▼  4  :  Shear  break-up  regime 


Images  for  the  LOx  /  Or^  air  i 

P  =  0,1  MPa,  DR  =  1200,  We  =  48,  Re 


Images  for  the  LOx  /  Hetium  couple 

P  =  2  MPa,  DR  =  350,  We  =  19,  Re  =  1350 


eV\A  ‘Jaqujnu  jaqa^ 
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Neuf  correlation  for  the  regime  transitions 


^  CO  M  T-  o 
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Pressure  (MPa) 


Normalised  break-up  initiation  time 

Break-up  initiation  time  measurements  with  images 
processing. 
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Weber  number,  We 


Normatised  hreahup  duration 


uoDBjnp  dn-)|eajq  pasiieuuoN 
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100 

Weber  number,  We 


Nombre  et  taille  des  gouttes  secondaires 


-  difficile  a  determiner  experimentalement 

-  resultats  des  differents  auteurs  peu  concordants 


Quelques  resultats  au  L.C.S.R.,  mais  plage  de  parametres  restreinte 
->  pas  de  correlation,  mais  tabulation  et  interpolation  possible 


OAygene  iigaide  /  heiium  (P  ss  0^  MPa  :  We  ^  137 :  Rc  ^  2570  :  Dq  *  /  39  mm) 


oxysene  liquide  f  helium  (P  =  OJ  MPa  ;  We  =  78 :  Re  *  1940  :  1J90  mm) 
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and  viscosity  ratio. 


Cone  t US  ions  (2) 
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ADVANCED  AIRBREATHING  AND  ROCKET  ENGINES 
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Subcritical  Droplet  Vaporization 


Supercritical  Droplet  Gasification 


•  Transient  Diffusion 

•  Thermodynamic  Non-ideality 

•  Property  Variation 

•  Ambient  Gas  Solubility 
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Thermophysical  Properties  and  Constitutive  Relations 
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EQUATIONS  OF  STATE 
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Temperature,  K 


Evaluation  of  Thermodynamic  Properties 

•  sensible  enthalpy  :  h(p,  T)  =  ho(T)  +  Ahexc  (p,  T) 

•  internal  energy  ;  u(p,  T)  =  uo(T)  +  AUexc  (p,  V 

•  specific  heat  :  Cp(p,  T)  =Op(T)  + ^Op^^^(p,  T) 

where 

h^(T),u°(T),Cp(T)  =  values  in  dilute -gas  limit 
Ahexc >  ^^exc^  ^^Pexc  corrections 


Pressure-explicit  type  of  EOS: 


Ah, 


'exc 


-  fp 
"  Jo  ‘ 


p 


2Ut. 


'PJ 


w+Rnz-^) 


AU, 


exc 


_  fp|  _P _ Lf ^ 

Jo 


p^UT/pj 


bp 


AC 


=  -T 

Pexc  Jo  p2 


^  ^2  _  ^ 


iZB 

dr 


bp+ 


Tidp/dTY 


p^{dpldp)T 


-R 
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Thermal  Conductivity,  W/(m  K) 


Transport  Properties  of  a  Pure  Substance 


TEMPERATURE.  K 


Mp»T)  —  ^(T)  +  AXgjjQ{p,T)  +  AX,Q|.j^(p,T) 

where  =  dilute  gas  limit 

AXgxc  =  dense  fluid  correction 

AX,crit  =  critical  enhancement 
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THERMAL  CONDUCTIVITY.  mW/cm  K 


SPECIFIC  HEAT,Cp,(Bto/lb®R)  THERMAL  CONDUCTIVITY,  Blo/h-hr 


THERMAL  CONDUCTIVITY  OF  OXYGEN 


TEMPERATURE,  K 


SPECIFIC  HEAT  OF  OXYGEN 


TEMPERATURE,  K 


SPECIFIC  HEAT,  Cp,  J/g 


DENSITY,  mol/L 


Critical  Enhancement  of  Transport  Properties 


•  Critical  Region  Proper :  I  AT  I  <  0.03,  I  Ap  I  <  0.25 

at  =  {T-Tc)/Tc  and  Ap*  =  (p  -pc)/pc 

•  Extended  Critical  Region 


•  Critical  Enhancement  of  Thermal  Conductivity 


Density,  g/cm 


PROPERTY  EVALUATION  FOR  MIXTURES 

•  Extended  Corresponding-State  Principle 

-  Density  based  scheme 

~  no^Poi"^ 0^^ n 

-  BWR  EOS  for  density  at  given  pressure  and 
temperature 

-  Mapping  function  determined  by 

•  critical  properties 

•  PItzer's  acentric  factor 

-  Reference  fluids  based  on  methane  or  propane 


126 


/’/(Dab^) 


Mass  Diffusivity  for  Multicomponent  Mixture 


•  Determine  Binary  Mass  Diffusivity 

—  Chapman-Enskog  theory  with  Lennard-Jones 
intermolecular  potential  energy  function 

—  Corrected  by  the  corresponding  state  approach 
ofTakahashi 


0.0  1.0  20  3.0  4.0  5.0  6.0 


Reduced  pressure,  Pr 


•  Apply  the  mixing  rule 

Di„,  =  (1-Xi)/:E(XjDij),  i=1 . N 
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Droplet  Vaporization  and  Combustion  in  Quiescent  Environments 
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THEORETICAL  FORMULATION 


•  Complete  Time-Dependent  Conservation  Equations 
for  Both  Gas  and  Liquid  Phases 

—  mass 

—  momentum 

—  energy 

—  species  concentration 

•  Finite-Rate  Chemical  Kinetics  for  Gas-Phase  Reaction 

•  Matching  of  Droplet  Behavior  with  Ambient  Gases 

—  Continuities  of  mass  and  energy 
—  Phase  equilibrium 
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VAPORIZATION  OF  LOX  DROPLET  IN 
HYDROGEN  FOR  PRESSURES  OF  5  -  250  atm 


CRITICAL  PROPERTIES  OF  OXYGEN  AND  HYDROGEN 


Tc  (K) 

Pc  (atm) 

Vc  (cm^/g) 

Zc 

OXYGEN 

154.6 

49.74 

2.294 

0.288 

HYDROGEN 

33.2 

12.96 

32.292 

0.306 
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Time  History  of  Droplet  Surface  Temperature 


Distribution  of  Mixture  Specific  Heat 


Reduced  Radial  Distance,  r  /  R^ 
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Droplet  Lifetime,  ms 


Effect  of  Pressure  on  Droplet  Lifetime 

500  <  <2000 K,  5<p^<  SOOatm 


•  Enthalpy  of  vaporization  decreases  with  pressure 

=>  Enhancement  of  vaporization  in  subcritical  regime 

•  Critical  mixing  temperature  decreases  with  pressure 

+  Gaseous  thermal  conductivity  increases  with  pressure 
=>  Enhancement  of  vaporization  in  supercritical  regime 
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Correlation  of  Droplet  Lifetime 


•  Processes  dominated  by  thermal  diffusion  in  supercritical  regime 

•  Simplified  analysis  based  on  thermal  diffusion  equation  in  terms  of 

-  thermal  diffusivity  a 

-  initial  droplet  radius  squared  R] 

-  transfer  number  (r„  -  T^)l {T^-T^) 

•  Correlation  of  reduced  droplet  lifetime 


■ 

^  T  ^ 

0.0115  +  0.542 

2  ^  oo  ^  c 

f 

/ 

T^-Tg  ^ 

where  / 


=  1  +  3.9 


f 


1  -  exp 


-0.035 


V 


1 

)) 


(T.-T,)/(T.-To) 
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Droplet  Lifetime,  ms  ^  Droplet  Lifetime,  ms 


Effect  of  Thermal  Diffusion  on  Droplet  Lifetime 
(Critical  Mixing  Temperature) 


Reduced  Pressure,  p  /  p^ 


of  Thermal  Diffusion  on  Droplet  Lifetime 
(Critical  Mixing  Composition) 


Reduced  Pressure,  p  /  p. 
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Temperature,  K 


Sp6ci6S  Composition  in  LOX/Hj/HgO  System 


Temperature  Distribution  in  LOX/Hj/HjO  System 
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Droplet  Lifetime,  ms 


Droplet  Lifetime,  ms 


Effect  of  Pressure  on  Droplet  Lifetime 

400  <T^<  2500  K,  5  <p„<  300  atm 


Subcritical  Regime 

•  Low  ambient  temperature 

Molecular  diffusion  decreases  with  pressure 
=>  Nitrogen  composition  decreases  at  the  interface 
=»  Surface  temperature  increases 
=>  Temperature  gradient  decreases  at  the  surface 
=>  Droplet  lifetime  increases  with  pressure 

•  High  ambient  temperature 

Enthalpy  of  vaporization  decreases  with  pressure 
=>  Enhancement  of  vaporization  with  pressure 

Supercritical  Regime 

•  Mixture  critical  temperature  decreases  with  pressure 

+  Gaseous  thermal  conductivity  increases  with  pressure 
Enhancement  of  vaporization  in  supercritical  regime 
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Correlation  of  Droplet  Lifetime 

•  High  ambient  temperature,  >  1000 /iT 

•  Various  initial  droplet  temperatures,  =  300, 325, 350K 

•  Reduced  droplet  lifetime  correlated  with  pressure 

- =  exp  ] 


where 

a  =  0.0161  for  n-heptane/air  system 
a  =  0.01  for  n-decane/air  system 
•  Average  deviation  <  2  % 
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Effect  of  Pressure  on  Vaporization  Rete  of  UDMH  Droplet 


S/_LULU  ‘lUBlSUOQ  UOUBZIJOdBA 
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Pressure,  MPa 


COMBUSTION  OF  n-PENTANE  FUEL  DROPLETS 
IN  AIR  FOR  PRESSURES  OF  5-140  atm 


CRITICAL  PROPERTIES  OF  CONSTITUENT  SPECIES 


n-C.sHip 

Critical  Temp. (k)  469.6 

Critical  Pressure  (atm)  33.3 


Q2 

CO? 

H?0 

126.2 

154.6 

304.2 

647.3 

33.5 

49.8 

72.8 

217.6 
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Time-Variations  of  Droplet  Surface  Temperature 


Time-Variations  of  Latent  Heat  of  Vaporization 
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Peak  Flame  Temperature  (K) 


Effect  of  Pressure  on  Peak  Flame  Temperature 


•  Thermal  diffusivity  decreases  almost  linearly  with  pressure, 

a  ~  =►  less  heat  transfer  from  the  flame  zone. 

•  Mass  diffusivity  decreases  relatively  fast  with  pressure, 

D  <  =►  much  less  fuel  supply  to  the  flame  zone. 

•  Lewis  number  increases  with  pressure, 

Le  =  -^  T  as  P  T 
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MILESTONE  TIMES  WITH  DROPLET  BURNING 

MECHANISMS 


Do  =  100 


0  20  40  60  80  100  120  140 

PRESSURE,  otm 


At  low  pressures,  gasification  dominates. 

At  high  pressures,  transient  gas-phase  diffusion  dominates. 


h“ 

CO 

3 

m 

s 

o 

o 

Q 


n  ^ 
2  Z 

&  2 
GC  Ij 


o 

0. 

< 

> 

h" 

UJ 

Q. 

o 

cc 

Q 

U. 

O 

<n 

LU 

</) 

z 

o 

CL 

(/) 

Lii 

o: 

o 


> 

o 


o 

CO 

o 

o 

u. 

1- 

z 

m 

m 

s 

< 

o 

H 


0 

(0 

c 

o 

Q. 

W 

0 

CE 

■D 

0 

Q. 

3 

O 

O 

I 

0 

D 

W 

« 

0 


0 

<Ji 

C 

o 

Q. 

W 

0 

CE 

TJ 

0 

Q. 

D 

O 

O 

i 

’o 

o 

0 

> 


144 


Droplet  Vaporization  Response 


M7M 

P7P 
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de  of  R 


Magnitude  of  Droplet  Vaporization  Response 


Phase  Angle  of  Droplet  Vaporization  Response 


Abnormal  Variations  of  Thermophysical  Properties  at  Critical  Point 
Cause  Sharp  Amplication  of  Droplet  Vaporization  Response 
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THEORETICAL  FORMULATION 
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,  m/s 


SUPERCRITICAL  DROPLET  GASIFICATION  IN 
CONVECTIVE  ENVIRONMENTS 

•  Droplet  Gasification 

•  Momentum  Transfer 

•  Droplet  Dynamics 

Key  Parameters 

•  Reynolds  number,  Re=p„\J^D//j„ 

•  Transfer  number,  B  =  (Ypc-Yp^j/d-Ypc) 

•  Lewis  number,  Le  =  a/D 

•  Ambient  density,  p„,  (or  ambient  pressure) 
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LOX  DROPLET  VAPORIZATION  IN  SUPERCRITICAL  HYDROGEN 
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LOX  DROPLET  VAPORIZATION  IN  SUPERCRITICAL  HYDROGEN 
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LOX  DROPLET  VAPORIZATION  IN  SUPERCRITICAL  HYDROGEN 
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DIFFERENT  MODES  OF  DROPLET  EVOLUTION 


Evolution  of  Critical  Surface 


155 


Critical  Mixing  Temperature  -  Critical  Mixing  Composition 


Evolution  of  Critical  Surface 


/  \ 
/  I 


156 


Supercritical  LOX  Droplet  Vaporization  in  Convective  Environments 
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Critical  Mixing  Temperature  -  Critical  Mixing  Composition 


Effects  of  Pressure  and  Velocity  on  Droplet  Lifetime 
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1/(1+0.09*Re^’*(p/p^g,) 


Droplet  Drag  Coefficient 
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Supercritical  Droplet  Combustion  in  Convective  Environments 
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REACTION  MECHANISMS 


O 


CM 


CM 


•D 

C 

CC 

w 

o 

X 

x“ 

d 

x’ 

o 


CM 


■D 

i 

O 

> 

c 

0) 

O 

'o 

o 

CL 

CO 

sz 

O) 

Lu 


O) 

c 

■  MM 
> 

§ 

c 

c 

o 

‘■p 

o 

cs 

0) 

oc 


+  X 
O 

t  o 
2  t 

+  X 

o  + 

+  O 
X  X 


o  5 

+  o 

o  + 

CM  W 

X  o 

t  t 

X  o 

+  + 

CNJ  CM 

o  o 

X  X 


CM 

o 

+  X 

O  CM 

V'O 

-1-  CM 

A  X 

O  CM 

X  X 

+  + 
CM  CM 

o  o 

X  X 


CM 

o 

X  X 

+  o 
o  + 

X  t 
O  X 

+  + 
O  CM  CM 

5  0  0 


o 

CM 

X 

O) 

c 

> 

o 

> 

c 

c 

o 


a> 

QC 


CM  CM 

X  X 


+ 
cmX 
^  O 
+  + 

o  5 

X  o 

t  f 

X  S 

+  + 

CM  CM 

o  o 

CM  CM 

X  X 


o 

(0 

o 

QC 

c 

o 

(0 

VS 

c 


X 

+ 
CM 

o 

X 

t 

CM 

S  O 

+  + 
CM  CM 

X  X 


+ 

X 

+ 

X 

t 


•o 

c 

(0 

z 

o 

O) 

c 

■  MB 
> 

o 

> 

c 

c 

o 

vs 

o 

(0 

0) 

a 

c 

'« 

o 


z 

o 

ID 

c 

(0 

z 

o 

O) 

c 

■■■» 

> 

§ 

c 


X 

X 

o 

c 

A 

X 

X 

+ 

+ 

■S 

‘J3 

+ 

+ 

CM 

o 

O 

o 

o 

CM 

X 

(0 

c 

CM 

X 

+ 

+ 

X 

o 

E 

t 

t 

o 

X 

f 

t 

1 

1 

t 

I 

I 

o 

T 

T 

X 

O 

1- 

+ 

+ 

CM 

o 

CM 

X 

o 

CM 

X 

c 

« 

X 

o 

+ 

+ 

+ 

CM 

+ 

+ 

+ 

X 

O 

X 

O 

O 

X 

o 

X 

o 

t 


CM 


X 

t 


+ 

+  X 

o  o 
+  + 
X  X 


161 


162 


163 


FILTERED  CONSERVATION  EQUATIONS 
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HIGH-PRESSURE  SPRAY  FIELD  DYNAMICS 
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GRID  DISTRIBUTION 
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BASELINE  FLOW  CHARACTERISTICS 


y.  m 


E 


50 

25 

0 


u-Velocity,  -20  <  u  <  250  m/s,  A  =  1 0  m/s 


mm 


1  atm,  0/F  =  4 


Spray  Distribution,  2576  Parcels,  6.3x10"^  Particles 


Vorticity,  0  <  co/co^gf  <  1 00%,  A  =  1 0%,  co^gj  =  -74x1 0  s 


100  atm,  0/F  =  4 


Spray  Distribution,  2680  Parcels,  2.8x10®  Particles 


Mass  Fraction,  0.4032  -  Yh2  <  1 ,  A  =  0.05 


Viscosity  Ratio,  0  ^  Vj/vi  <  35.52,  A  =  5 


400  atm,  0/F  =  4 


Spray  Distribution,  1363  Parceis,  5.1x10®  Particles 
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SHEAR  COAXIAL  UNI-ELEMENT  INJECTOR 
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Table  5:  Inlet  conditions  employed  for  the  analysis  of  high-pressure  hydrogen-oxygen  mixing  and 
combustion  processes.  Subscript  1  denotes  the  hydrogen  stream,  subscript  2  denotes  the  oxygen 
stream.  -  - - 


1  2  3  4  5  6 


p  [atm] 

100 

100 

100 

100 

100 

100 

Ti  [A'] 

150 

150 

150 

300 

300 

300 

T2  [K] 

100 

100 

100 

300 

300 

300 

T2/T1 

0.6667 

0.6667 

0.6667 

1 

1 

1 

Pi  [kg/m^] 

15.37 

15.37 

15.37 

7.752 

7.752 

7.752 

P2  [kg/mp] 

1118 

1118 

1118 

134.7 

134.7 

134.7 

P2/P1 

72.75 

72.75 

72.75 

17.37 

17.37 

17.37 

Ui  [m/s] 

500 

250 

125 

250 

125 

62.5 

U2  [m/s] 

30 

30 

30 

30 

30 

30 

u2/m 

0.06 

0.12 

0.24 

0.12 

0.24 

0.48 

P2U2/P1U1 

4.365 

8.730 

17.46 

2.084 

4.169 

8.338 

P^Vl/PiVl 

0.2619 

1.048 

4.190 

0.2501 

1.001 

4.002 

Ml 

0.4828 

0.2414 

0.1207 

0.1786 

0.08931 

0.04465 

M2 

0.04572 

0.04572 

0.04572 

0.09203 

0.09203 

0.09203 

Res, 

1424000 

715700 

359400 

243200 

121900 

60990 

Res^ 

245000 

245000 

245000 

201700 

201700 

201700 

Pri 

0.7842 

0.7842 

0.7842 

0.6913 

0.6913 

0.6913 

Pr2 

1.917 

1.917 

1.917 

0.7832 

0.7832 

0.7832 

Sci 

1.560 

1.560 

1.560 

1.365 

1.365 

1.365 

SC2 

89.09 

89.09 

89.09 

0.2207 

0.2207 

0.2207 

Lei 

1.964 

1.964 

1.964 

1.975 

1.975 

1.975 

Le2 

46.48 

46.48 

46.48 

0.2817 

0.2817 

0.2817 

Res  =  UeS/i>  where  Uc  represents  the  edge  velocity  and  S  the  channel  half  width. 
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Density,  kg/rn"  H  Velocity  Vectors 
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Figure  19:  Contours  of  density,  temperature,  velocity,  and  vorticitydn  the  near-field  region  of  the  splitter  plate  for  inlet  conditioii-3dn 


Density,  kg/m  ■  Velocity  Vectors 
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figure  21:  Contours  of  density,  temperature,  velocity,  and  vorticitydn  the  near-field  region  of  the  splitter  plate  for  inlet  condition-^dn 


Mass  Fraction  ■  HoO  Mass  Fraction 
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Figure  22:  Contours  of  H2,  O2,  H2O,  and  OH  mass  fraction  in  the  near-field  region  of  the  splitter  plate  for  inlet  condition  5  in  Table-5u 


PLANAR  COAXIAL  TRI-ELEMENT  ARRAY 
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INTRODUCTION  : 
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Lahoratoire  E.M2.C 


SPRAY  MODEL  :  KIVA3  adapted  models 
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Lahoratoire  E.M2.C. 


SPRAY  MODEL  :  KIVA3  adapted  models 
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COMBUSTION  MODEL  :  Coherent  Flame  Model 
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Laboratoire  E.M2.C. 


COMPUTATIONAL  DOMAIN  :  Configuration  and  mesh 
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COMPUTATIONAL  DOMAIN  :  Configuration  and  mesh 
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Laboratoire  E.M2.C. 


pCULATIpN  qONDITIONS!;  general  inflow  conditions 


SIMULATION  CONDITIONS  :  4  different  inflow  conditions 
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14.5 


COMPARISON  WITH  EXPERIMENT  :  LIF/simulation 


COMPARISON  WITH  EXPERIMENT  :  LIF/simulation 


Laboratoire  E.M2.C. 


RESULTS  :  Point  A,  with  recess 


4^ 

CENTRALE 


Lahoratoire  E.M2.C. 

195 


RESULTS  :  Point  A,  with  recess 


Point  A 


Water  mass  fraction 


2 
1 
0 
-1 
-2 

Oxygen  mass  fraction 


CENTRALE 


Laboratoire  E.M2.C. 
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RESULTS  :  Point  A,  with  recess 


Point  A 

Radial  velocity  (cm.s”^) 


CENTRALE 


Lahoratoire  E.M2.C. 
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Point  A  with  recess  Point  A  without  recess 
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Resultats  numeriques  :  Taux  de  reaction. 


RESULTS  :  Point  A,  without  recess 


Laboratoire  E.M2.C. 


COMPARISON  WITH  EXPERIMENT  :  PDA/simulation 


Laboratoire  E.M2.C. 


COMPARISON  WITH  EXPERIMENT  :  PDA/simulation 
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COMPARISON  WITH  EXPERIMENT  :  PDA/simulation 
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COMPARISON  WITH  EXPERIMENT  :  PDA/simulation 


Laboratoire  E.M2.C. 


DISCUSSION  :  reaction  rate,  spra 
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Thin  Film  Dynamics 
and  the 

Physics  of  Atomization 


Dirk  Meinkohn 

Deutsches  Zentrum  fiir  Luft-  und  Raumfahrt  DLR 
D-74239  Hardthausen 
e-mail:  dirk.meinkoehn@dlr.de 


Liquid  Rocket  Combustion  Chamber  Flow  Dynamics 
International  Workshop  CNES  (Paris),  May  27-28  (1999) 
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General  Remarks 


•  liquid  propellants:  combustion 

-  gasification 

-  mixing  processes 

-  combustion:  efficiency 

-  motor  operations:  stability 

•  spray  states:  intermediate,  important 

•  feedback  loop: 

atomization  ^  spray  combustion 


•  atomization  problematics:  no  comprehensive  model  available 

•  atomization  hydrodynamic  process 

-  liquid  membranes/ligaments:  flow  processes 

-  liquid  edges:  formation/dynamics 

-  droplet  detachment:  neck  formation/pinching/rupture 

•  asymptotic  descriptions 
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Atomization  Processes 


•  forcing  process: 

-  nozzle  flow  — jet 

-  air  assist/electrostatic/ultrasonic/ . . . 

-  generation  of  liquid  membranes/filaments 

•  decomposition  of  membranes/filaments 

-  droplets/spray  formation 

•  investigation:  applications  oriented/empirical  approach 

-  forcing  process:  technicalities/devices 

-  spray  properties:  empirical 

-  feed-back  from  downstream  combustion:  neglected 

•  investigation:  physics  orientated/predictive  approach 

-  liquid  flow  processes  in  membranes/filaments:  2-d/l-d 

-  interaction  with  3-d  environment: 

*  formation/dynamics  of  liquid  edges 

*  criticality  static  equilibrium  of  liquid  edges) 

*  propagation  of  perturbations  in  liquid  membranes 

*  rupture/necking/pinching 

*  feedback  processes 

•  deductive  approach  ->  predictive  approach 
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147.  Kelvin-Helmholt:  waves  on  a  thin 
liquid  sheet.  Water  is  ejected  downward  in 
the  form  of  a  thinning  sheet  from  a  fen- 
spray  nozzle  that  is  vibrated  sinusoidally 
normal  to  the  plane  of  the  sheet  at  a  fre¬ 
quency  giving  wave  instability.  A  side  view 
shows  the  initial  exponential  growth  pre¬ 
dicted  by  two-dimensional  theory  and  the 
development  of  a  complex  wave  form  as  the 
amplitude  increases.  The  stream  of  drops 
along  the  axis  at  the  bonom  is  seen  in  the 
plan  view  to  result  from  collapse  of  the  rims 
of  the  sheet.  Cropper,  Dombrouiski,  Jepson  & 
Pyott  1973 


ctnrjpJUA'  ^ 


210 


Announcement 


Daimler-Benz  Aerospace 


Raumfahrt-^nfrastaiktur 


Third  French-German  Colloquium  on 
Research  on  Liquid  Rocket  Propulsion 

4  L'  <Y^'  ^ 


Topics  ; 

•  Heat  Transfer 

•  Diagnostics 

•  Modelling 

•  Atomisation 


5-6  June  1997 
Marseilles  (France) 


3-D  Simulation  of  liquid-gas  shear  layer.  S  ^eski  &  J.Li  ^ 

Ctt 

IZ 

(  prepared  by  the 

^  CNES  SEP  DASA  DLR  joint  programme  group 

organised  by  the  CNRS  (IRPHE)  in  partnership  with 
the  Universite  de  la  Mediterrannee  and  the  city  of  Marseilles 
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149,  Breakup  of  a  liquid  sheet.  Studies  of  drop  forma¬ 
tion  in  sprays  show  that  the  production  of  filaments  gener¬ 
ally  constitutes  an  intermediate  stage.  Controlled  hiamen: 
formation  should  therefore  provide  a  means  of  regulating 
the  size  of  drops.  One  approach  is  to  subject  the  nozzle  to 


forced  vibrations.  This  photograph  shows  a  flat  laminar 
sheet  of  water  issuing  from  a  fan-sprav  nozzle  chat  is  oscil¬ 
lated  a.xially  at  resonant  frequence.  P.horogrurh  by  ,V. 
Dorr/i:rou‘ski 
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Physical  Characterization 


/ 


•  flow  processes  ^  hydrodynamics 

1.  liquid  membranes/filaments:  deformation  propagation  of  pres 
sure  perturbations 

2.  propagation  of  existing  free  edges 

3.  rupture  processes  for  membranes/filaments 

4.  formation  of  free  edges 

•  long  wave  theory: 

reference  wave  length  X  large  ->  reduction  in  dimensionality 

**  2-d  membranes  /  1-d  filaments  ** 

•  characteristic  dimensions:  cm  ,  sec  ,  g 

•  characteristic  variables 


(relative)  flow  velocity 

u 

[cm/ sec] 

characteristic  length 

h 

[cm] 

density 

P 

[g/cmr'] 

kinematic  viscosity 

V 

[cm// sec] 

surface  tension 

o 

[g/secp-] 

•  characteristic  groups 
Reynolds  Re  =  uh/\ 

Capillary  Ca  =  upv/G 

Weber  We  =  c/hpu^  = 

Ohnesorge  Oh  =  v-y/p/o/i  = 

Bond  Bo  =  pgh}/G 
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Hydrodynamic  Models 


•  Navier-Stokes  equations  ^  continuum-theoretic 


V  •  =  0  ,  p  =  const 


3^'6  +  ('6- V)a)  =  iv/>-f-v(V- V)!) 


•  thin  liquid  membranes/ligaments 

p  =  OK 


•  2  limit  cases  (no  gravity): 


case  1 : 


Weber-number  hydrodynamics 


-  no  viscosity/no  dissipation 

-  propagation/dynamics: 

*  pressure  perturbations 

*  free  edges 


case  2: 


Ohnesorge-number  hydrodynamics 


-  viscosity  important/strong  dissipation 

-  rupture  processes: 

*  continuum-theoretic  asymptotics 

*  needle’s-eye-role  of  microphysics 
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Weber  Number  Hydrodynamics 


•  static  criticality;  drag  force  D  —  surface  tension  force  T 


•  typical  estimate:  diameter  of  maximal  drop  dmax 

reference  diameter  do 


•  dynamics:  propagation  of  perturbations 


•  velocity  v  =  {fCc}  =  —  V$  -4  velocity  potential  $ 

•  wave  equation  (Taylor): 

•  antisymmetric  waves: 

U  =  aWasI^  Smh{Kz)  COs(k[;C  —  Wagt]) 

c  =  aWasK'  cosli( k::) CO .s(k[x  — 


•  wave  speed  Was  — >  dispersion  relation 
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Geometric  Acoustics:  Rays 


Liquid  Sheets:  Edge  Dynamics 


•  cylindrical  edge; 

-  retraction  velocity  =  v 

-  tension  force  (per  unit  length  of  edge): 

T  =  f  (/i/2)  dor  2a-/h)  sin(a)  =  2(7 
J  0 

-  mass  (per  unit  length  of  edge):  dM  =  phdx  — ^  M  =  phx 


•  stationary  edge  — >  geometric  acoustics 


•  conclusion: 

1.  importance  of  Weber  number  for  edge  dynamics 

2.  analogy  Weber  number  with  Mach  number  (gas  dynamics) 

3.  ray  tracing,  method  of  characteristics 
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Droplet  Detachment:  Necking/Pinching  Phenomenon 


•  acoustic  theory:  inertia,  surface  tension,  no  viscosity 


characteristic  time: 


W  =  h^pl(T  <  1 


•  droplet  separation:  necking/pinching 


_ . .  1 

\  >  / 

I 

.1 - 1 

. . . ^ 

. . . r 

i 

1  ^ 

- 


—  large  curvature  k  ~  l//i,  large  surface  tension  forces 

—  well-posedness:  dissipation/ viscosity  necessary 

—  Ohnesorge  hydrodynamics 


•  modelling:  full  Navier-Stokes  +  continuity  — >■ 

inertia  +  surface  tension  +  viscosity 


viscous  time  scale: 


•  (Keller/Miksis  1983,  Eggers/Dupont  1994)  pinching  — > 

asymptotic  phenomenon 

—  long  wave  approximation 

—  large  neck  curvature  — >  boundary  values  negligible  self-similar  solution 

—  new  scales  for  length  and  time:  surface  tension  as  well  as  viscosity 
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TA.  Kowalewski/ Fluid  Dynamics  Research  17  (1996)  121-145 


t  —  4.0/i5  t  =  IZfis 


Fig.  3.  Separation  of  the  droplet  from  the  liquid  jet  of  radius  rj  =  98.5  yum.  Time  is  given  relative  to  the  moment  of  the  jet 
breakup.  Liquid  medium  of  relatively  low  viscosity  (MIXD).  The  width  of  one  frame  corresponds  to  1  mm 
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J.  Eggers  and  T.  F.  Dupont 


Figure  2  Simulation  of  a  drop  of  water  suspended  from  a  circular  orifice  of  radius  =  0.26  cm.  This 
aives  parameter  values  =  0.992  and  /..//■„  =  4.S9  x  10  *.  compare  table  1.  The  time  distance 
between  profiles  is  0.4(r*p/7)^  starling  from  a  point  where  the  drop  is  already  falling.  We  also 
superimpose  a  profile  at  the  snap-off  and  the  corresponding  e.xperimental  picture  from  Peregrine 
et  at.  (1990).  There  is  no  adjustable  parameter  in  the  comparison.  To  enhance  contrast,  we  erased  the 
background  in  the  original  photograph. 
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Conclusions 


•  physics  of  atomization: 

—  deformation  processes/instability  of  liquid  membranes/ligaments: 

Weber  number  hydrodynamics,  acoustics  hyperbolic  PDE 

—  pinching/rupture  of  liquid  membranes/ligaments: 

Ohnesorge  number  hydrodynamics  parabolic  PDE 

—  post-rupture '  processes:  edge  dynamics/secondary  rupturing/satellite 
droplets 

•  modelling  efforts; 

—  no  comprehensive,  all-encompassing  model  available  as  yet 

—  available:  partial  models  covering  asymptotic  aspects 

—  interaction/coupling  with  external  acoustics:  partially  possible 

•  goal:  physical  characterization  of  atomization 

—  process  duration 

—  droplet  size  distribution 

—  fully  coupled  wdth  external  acoustics  — >  elliptic  PDE 

•  compare:  state-of-the-art  in  phase  transition  modelling 

—  no  closed  comprehensive  model 

—  two  asymptotic  partial  models: 

*  nucleation 

*  late-stage  ripening/Ostwald  ripening 
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•Resulting  "spray  flame"  is  controlled  by  the  turbulent  scalar  mixing 
and  the  turbulent  dispersion  of  droplets. 


no  effep-i 


Laser  Mie  scattering  of  the  reacting  spray.  The  green  and  blue  colors  correspond  to  the 
droplet  scattering  and  the  flame  emission  respectively.  (Stepowski  et  al.,  1992) 


•The  random  positions  and  properties  of  evaporating  droplets 
intervene  in  the  turbulent  micro-mixing  in  the  gas  mixture. 

•It  modifies  the  probability  distribution  of  temperature  and 
concentrations  and  eventually  the  mean  local  rate  of  chemical 
reactions. 

•Models  of  turbulent  gaseous  flames  that  accounts  for  turbulence- 
chemistry  interaction  are  often  based  on  the  PDF-equation  method. 
The  objective  here  is  to  extend  the  PDF-equation  approach  to  the 
spray  combustion. 
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PDF-Equation  Approach  in  Multicontinua  Formulation: 
Example  for  vapor  and  liquid  mass  fractions 


Vapor  mass  fraction  signal 


Results  in 
joint  PDF 
P(Yv.YJ 


t 


Conditional  (Y^-0)  vapor  mass  fraction  signal 


t  |Yl=0 


Results  in 

conditional 

PDF: 

P(Yv,YJ  = 
f(Y^)*6(Y,-0) 
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From  the  stochcistic  differential  equation  to  the  PDF-equation 


Yt 


A  • 


unconditional  PDF  P{Xv) 


time 


•  Fluctuations  (l^-,  T)  in  the  gas  are  caused  by  the  random  velocity  field 
and  by  the  random  rate  of  vaporization/condensation.  The  set  of  all 
events  constitutes  a  sample  space. 

•  Instantaneous  continuity  for  the  vapor  mass  fraction  : 

dY 

p-^  +  •  vn  =  V  (DpVK)  +  pci,,  (i-y„)  +  pci,j,„ 


•  Instantaneous  continuity  for  the  gas  mixture  density  : 

|  +  V.(pf)=M. 

•  Vaporization  rate  :  uig  =  PvVvap,s-^', 

-4  E  =  |Vyi,|  (J(l^  —  Y^s)  -  instantaneous  interface  per  unit  volume 
yvap,s  -  vapor  velocity  relatively  to  the  displacement  of  interface 
(jig  =  0  if  Kap,»  =  0  (or  Yy  Yys) 

‘jckem  -  chemical  reaction  rate 
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•  Unconditional  PDF  balance  equation  : 


$/»#CC. 


+  V  [(/)f  |y„=y;)/^n)]  =  {pu,s\y„=y„)p{y„) 

-Jr  [fl(y.)(w,iK,=n>  (i-n)] 

-Jr  [i^n)(/**’cie„|K=v;>] 


•  Integration  over  all  spectrum  of  leads  to  the  averaged  continuity 
equation  : 

dt 


+  V  -  (pF)  =pUJlg 


•  Vaporization/condensation  rate  is  explicitly  available  from  Lagrangian 
tracking  of  the  droplets 

•  Non-linear  chemical  reaction  term  is  presented  in  exact  form 
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Unconditional  joint  PDF-equation  for  gas  mixture 
(gaseous  species  and  gas  enthalpy)  : 


^{p)P{Yu...,YM,h)  +  v[{pV)P(Yu...,Y„,h)\  =  {pu,,)P{Yu...,YM,h) 

-  [((V  (pDVYi)}  -  {pu„)Yi  +  p(Yu  ...,Ym,  ft)] 

,=i  oYi'- 

-  ^  [((V  ipDVh))  +  ipq,,))  />(k„  y„,  ft)] 

g^^P^cherT^P{^\.i '"i^M 


•  Enthalpy  source  (or  sink)  due  to  the  droplets  qig  is  explicitly  available 
from  Lagrangian  tracking  of  the  droplets  ;  9/^  =  0  if  the  heat  exchange 
at  the  interface  =  latent  heat 

•  Alongside  with  the  CFD,  the  modeling  of  PDF-equation  supplies  each 
cell  with  the  probability  distributions  in  order  to  calculate  the  mean 
combustion  rate 
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Micro-mixing  modeling 
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Laminar  and  turbulent  parts  in  the  mixing 


mtxt 


VS 


At  what  conditions  entirely  turbulent  mixing  occurs? 


For  stagnant  droplet: 

^  P  liq  .  k, 

^evap  ^  ’  ^tur  ^  > 

SD  Pgas 

^evap  P  liq 

- cc - —  - - 

^tur  8t|  P  gas  k 


AU 


ocRe 


-1/4 


OC 


1 


P  liq 
P  gas 


JlV  ^  1  P  liq  ^tur  . 

SRe^'^^  P  gas  ^evap 


do>n 


mixing  is  turbulent. 


For  mobile  droplet: 


Ud  -  u 


gas 


mixing  is  turbulent. 
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Simulation  of  convective  transport 
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Simulation  of  diffusive  transport 
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ray  in  the  constant  volume  bomb,  (inert  case) 
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Experimental  and  numerical  inlet  setup 
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Computations  of  temperature  and  vapor  distributions 
using  Eddy-Break-Up  combustion  model  (middle)  and  PDF  approach  (right) 


239 


Computed  results  Stepowski  et  al.,  1994; 

~  ,  „  .  I  OH  planar  laser-induced  fluorescence  images 

.  ^  mixing  controlled  zone  ® 

mean  combustion  of  the  reacting  spray.  Images  represent  two 
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Computed  probable  zone  of  stabilized  Stepowski  et  al.,  1994 


■g 
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OVERVIEW 


Rocket  Injector  Design  Criteria 

Propellant  Combination  j 

=>  Performance,  Stability  and  Reliability  Aspects 
Low  Cost 

Focused  Application-oriented  Injector  Research  at  Penn  State 

=>  LOX-rich  Preburner  -  LOX/GH2  Propellants 

Orbital  Maneuvering  System  (OMS)  Upgrade  -  LOX/EthanoI  Propellants 
=>  Gas/Gas  Injector  Technology  -  GO^/GH,  Propellants 

Injector  Flowfield  Research  at  Penn  State 

=>  Dynamics  of  Impinging  Jet  Injectors 

=>  Chamber  Mach  Number  Effects  on  Injector  Flowfield  -  GO2/GH2  and 
LOX/GH2  Propellants 

Summary/Future  Plans 
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ROCKET  INJECTOR  DESIGN 


Application  Defines  Propellants/Injector  Design 


LOX/GH2 

•  SSME 
•X-33 

•  LOX  Rich 
Preburner 


PROPELLANTS 


•  L0X/GH2 

•  GO2/GH2 

•  LOX/Ethanol 

•  LOX/RP-1 

•  LOX/H,0, 


LOX/H2O. 

•  Upper  Stage  Engine 


LOX/RP-1 

•X-34 


GOj/GH, 

•  Full  Flow  Engine  Cycle 

•  RBCC  Rocket 


LOX/Ethanol 
•  OMS 
•RCS 


pennState 
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INJECTOR  RESEARCH  AT  PENN  STATE 


•  Focused  Application-oriented  Injector  Research  at  Penn  State 

=>Gas/Gas  Injector  Technology  -  GO2/GH2  Propellants 
— >  Shear  Coaxial  Injector 
— >  Swirl  Coaxial  Injector 
Impinging  Injector 

=>  LOX-rich  Preburner  -  LOX/GH2  Propellants 
->  Swirl  Coaxial  Injector 

Orbital  Maneuvering  System  (OMS)  Upgrade  -  LOX/Ethanol  Propellants 
— >  Pintle  Injector 

Bi-centrifugal  Coaxial  Swirl  Injector 
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GAS/GAS  INJECTOR  TECHNOLOGY 


Motivation 

•  Data  Base  on  Gas/Gas  Combustion  for  Rocket  Application  Is 
Limited  To  Studies  of  Global  Characteristics 

-  Aerojet  Report  (1973) 

•  Gas/Gas  Injector  Design  Requires  a  Data  Base  On  Local 
Phenomena 

-  Mixing  and  Combustion  Characteristics 

—  Injector  Face  and  Wall  Heat  Transfer  Issues 

•  Laser-based  Diagnostic  Techniques  and  CFD  Models  Have 
Matured,  and  Can  Be  Used  As  Design  Tools  for  Next-generation 
Technology  (RLV)  Development 


pennState 
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GAS/GAS  INJECTOR  TECHNOLOGY 


Raman  Species  Measurements  for  60°  Swirl  Injector 

O/F =6  H,  —  H 


Hj  Instantuncous 

H,  Average  of  8  RcaMzalions 

ufm^9wmmm 


HjO  Average  of  1 1  Realizations 


P^.=993  psia;  x=3.5  in. 


•Ur 


-'m  ,  .  'i  ^ 


H,0  Instantaneous 


I V  \  1 

..-h  ! 

'.V-.'  I  '  '  -I 

I  i 


.;()  ill  '»  HI  -H 

r(rrm) 

X=89  tnm 


Near-Complete  Mixing  and  Combustion  Within  90  mm 

pennState 


248 


GAS/GAS  INJECTOR  TECHNOLOGY 


Ht/Oj  Velocity  Ratio=8.0 


Hj/Oj  Velocity  Ratio=5.9 


pennState 


HVOj  Velocity  Ratio=6.8 


UJOj  Velocity  Ratio=3.9 


GAS/GAS  INJECTOR  TECHNOLOGY 

Summary 


Swirl  Coaxial  Injector 

-  Mixing  and  Combustion  Is  Extremely  Fast 

-  Heat  Release  Achieved  Within  a  Short  Distance 

-  Face  Temperature  Increases  With  Swirl  Angle 

-  Use  of  Swirl  Coaxial  Injector  Will  Require  Effective  Face  Cooling 
Strategies 

Shear  Coaxial  Injector 

-  Mixing  and  Combustion  Is  Slower  Than  Swirl  Coaxial  Injector 

-  Mixing  and  Combustion  Increases  With  GH2  to  GO2  Velocity  Ratio 
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LOX-RICH  PREBURNER 


Motivation 

•  Improve  Understanding  of  LOX/GH2  Combustion  at  LOX-rich 
Conditions 

•  Conduct  Detailed  Study  of  Swirl  Coaxial  Injector  Performance 

•  Data  Base  to  Validate/improve  Computational  Combustion  Tools 


LOX-RICH  PREBURNER 


RAMAN  SETUP 


ICCD 


Jl 


105  mm 
Nikkor  Lons 
f/#=1.8 


Holographic 

Imaging 

Spectrograph 

f/#=:1.8 


A  Beam 

I 


Stop 


IR 

Filler 


f=0.75  m 
Focusing  Lens 


Laser 

Window 

Imaging 

Mirror 


11 


Rocket 


INSTRLMENTCELL 


o 


Personal  Computer 


Gate/Pulse  Generator 


Data  Acquisition  System 


1 
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LOX-RICH  PREBURNER 

X=I27  mm,  X/d=37  (X=axial  distance,  d=LOX  post  ID) 
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LOX-RICH  PREBURNER 


X-127  miriy  X/d=37  (X^xial  distance,  d=LOX post  ID) 

Q9-1 - - - - - 

Q8  ^  ♦ 

|“-  ♦  ♦♦♦  ♦  ♦♦  ♦ 

I  “  V 


♦  Oxvaen  • 


Racial  lixation  (inch) 
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LOX-RICH  PREBURNER 


Summary 

Single  Shot  Linewise  Multi-species  Raman  Images  Of  LOX/GH, 
Rocket  Combustion  Obtained  For  Oxidizer-rich  Conditions 
Resultant  Temperature  And  Major  Species  Profiles  Yield  Expected 
Values 

Swirl  Coaxial  Injector  Mixedness  Is  Further  Evaluated 


pennState 
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OMS  UPGRADE 


Motivation 

NASA  Is  Investigating  the  Feasibility  of  a  Non-toxic  O.MS/RCS 
Upgrade  for  the  Space  Shuttle  to  Eliminate  the  Safety  and  Cost  Issues 
Associated  With  Hypergolic  Propellants 
Broad  Objectives 

To  Provide  NASA  and  Contractors  With  an  Objective  Base  From 
Which  Future  Injector  Design  Decisions  Can  Be  Made  for  OMS  and 
Future  HEDS  Applications 
Specific  Objectives 

To  Measure  and  Analyze  Detailed  Combustion  Characteristics  of 
LOX/Ethanol  Propellant  Combinations  at  Representative  Thrust 
Chamber  Conditions  for  Various  Candidate  Injectors 
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OMS  UPGRADE 
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OMS  UPGRADE 


OMS  Engine  Design 
Thrust;  27  kN  (6000  IbO 
Isp:  326  sec 
Mixture  Ratio:  1.7 

Propellant  Flowrate:  8.3  kg/s  (18,4  Ibm/s) 
-Oxygen:  5.26  kg/s  ( 1 1 .6  Ibm/s) 
-Ethanol:  3.08  kg/s  (6.8  Ibin/s) 


Sub-Scale  Engine  Design 
1/60  of  Full  Scale-mass  tlowrate 
Mixture  Ratio:  1.7 

Propellant  Flowrate:  0.14  kg/'s  (0.307  Ibm/s) 
-Oxygen:  0.088  kg./s  (0.193  Ibm/s) 
-Ethanol:  0.052  kg/s  (0. 1 14  Ibni/s) 


Case  # 

1 

2 

■Bn 

4 

5 

6 

P,  (MPa) 

l.IO 

1.10 

mmm 

2.07 

2.07 

2.07 

O/F 

2.0 

1.7 

1.4 

2.0 

1.4 

LO.X  tlowrate  (ks/s) 

0.095 

0.088 

0.082 

0.177 

0.155 

ethanol  tlowrate  (k2/s) 

0.048 

0.052  i 

0.059 

0.088 

0.096 

HuDi 

GN2  flowrate  (k 2/s) 

0.014 

0.014 

0.014 

0.014 

0.014 

0.014 

Nozzle  Diameter  (mm) 

16..M 

16..51 

16.51 

16.51 

16.51 

16.51 
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OMS  UPGRADE 


Over  100  LOX/Ethanol  Tests  Were  Conducted  on  Pintle  Injectors 


94  -j 


0 

0 


B  0° 


*  D3  Pintle  j 
0  K1  Pintle ' 


1  1.1  1.2  1.3  1.4  1.5  1.6  1.7  1.8  1.9  2  2.1  2.2 

0/F 
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OMS  UPGRADE 


Bl-CENTRIFUGAL  NEAR  INJECTOR  FLOWFIELD 


Outer  LOX  flow  is  separated  from  inner  ethanol  cone  by  layer  of  vaporized  oxygen 
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OMS  UPGRADE 


MODIFIED  BI-CENTRIFUGAL  INJECTOR  WITH  LARGE  RECESS 

•2M  w/Iarge  recess  face  plate  exhibits  good  mixing  and  combustion  over 
wide  range  of  operating  conditions  and  is  insensitive  to  LOX  quality 
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OMS  UPGRADE 


PERFORMANCE  OF  Bl-CENTRIFl'GAL  INJECTOR 


O  1  M  (90  inj/fire  plate)  low  P 

•  1  M  (90  inj/fire  plate)  high  P 

Lj  2  M  ‘,60  inj/fire  plate)  low  P 
■  2  M  (60  inj/fire  plate)  high  P 

9  2  low  P 

A  2  Delta  (narrow  gap)  low  P 
A  2  Delta  (narow  gap)  high  P 
5  2  M  w/big  recess  low  P 

□  2  M  w/big  recess  high  P 

O  2  0  (open  w/big  recess)  low  P 

♦  2  0  (open  w/big  recess)  high  P 

0  3  Delta  A  high  P 

V  3  B  (reverse  scheme)  low  P 
T  3  B  (reverse  scheme)  high  P 


0.7  0.9  1.1 


1.3  1.5  1.7  1.9 


LOX  Quality  (Pc/Py  at  T3) 
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OMS  UPGRADE 


Summary 

PINTLE  INJECTOR 

•  Pintle  Injector  Designed  for  OMS  Conditions  Operated  in  a  Stable 
Manner 

•  For  Pintle  Injector,  c*  Efficiency  Increase  With  O/F 
BI-CENTRJFUGAL  SWIRLER 

•  Marginal  LOX  Quality  Influenced  Combustion  Efficiency  for  All  Injector 
Configurations 

•  Shadowgraph  Imaging  Was  Instrumental  in  Assessing  Flow  Phenomena 
Which  Caused  Lower  Than  Expected  Performance 

•  Use  of  Recess  Is  Reasonable  Solution  to  Spray  Cone  Separation  Problem: 
-  Recess  Insensitive  to  LOX  Quality  (of  Significance  to  OMS  Upgrade), 
Stable  Under  .All  Conditions  Tested,  LOX  Cooling  of  Injector  Face 
Effective  Here 
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INJECTOR  RESEARCH  AT  PENN  STATE 


Injector  Flowfield  Research  at  Penn  State 

=>  Dynamics  of  Impinging  Jet  Injectors 
— >  Like-on-like  Impinging  Injectors 

=>  Chamber  Mach  Number  Effects  on  Injector  Flowfield  -  GOt/GH,  and 
LOX/GH,  Propellants 
->  Shear  Coaxial  Injectors 


pennState 


DYNAMICS  OF  IMPINGING  JET  INJECTORS 


Motivation 

•  Determine  Conditions  Under  Which  Impinging  Jet  Injector 
Periodic  Atomization  Can  Cause  Combustion  Instabilities 

=>  Longitudinal  Mode  Excitation  for  GOj/ethanoI  Propellants  in 
Model  Rocket  Chamber 

=>  Electromechanical  Drivers  Controlling  Individual  Ethanol  Jets 
(of  Impinging  Jet  Injector) 


pennState 


DYNAMICS  OF  IMPINGING  JET  INJECTORS 


•  Ethanol  Impinging  Jet  Injector 
(Individual  Electromechanical  Drivers 

•  Annular  Flow  of  GO2 

pennState 


258 


259 


DYNAMICS  OF  IMPINGING  JET  INJECTORS 

Summary 


•  Successful  Demonstration  of  Combustion  Instability  by  Periodic 
Atomization  for  Impinging  Jet  Injector 

•  Technique  Has  Potential  As  Active  Control  of  Combustion  Instability 


pennState 


CHAMBER  MACH  #  EFFECT  STUDIES 

Motivation 


•  Document  Effect  of  Chamber  Mach  #  on  Mixing  and  Combustion 
Characteristics  of  Shear  Coaxial  Injector  for  Both  GOi/GH,  and  LOX/GHj 
(Currently  in  Progress)  Propellants 

•  Obtain  Data  Base  for  CFD  Code  Validation 


pennState 
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CHAMBER  MACH  #  EFFECT  STUDIES 


•  GOVGHj  Study  Injector  ^ - 

•  O/F  =4  GOVGHj  Combustion  , 

•  Background  Air/GH,  Flow 

•  1200  K  Background  Flow 

•  Chamber  Pressure  -1.3  MPa 

(190  psia) 

•  Chamber  Mach  #  Varied 

=>for  M#  of-0,  0.1,  0.2,  0.3 


-.coo;  o c o 


Air  Holes 


I  coocooocooo 
1 cOocOo  oCo  oC 
cocoooccoc 


c  GO^/GHj  Shear  Coaxial  Injector 
,  Swirlers 
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CHAMBER  MACH  #  EFFECT  STUDIES 


1 

0.1 

0.093 

0.024 

0.730 

2 

0.2 

0.093 

0.024 

1.887 

0.016 

3 

0.3 

0.093 

0.024 

3.197  ^ 

0.030 

Case 

0/F  Main 

0/F  Background 

0/F  total 

pennState 
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CHAMBER  MACH  #  EFFECT  STUDIES 


Focused  Laser  Beam  Passes 
Vertically  Through  Rocket 
Chamber 


ICCD  Camera 


Imaging  Spectrograph 


pennState 


CHAMBER  MACH  #  EFFECT  STUDIES 


M#  =  0  M#  =  0.1 


pennStaje. 
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CHAMBER  MACH  #  EFFECT  STUDIES 
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Recent  Studies  on  Prebumer  Injector  Characteristics 
for  the  Hydrogen-Oxygen  Engine 

J.Wang,  S.Q.Cheng,  L.Z.Song,  RX.Hu 
11th  Research  Institute  of  China  Aerospace  Corporation 


Abstract  In  order  to  develop  the  technology  for  high  performance  hydrogen- 
oxygen  staged  combustion  rocket  engine,  four  kinds  of  preburner  injector 
have  been  studied  theoretically  and  experimentally.  By  using  the  Phase 
Doppler  Particle  Analyzer(PDPA),  the  spray  characteristics  of  these  injectors 
are  investigated,  through  varying  liquid-to-air  pressure  drop  ratio,  liquid 
mass  flow,  and  back  pressure  in  the  high  pressure  rig.  Fire  tests  have  also 
been  completed  for  demonstration.  The  spray  performance,  gas  temperature 
distribution,  combustion  stability  and  ignition  performance  for  different 
injector  configurations  are  revealed. 

1.  Introduction 

The  prebumer  is  a  key  component  for  the  staged  combustion  cycle  engine,  owing 
to  it's  higher  pressure  and  flowrate  than  that  of  the  gas  generator  cycle  engine.  It's 
combustion  performance  and  stability,  gas  temperature  distribution  at  exit  are  closely 
related  to  it's  injector  configurations.  In  order  to  develop  high  performance  prebumer 
injector,  recently  we  studied  four  kinds  of  injector  configurations  theoretically  and 
experimentally. 

Four  kinds  of  injector  configurations  are  depicted  in  fig.  1-3.  The  coaxial  injector 
with  swirler  of  liquid  oxygen  in  centre  can  get  high  performance  and  commonly 
used  in  china  Hydrogen  Oxygen  Rocket  Engine,  The  coaxial  injector  with  straight 
stream  of  liquid  oxygen  in  centre  use  the  successful  experiences  of  the  SSME/LE-7 
prebumer  injectors  and  it  has  a  well-earned  reputation  for  it's  reliable  operation  and 
processing;  The  bi-coaxial  injector  is  a  new  kind  of  injector  configuration  we 
designed  specially.  It's  hydrogen  sleeve  consist  of  two  segments,  hydrogen  and 
oxygen  bum  at  near  stoichiometric  mixture  ratio  in  the  first  segment.  In  the  second 
segment,  is  injected  the  remainder  of  hydrogen  which  mix  and  cool  the  stream  exited 
from  first  segment.  Turbulent  diffusion  combustion  is  transformed  into  premixture 
combustion  for  this  injector,  so  that  it's  combustion  rate  increases.  It's  spray  and 
mixture  performance  are  as  same  as  that  of  the  shear  coaxial  injector. 

Prebumer  injector  unit  is  made  of  injector  face  plate,  central  ignition  port  and  four 
injector,  as  shown  in  Fig.4.  The  ignition  port  diameter  is  16mm.  Four  injectors  with 
same  configuration  are  equally  located  at  a  circle  of  diameter  32mm.Four  injector 
configurations,  namely  A,B,C,D  are  used  respectively  in  their  injector  elements  for 
testing. 
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we  employed  the  Phase  Doppler  Analyzer(PDPA)  to  study  spray  droplet  size  and 
flowrate  distribution  for  different  injector  configuration,  using  air  and  water  as 
imitation  agents,  varying  liquid-to-air  injector  pressure  drop  ratio,  liquid  mass 
flowrate  and  back  pressure  in  the  high  pressure  rig.  Theoretical  approach  is  to 
calculate  the  atomization  combustion  model  which  consider  gaseous  phase,  droplet 
phase,  turbulence  and  combustion  effects.  Changing  the  measured  spray  data  into 
inlet  conditions  of  the  imitation  calculation,  we  can  get  the  spray  combustion 
efficiency  and  gas  temperature  distribution  and  other  parameters  for  the  four  kinds  of 
injector  configuration.  The  injector  elements  are  also  undertaken  fire  testings,  it's 
chamber  pressure  ranges  6.85-7.73  MPa,  it's  mixture  ratio  ranges  O.639-0.845. 
Combustion  efficiency,  gas  temperature  distribution  and  vibration  are  measured  . 


[21  mm 

Fig.3  bi-coaxial  injector 


Injector  arrangement 
Fig. 4  injector  unit 


2.  Test  Facilities 

Injector  spray  testing  goes  on  in  the  high  pressure  rig,  as  shown  in  Fig.  5.  Normal 
pressure  in  the  rig  is  3.5  MPa,  it's  internal  diameter  320mm,  highness  1200mm. 
There  are  four  windows  on  it's  flank,  inserted  with  plate  glass  of  thickness  25mm 
and  diameter  100mm,  used  for  PDPA  light  emission,  receive  and  monitor.  Safety 
valve  and  pressure  gauge  are  mounted  on  the  upper  cover,  the  regulator  on  bottom  is 
used  to  keep  chamber  pressure. 
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Fire  test  system  for  prebumer  injector,  in  Fig.6,  consist  of  hydrogen  flow  lines,  liquid 
oxygen  flow  lines,  nitrogen  pressurization  system,  nitrogen  cutoff  system  and 
measurement  system.  The  venturi  orifice  and  sonic  nozzle  are  respectively  set  on 
liquid  oxygen  flow  line  and  gaseous  hydrogen  flow  line.  The  measured  data  include 
propellant  flowrate,  inlet  injector  pressure,  gas  average  pressure  and  fluctuation 
pressure,  vibrations  corresponding  to  axial,  radial,  tangential  direction  at  injector 
element  head.  Three  temperature  gauges  are  mounted  on  exhaust  duct  to  measure  gas 
temperature  distribution. 


Fig.5  spray  measurement  system  Fig.6  fire  testing  system 


3.  Theoretical  Model 

the  spray  combustion  model  in  [4]  is  adopted  for  theoretical  calculation.  The  model 
consider  gas  phases  combustion,  as  well  as  droplet  behavior  in  gas  phases  and  their 
interaction.  These  processes  are  coupled  each  other,  so  that  the  calculation  of  the 
spray  combustion  model  should  be  done  by  solving  equation  system  of  gas  phase 
equation,  droplet  equation,  combustion  and  turbulence  equation. 


3.1  Gas  Phase  Governing  Equations 

We  can  get  the  gas  phase  governing  equations  concerning  the  continuity, 
momentum,  energy  and  composition  in  the  following  average  form: 


dp  d  .  - 


pm 


d  , - >,  d  f - .  d  ,  —  /  -/  X  ^ p  ,  ^ ^ ij  ,  “Tr  C 

>>  =  -  ax,-  dxi  +  ^ 

^  /”rr\  ^  / - rj  , - ^  ^  X 

+  =  17^  dJ.^Pr  dl^  ^ 

d  -  -  -  dj^  dq.- 

+  +S,H 
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(2) 


Turbulence  model  employs  k  -  e  equations,  namely 

■  j-ypk)  *  ||j  +  G»  -  ^ 

Assuming  that  stress  is  proportional  to  average  strain,  the  expression  of  Reynold's 
strain  is  as  following: 


- ; — 7-  r^ii;  -i  A  l-v 

-  puiu^=  /..•  +  3^,  J  -  3  pko,, 


(3) 


3. 2  Droplet  Governing  Equations 

Considering  turbulence  diffusion  effects,  droplet  movement  equation  is  written  as: 


{  .  dU^ 
■  dV, 

m-i 


=  Co^^{V-iV, 


wl 


-  '-D  2  D  +  f»d  — 


,  •  dWj  Aj  juRef„r  TT’-  ^  • 

Coy^Cw- + 


(4) 


When  droplet  collide  with  wall,  the  droplet  rebounds  with  same  axial  velocity  and 
opposite  radial  velocity.  It's  trajectory  equation  is  : 


d^d  r  j  dr*^/  d^  ..tt 

=  d7  = 


dt 


dt 


(5) 


3. 3  Droplet  Vaporization  Model 
Droplet  vaporization  rate  can  be  calculated  by  following  formula: 


_  27rr^OgDgN  urn  I  r  1  -  y,.co. 

^  Inl  J 


here  A  =  1  + 


B  =  1  +  1) 

dt  A 


(6) 


The  droplet  radius  variation  can  be  written  as  the  following: 
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dt 


(7) 


my.s  D  dpi 
nD^pi  ^pi  dr 


Droplet  energy  equation  is: 


here  2^ 


dT/  _  3a^,N„jZi  -  Ti  _  AHz' 
dr  “  ICpinpi  e^.  -  1  Z2 


m. 
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Z,  =  AC,,  +  {A  -  DC 


pe 


Pb 

Pe.s 


(8) 


3. 4  Combustion  Model 

We  use  Spalding's  EBU  model  to  imitate  droplet  combustion,  assuming  that 
turbulent  combustion  zone  consist  of  large  quantity  of  flammable  gas  messes  and 
burned  gas  messes.  Chemical  reaction  take  place  at  interfaces  between  messes. 
Chemical  reaction  rate  depends  on  flammable  gas  messes  break  rate  which  is  directly 
proportional  to  turbulence  fluctuation  energy  decay  rate. 

X 

Rf,  If  =  -  min  [  I  jR/  I ,  I  I  ] 

.  Rf‘-  E/RT)  (5, 

here  Ck=1.07;  A  stands  for  reaction  constant; 

4. Results  and  Discussion 
4.1  Spray  Measurement 

The  spray  is  measured  downstream  at  a  section  of  110  mm  distance  from  nozzle 
exit.  At  different  radius  ,we  measure  droplet  density  Nd,  average  volume  diameter 
D30, average  area  diameter  D20  ,then  calculate  average  droplet  diameter  D32  by 
following  formula: 

rN/r)D3o(r)dr 

:  ^32  =  % -  (10) 

NAr)D2Q(r)dr- 

Varying  gas/liquid  injector  pressure  drop  ratio,  back  pressure  in  high  pressure  rig  and 
liquid  flowrate,  we  can  obtain  the  relation  between  D32  and  A  PG  /  A  PG  ,PB  ,m  L, 
as  seen  in  Fig.7-9.  The  Fig.  10  show  liquid  flowrate  distribution  v.s  radius.  It  can  be 
seen  that  it's  average  spray  droplet  diameter  is  smallest  for  centrifugal  coaxial  injector; 
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Average  diameter  Djj/  n  m  Average  diameter  D32/  u  m 


For  shear  coaxial  injector  and  bi-coaxial  injector,  their  spray  droplet  diameters  are 
closely;  For  bi-coaxial  injector,  when  it's  oxygen  post  recess  increases,  it's  average 
spray  droplet  diameter  will  decreases.  For  all  four  injector  configurations,  their 
average  spray  droplet  diameters  will  decrease  when  injector  pressure  drop  ratio 
increase,  and  their  average  spray  droplet  diameters  will  increase  when  back  pressure 
and  liquid  flowrate  increase.  Fig.  10  show  that  for  shear  coaxial  injector  and  bi¬ 
coaxial  injector,  their  flowrate  distribution  concentrate  on  spray  centre;  for  swirl 


coaxial  injector,  it's  flowrate  distribution 
ratio  dustribution. 


0.4  0.6  0.8  1  1.2  1.4  1.6 


aPg/ap, 

Fig.7  gas/liquid  injector  pressure 
drop  ratio  effects  on  D32 


0. 026  0. 028  0. 03  0. 032  0. 034  0. 036  0. 038 

Water  flow  /kg.s 

Fig.9  flowrate  effects  on  D32 


is  well  distributed,  that  benefits  mixture 


A  Pg/MPa 


Fig.8  back  pressure  effects 
on  D32 


R/mm 

Fig.  10  flowrate  distribution 


4.  2  Theoretical  Calculation 

Cartesian  coordinates  grid  system  are  used  for  theoretical  calculation.  We  will 
change  the  non-regular  physical  plane  into  regular  rectangular  plane  for  calculating. 
The  grids  distribution  is  made  exponentially.  In  order  to  ensure  calculation 
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convergence  and  accuracy,  we  make  denser  grids  at  near  injector  face  and  convergent 
section,  where  there  are  great  variation  in  chemical  reaction  and  flow  conditions. 


Fig.  11  Calculation  domain  and  grid  plane 


We  make  the  equations(l)  and  (2)  discretization.  The  liquid  phase  governing 
equations  (4)  and  (5)  are  coupled  with  gas  phase  governing  equations  through 
diffusion  velocities.  We  can  directly  integrate  the  liquid  phase  governing 
equations  ,while  assuming  that  gas  phase  field  is  known.  The  time  interval  required 
for  droplet  movement  from  position  k  to  k+1  is  decided  by  following  expression: 


min| 


Uj.t 


(11) 


here  U  denote  droplet  axial  and  radial  velocities  at  position  k.  The  SIMPLE 

was  used  to  solve  the  equation  system.  The  pressure  field  is  indirectly  given  by 
continuity  equation,  so  that  we  can  put  the  pressure  field  into  momentum  equation  to 
get  velocity  field.  As  soon  as  the  velocity  field  meets  continuity  requirements,  the 
calculation  will  goes  to  the  end,  otherwise  the  calculation  repeats  after  revising 
pressure  field  again.  In  calculation,  we  assume  that  all  flux  at  symmetric  axis  are 
equated  to  zero,  the  wall  is  adiabatic ,  fixed  and  can't  be  penetrated,  and  there  is  sonic 
velocity  at  prebumer  exit.  The  flow  conditions  are  greatly  coupled  with  throat  choke 
conditions.  The  inlet  parameters  are: 

liquid  flowrate  0.0885  kg/s;  liquid  oxygen  temperature  90K;  gaseous  hydrogen 
temperature  300K. 

Mixture  ratio  variation  can  be  realized  by  changing  gaseous  hydrogen  flowrate.  For 
centrifugal  injector,  liquid  oxygen  tangential  velocity  should  be  given.  For  bi-  coaxial 
injector,  gaseous  hydrogen  velocity  at  different  entrance  is  given  to  corresponding 
entrance  grids. 

Initial  spray  droplet  size  is  obtained  by  spray  testing  data.  Dividing  the  measured 
droplet  size  into  10  groups,  inputing  diameter  and  mass  percentage  of  each  group,  the 
calculation  begin  and  the  results  will  display  on  plane. 

For  four  kinds  of  prebumer  injector,  the  pictures  are  almost  same.  Now  we  take 
straight  coaxial  injector  as  example  to  explain  the  results.  Velocity  vector  field  is 
depicted  in  Fig.l2.  It  show  that  there  is  a  recirculation  zone  near  injector  face.Fig.13 
show  temperature  field,  it  imply  that  there  is  great  temperature  variation  in  the  range 
of  30-70mm  downstream  injector  face,  where  exists  strong  combustion  reaction. 
Liquid  oxygen  droplets  injecting  from  injector  move  forward  along  it's  trajectory. 
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owing  to  their  vaporization,  aroimd  them  exist  strong  combustion  zones.  The  oxygen 
vapor  immediately  reacts  with  gaseous  hydrogen  before  it's  sufEcient  diffusion,  so 
that  exists  higher  local  mixture  ratio.  The  highest  temperature  can  reach  about 
3000K.The  chemical  reaction  basically  stops  in  excess  of  70mm  downstream  injector 
face.  Gas  temperature  gradually  become  closely  and  the  final  temperature  is  about 
lOOOK.  At  different  axial  distance,  the  temperature  variation  relations  with  radius  are 
depicted  in  Fig.l4.  The  temperature  near  injector  face  (X<26.5mm)  is  under  800K, 
that  imply  combustion  don't  take  place.  The  combustion  begins  until  the  distance  is  in 
excess  of  26.5mm  and  stops  about  at  X=78.9.The  low  temperature  zones  near  injector 
face  are  produced  by  recirculation.  After  combustion  begin,  exists  initial  temperature 
peak  in  radial  direction  about  at  Y=5mm,  then  the  peak  move  gradually  toward  the 
centre.  In  excess  of  X=79.9mm,  gas  temperatures  are  almost  equal.  Radial  distribution 
of  axial  velocities  is  shown  in  Fig.  15. It  imply  that  recirculation  zones  exist  at  0-43mm 
axial  distance  and  the  flow  finally  become  stable  accompanying  with  the  end  of 
combustion.  Fig.  16  show  that  radial  velocities  in  combustion  zone  are  relatively 
higher,  owing  to  turbulence  mixture  and  then  they  decrease  gradually.  The  radial 
velocities  are  almost  equal  to  zero  before  the  convergent  section(X=173.9).  That 
means  diffusion  and  mixture  are  very  weak,  the  flow  is  in  equilibrium. 


Fig.  1 2  velocity  vector  field 
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Fig.  13  isothermal  field 
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Fig.  1 4  gas  temperature  radial  distribution  Fig.  1 5  radial  distribution  of  axial 
at  different  sections  velocity  at  different 

sections 


H  (m) 

Fig.  1 6  radial  distribution  of  radial 
velocity  at  different  sections 


4.3  Fire  Testing  Results 

Six  fire  testings  for  prebumer  have  been  finished.  The  testing  measured  data  are 
summarized  in  table  1 . 


Table  1  testing  data 


mode 

YROl 

YR02 

YR03 

YR04 

YR06i 

YR06,i 

prebumer 

pressure 

Pf,  MPa 

7.400 

7.585 

7.489 

7.080 

7,732 

6.848 

0.376 

0.375 

0.369 

0.367 

0.322 

275 


H2  flow 

Gr,  kg/s 

0.438 

0.448 

0.461 

0.461 

n 

0.470 

mixture  ratio,  r 

0.845 

0.838 

0.813 

0.799 

n 

0.673 

characteristic 
velocity  C,  m/s 

2408 

2269 

2308 

2272 

gas  temperature 
Tfl  ,K 

n 

/ 

/ 

1003 

900 

gas  temperature 
Tf2,K 

■a 

1091 

1058 

953 

866 

gas  temperature 
Tf3,K 

1134 

1103 

1056 

1054 

989 

896 

theoretical 

temperature 

Tfn,,K 

1118 

1108 

1082 

1072 

1044 

954 

biggest 
temperature 
difference  A  Tf,K 

28.98 

20.09 

35.6 

4.2 

50.23 

34.73 

A.  Combustion  Efficiency.  It  can  be  calculated  by  the  following; 

nc=C*„/C*,  C*.=PfAt/G,  (12) 

here  is  experimental  characteristic  velocity; 

C*,h  is  theoretical  characteristic  velocity,  corrected  by  supply 
gas  hydrogen  temperature 
At  is  throat  area  of  exhaust  duct 

The  other  measured  data  are  all  close  to  100%,  this  is  owing  to  having  bigger 
hardware  volume.  Theoretical  calculation  show  that  the  propellant  combustion  goes 
on  to  the  end  in  100mm  downstream  injector  face,  combustion  efficiencies  greater 
than  100%  result  in  the  measurement  error.  Comparing  four  prebumer  configurations, 
we  found  that  option  D  has  higher  combustion  efficiency  and  lower  for  B. 

B.  Gas  Temperature  Distribution 

Gas  temperatures  are  measured  at  different  locations  in  a  same  exhaust  duct 
section.  The  increase  rate  of  measured  temperature  with  mixture  ratio  is  more  rapid 
than  that  of  theoretical  calculation.  Comparison  of  the  temperature  difference  show 
that  the  gas  temperature  is  well  distributed  for  D  ,  A  Tf=20.09,  and  it's  temperature 
distribution  is  poor  for  B,  A  Tf=50.23. 

C.  Vibration 

Vibration  gauges  were  mounted  at  the  head  cover  of  injector  elements  in  axial, 
radial  and  tangential  directions.  The  measured  vibration  data  are  given  in  table  2,  the 
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peak  frequency  and  it's  amplitude  are  given  in  table  3.  Comparing  their  vibration 
data  of  different  injector  configurations  ,  there  is  lower  vibration  for  D,  it's  complex 
vibration  is  325.4  m/s^2,  and  it's  component  vibration  is  36.8  m/s^2  (3025  Hz). 
The  option  B  has  higher  vibration,  it's  complex  vibration  is  686  m/s^,  component 
vibration  is  196.4(5575  Hz).  The  measured  gas  pressure  fluctuation  show  that  there 
is  no  combustion  instability,  the  peak  to  peak  pressure  fluctuation  is  less  than  4%  of 
rated  chamber  pressure. 


Table  2  mechanical  vibration 


parameter  name 

YROl 

YR02 

YR03 

YR04 

YR06 

option 

C 

D 

A 

B 

B 

axial  acceleration 
a,,  m/s^2 

306.3 

261.7 

514.5 

196 

356.7 

radial  acceleration 
a2,  m  ls'^2 

414.5 

279.6 

676.2 

686 

541 

tangential  acceleration 
a3,m/s^2 

22,2 

325 

/ 

174 
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Table  3  mechanical  vibration 


parameter  name 

YROl 

YR02 

YR03 

YR04 

YR06 

axial  frequency 

Hz 

msm 

5650 

axial  amplitude 

m/s''''2 

TX6 

32 

Hz 

5675 

5575 

radial  amplitude 

m/s^2 

mm 

tangential  frequency 

nai 

mam 

iggn 

5575 

tangential  amplitude 

15.1 

mm 

D.  Ignition  Performance 

We  employ  pyrotechnic  ignitor  AN8100-00B  which  has  150mm  flame  length  and 
about  four  second  burning  duration.  Pressure  peaks  occur  at  about  0.6  s.  After  liquid 
oxygen  and  gaseous  hydrogen  inject  into  prebumer,  the  chamber  pressure  rapidly 
increases  up  to  peak  value,  there  is  a  time  interval  between  hydrogen  entering  time 
and  ignition  time,  which  were  given  in  table  4.  The  ignition  performance  is  closely 
related  to  the  time  interval,  and  the  shorter  interval  means  well  ignition  performance. 
The  data  show  that  option  A  has  shorter  time  interval  and  it  is  longest  for  option  B. 
this  may  be  related  to  mixture  ratio  distribution  (Fig.lO).  For  swirl  coaxial  injector, 
oxygen  injection  flow  has  strong  transversal  penetrability  into  hydrogen,  so  that  it 
can  improve  mixture  ratio  distribution  and  ignition  performance. 
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Table  4  Ignition  delay 


parameters 

name 

YROl 

YR02 

•  YR03 

YR04 

YR06 

option 

C 

D 

A 

B 

B 

hydrogen 
enter  time,  s 

0.4 

0.4 

0.7 

0.4 

0.4 

ignition 
time ,  s 

1.6 

1.6 

1.6 

1.8 

1.8 

time 

interval,  s 

1.2 

1.2 

0.9 

1.4 

1.4 

S.ConcIuding  Remarks 

A.  Gas/liquid  coaxial  injector  spray  droplet  diameters  decrease  with  the  increase  of 
gas/liquid  injector  pressure  drop  and  increase  with  the  increase  of  liquid  flowrate  and 
back  pressure. 

B.  Average  spray  droplet  diameter  for  swirl  coaxial  injector  is  smaller  than  that  of 
shear  coaxial  injector.  It  has  good  transversal  penetrability  and  ignition  performance. 

C.  The  bi-coaxial  injector  (option  D)  has  the  advantage  of  high  combustion 
efficiency  and  stability,  uniform  temperature  distribution  but  it  needs  to  further 
improve  it's  ignition  performance. 

D.  For  gas/liquid  coaxial  injectors,  there  are  recirculation  zones  in  the  range  of  0- 
43mm  downstream  injector  face,  and  the  combustion  occurs  in  range  of  30-70mm 
from  injector  face. 


Reference 

1.  Sankar  S  .V,  Zhu  J.Y,  Bachalo  W.D.  Experimental  studies  on  the  behavior  of 
coaxial  rocket  injector  spray.  Areome  Trie,  Sunnyvale,  CA,  U.S.A.  1994 

2.  Lefebvre  A.H.  Atomization  and  spray  Hemisphere  Pub.  XO,  1989 

3.  Yost  M.C.  Prebumer  of  staged  combustion  rocket  engine,  NASA  CR-1 35356 

4  Z.G.Wang  Combustion  Processes  imitation  calculation  for  liquid  rocket  engine. 
National  Defence  Science  and  Technology  University,  China,  1995. 
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Calculation  of  MASCOTTE 
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Equations  for  the  Lagrangian  tracing  of  particies 
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Equations  solved  by  MSD 
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MASCOTTE  simulation 


MASCOTTE  simulation  at  1  Bar  (A-l) 

GH2  mass  flow  rate  :  15  g/s 
LOX  mass  flow  rate  :  50  g/s 
k-e  model 


M  ASCOTTE  simulation  at  1  Bar  (A-1 ) 


M ASCOTTE  simulation  at  1  Bar  (A-1 ) 


M ASCOTTE  simulation  at  1  Bar  (A-1 ) 
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M  ASCOTTE  simulation  at  1  Bar  (A-1 ) 
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MASCOTTE  simulation  at  10  bars  (A-10) 

GH2  mass  flow  rate  :  23.7  g/s 
LOx  mass  flow  rate  :  50  g/s 
EBUP  =  2500 
k-e  model 
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Conservation  equations  for  the  continuous  phase 
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Work  to  be  done 
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RESEARCH  STATUS  AND  PERSPECTIVES  IN 

LIQUID  ROCKET  COMBUSTION  CHAMBER  FLOW  DYNAMICS 


Basic  Research  in  Liquid  Rocket 
Combustion  at  the  Air  Force 
Research  Laboratory 


Doug  Talley 
Propulsion  Sciences 
and  Advanced 
Concepts  Division 

27-28  May  1999 


Required  Injector  Characteristics 


Complete  combustion  in  the 
shortest  possible  length 

-  Main  injectors: 
performance  vs  weight 
tradeoffs 

-  Preburners/GG’s: 
downstream  component 
interactions,  eg,  turbine 
blades,  etc 

Acoustically  stable 

-  Chamber  modes 

-  Feed  system  coupling 
Chamber/wall  compatibility 

-  Heat  transfer/cooling 

-  Oxygen  blanching 


•  Minimize  pressure  drop 

•  Throttling 

•  Ignitable;  minimum  ignition 
transients 

•  Cost,  weight 

•  The  “ilities:” 

-  Reliability 

-  Maintainability 

-  Manufacturability 

-  Durability 

-  Operability 


Subscale  hot  fire  facility 

Fuel 

H2(g).  CH4(g) 

Oxidizer 

02(g) 

Purge  gas 

N2(g).  He(g) 

H2  mass  flow  rate 

.15  ibm/s  (.07  Kg/s) 

CH4  mass  flow  rate  .25  Ibm/s  {.  1 1  Kg/s) 

02  mass  flow  rate 

l.O  Ibm/s  (.45  Kg/s) 

N2  mass  flow  rate 

.5  Ibm/s  (.23  Kg/s) 

Water  flow  rate 

161bnVs(7  Kg/s) 

Max.  system  press 

2640  psi.  (179  atm) 

128  ch. 

200  kbs  scanning  .A/D 

16  ch,  2  MHz  per  ch  A/D,  independently  controlled 

Central  laser/optics  room 

.  2000  psi  liquid  hydrocarbon  capability 

20QQ 

LOX  capability 

i} 


4  Ihp 


Gas/gas  hardware 


311 


312 


314 


^  Overall  Objectives 

Determine  the  impact  of  thermodynamic  state 
relative  to  the  critical  point  on  the  transport  and 
mixing  of  droplets,  jets,  and  sprays 


Issues 

•  Vanishing  surface  tension  and  heat  of  vaporization 

•  Equivalent  gas/liquid  densities 

•  Enhansed  gas/liquid  solubility 

•  Mixing  induced  critical  point  variations 


\v 
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f  Low  Reynolds  Number  N2  Jets 


Woodward  and  Talley,  AIAA  96-046S 


y; 


N2  N2 

Pr=0.83  Pr=1.03 


N2  N2/He.  MR=3.9 

Pr=2.03  Pr=2.03 


Raman  images  of  LN2  jet  at  69  atm  into  (a) 
gaseous  nitrogen  and  (b)  helium 


I  T 


t  «ip  t  SIlL't!! 


-4 
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Characteristic  Times 


Characteristic  bulge  formation  time  {x^  at  the  jet 
interface  (Tseng  et  al.):  (p|L^/o)''*  ;  P|  ,  L,  a  are  lic^uid 
density,  characteristic  dimension  of  turbulent  eddy,  and 
surface  tension,  respectively. 

Characteristic  time  for  gasification  (x^)  (D-sc^uare  law): 

D  and  K  are  drop  diameter  and  vaporization 
constant. 

A  Hypothesis:  If  these  two  characteristic  times  (calculated 
for  appropriate  length  scales)  are  comparable  then  an 
interface  bulge  may  not  be  separated  as  an  unattached 
entity  (onset  of  the  gas-jet  behavior  at  supercritical 
condition) 


•  Theoretical  isothermal  lic^uid  spray  growth  rate  (G^)  based 
on  0rr-6ommerfeld  equation  and  stability  analysis  to  find 
the  wavelength  of  the  most  unstable  interface  wave: 
0=O.27[O  +  (p/p,)°-5] 

•  Papamoschou/f^ashko  theory  for  incompressible  variable- 
density  gaseous  mixing  layer/jet: 

ep,,=  o.i7  [1  +  (p/Pi)^'^] 

•  Dimotakis  theory  for  incompressible  variable-density 
gaseous  mixing  layer/jet: 

GpS  0.212  [0.59 +  (p,/p,)0-5] 

•  ALL  HAVE  THE  SQUARE  ROOT  OF  DENSITY  RATIO  AND 


THE  SAME  EQUATION  FORMAT 


f  Correlation 

of  the  information  of  the  previous  sli^e  the  following  correlation 
is  proposed  for  both  sub-  and  supercritical  measured  growth  rates: 


Qa,=  0.27  [  (T,/(V  T,))  +  (p/p,)°'®] 


N2  into  N2 


GchS  0.27  [  (V(T,+  Tj)  +  {(P/P,)-0.6(P/P,)} 


N2  into  He 


For  isothermal  \\<\\S\d  case:  >>  and  x^  — ^  It  then  collapses  to 

the  isothermal  spray  case. 

For  subcritical  the  (Xi/(Xt,+  X^))  is  calculated  until  it  reaches  0.5.  After 
that  it  is  maintained  constant  at  0.5  for  supercritical  jet.  The 
transition  point  is  found  to  be  approximately  when  (Xb/(Xi,+  X^))  =  0.5. 
Only  a  horizontal  axis  variable  transformation  is  needed  to  fit  with  the 
N^-into-He  experimental  growth  rate  data. 


Empirical  Correlation  Results 


O  N2  J«t  into  N2  L/D>200  (*) 

♦  Cold  Ho  Jot  into  N2;  L/D>200  (*) 

■  Cold  N2  jot  Into  Ho;  L/Da200  (*) 

M -  Brown  &  Bochko  (Ho/N2)lncomprocsiblo  variablo-donctty  Mixing  Layor 
A  02  jot  Into  N2:  L/D>200  (*) 

- Incomproaaiblo  Varlablo-Donsity  (PapamocchouiRochko  thoory) 

- Dimotakis  (thoory) 

■*~~Propoaod  Modol  (Chohroudl  ot  al.) 

■  N2  jot  into  Ha;  L/0>200  (*} 

X  N2  into  (0.1  CO+0.9N2):  L/D-200  (•) 

▲  N2  Into  (0.5CO.^0.SN2);  L/D»200  (*) 

^  N2  Into  (0.5HO-I.0.5N2)  (*) 

Proposod  Modol  (Chohroudl  ot  al) 


N2  into  He 


N2  into  N2 


■  Correlation 


Chamber/Injectant  Density  Ratio 
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Fractal  Dimension  vs  Reduced  Pressure 

Chehroiidi  et.  al.,  AIAA  99-2489 


f  Summary  and  Conclusions 


•  Structural  differences  in  cryogenic  jets  have  been  observed  below 
and  above  the  thermodynamic  critical  point. 

•  Liquid-Jet  like  appearance  occurs  up  to  near  the  critical  point, 
similar  to  second  wind-induced  liquid  jet  breakup  regime. 

•  Gas-jet  like  appearance  occurs  above  the  critical  point.  No  drops 
are  observed. 

-  Supercritical  spreading  rate  measurements  agree  quantitatively  with 
incompressible  variable  density  mixing  layer  experiments  and  theory. 

-  Supercritical  fractal  dimensions  agree  quantitatively  with  gas  jet 
measurements. 

•  New  and  existing  mixing  layer  growth  rate  experiments  and  theory 
have  for  the  first  time  been  consolidated  into  a  single  plot  as  a 
function  of  density  ratio,  where  the  density  ratio  spans  three  orders 
of  magnitude. 

•  A  physical  mechanism  and  correlation  have  been  proposed  to 
describe  the  transition  from  spray  to  gas  jet  behavior. 


Opportunities  for  Future  Research 

1990’s  at  AFRL 

-  High  pressures,  “steady.” 

-  Gas/gas  injection 

Future:  1990’s  topics  remain  relevant.  Add: 

-  Injector/chamber  interactions 

-  Transients 

•  Organized 

•  Un-organized 

-  Revolutionary  cycles 

•  Pulsed  detonation  propulsion 

•  Combined  cyole 
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THE  STRUCTURE  OF  CRYOGENIC  JET 
FLAMES  INFERRED  FROM 
EXPERIMENTS 


A.  Tripathi,*,  M.  Juniper, ^  P.  Scouflaire,^ 
D.  Durox,^,  C.  Rolon,^,  S.  Candel,^ 

EM2C  Laboratory,  CNRS,  Ecole  Centrale  Paris 
F-92295  Chatenay-Malabry,  France 

L.  Vingert,®  and  M.  Habiballah,"" 

ONERA,  91120  Palaiseau,  France 


Abstract 


This  article  reports  new  experimental  results  on  cryogenic  jet  flames  formed 
by  a  coaxial  injector.  This  element  fed  with  liquid  oxygen  and  gaseous  hydro¬ 
gen  is  placed  in  a  square  combustion  chamber  equipped  with  quartz  windows. 
The  flame  is  studied  with  a  set  of  optical  methods  :  light  emission  from  OH 
radicals  and  from  H2O  molecules,  laser  induced  fluorescence  of  OH,  elastic 
scattering  from  the  liquid  oxygen  jet,  shadowgraphy.  These  techniques  are 
used  to  examine  flame  spreading  near  the  injection  plane.  The  instanta¬ 
neous  flame  location  was  previously  described  with  simultaneous  light  elastic 
scattering  images  and  laser  induced  fluorescence  of  OH.  More  recently,  the 
mean  flame  structure  was  extracted  from  average  emission  images  treated 
with  Abel’s  transform.  This  numerical  transformation  provides  the  local 
volumetric  light  emission  of  the  species  being  imaged.  Previous  studies  con¬ 
cerned  radiation  from  OH  radicals.  New  measurements  have  been  carried  out 
by  imaging  a  band  of  wavelengths  which  characterizes  H2O  radiation.  The 
transformed  images  of  OH  yield  a  structure  which  is  indicative  of  the  loca¬ 
tion  of  reaction  while  the  transformed  images  of  H2O  radiation  portray  the 
distribution  of  products  or  temperature  in  the  flame.  This  combined  imag¬ 
ing  provides  further  information  on  the  flame  structure.  Additional  clues  on 
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the  flow  near  the  injector  are  provided  by  shadowgraph  images.  System¬ 
atic  experiments  have  been  carried  out  by  varying  the  operating  parameters 
and  the  cliamber  pressure.  Effects  of  LOX  tube  recess  have  also  been  ex¬ 
amined  at  various  pressures  and  for  a  .set  of  operating  conditions.  The  data 
provide  a  unified  picture  of  the  flame  in  the  vicinity  of  the  injection  plane. 
Results  may  be  used  to  describe  effects  of  pressure  and  analyze  the  modifica¬ 
tion  in  the  flame  pattern  induced  by  recess.  The  experimental  investigation 
serves  as  a  basis  for  modeling  efforts  described  in  a  companion  pap&r  and  in 
a  previous  article  cited  in  reference. 

This  work  was  carried-out  with  support  from  CNES,  SEP,  OhiERA  and 
CNRS. 


Figure  1:  Cryogenic  combustion  facility 


1  Experimental  Facility 

The  experimental  facility,  called  Mascotte,  shown  in  Figure  1,  is  situated 
at  ONERA,  Palaiseau.  It  consists  of  a  coaxial  injector  in  which  an  annulus 
of  hydrogen  at  high  velocity  is  injected  around  a  jet  of  liquid  oxygen  at  low 
velocity.  The  dimensions  of  the  injector  are  in  the  same  range  as  those  used  in 
the  Vuicain  rocket  motor  of  Ariane  5.  A  pressure  of  60  bar  can  be  maintained 
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Figure  2:  Low  pressure  injector  High  pressure  injector 

within  the  chamber,  which  is  square.  The  first  few  centimetres  of  the  flame 
can  be  examined  through  quartz  windows  on  each  side  of  the  chamber. 

Cold  flow  studies  have  shown  the  importance  of  the  momentum  ratio, 
defined  as  J  =  in  determining  jet  break-up  mechanism.  This 

parameter  is  kept  constant  on  changing  pressure,  which  requires  a  reduction 
of  the  hydrogen  flow  cross-section  at  high  pressure.  The  injector  geometries 
are  shown  in  Figure  2. 

Our  laboratory  has  carried  out  five  series  of  experiments  on  Mascotte,  the 
details  of  which  are  summarised  in  Table  (1).  Versions  VOl  and  V02  of  the 
combustion  facility  permitted  pressures  up  to  10  bar.  At  these  pressures,  the 
flame  position  was  determined  using  planar  laser-induced  fluorescence.  The 
liquid  surface  position  was  found  simultaneously  by  observing  the  elastic 
scattering  of  the  laser  sheet.  It  is  known  that  LIF  and  elastic  scattering 
are  less  effective  at  high  pressure,  as  well  as  being  technically  complicated. 
It  is  more  suitable  to  use  visualisation  techniques  based  on  emission  and 
shadowgraphy.  These  techniques  are  simpler,  but  quite  effective. 


2  Experiments  at  1,  5  and  10  bar 

A  diagram  of  the  apparatus  used  for  planar  laser-induced  fluorescence  is 
shown  in  Figure  3.  The  laser  sheet  excites  OH  radicals  in  the  Qi{6)  band, 
which  decay  emitting  a  photon  at  a  different  wavelength.  The  emission  in¬ 
tensity  is  found  to  be  almost  independent  of  temperature  in  a  range  between 
1500  and  3000K.  The  signal  then  mainly  depends  on  the  concentration  of 
hydroxyl  (OH)  radicals.  These  radicals  are  an  intermediate  in  the  global 
reaction  of  H2  with  O2  and  regions  of  high  level  of  OH  mark  the  flame  front. 
Filters  exclude  all  but  the  fluorescence  wavelength  and  thus  the  flame  po¬ 
sition  is  recorded  by  an  intensified  CCD  camera.  Simultanously,  a  second 
CCD  camera  records  the  scattering  of  the  laser,  sheet  by  the  liquid  surface. 
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Version 
of  rig 

Dates 

Pressure 

Visualisation  technique 

3  series 
of  tests 

1,5,10  bar 

Laser-induced  fluorescence 
Elastic  scattering 

OH  emission 

V03 

June  1998 
March  1999 

10,60  bar 

OH  emission 

H20  emission 

Spectroscopy 

Shadowgraphy 

Table  1:  Experiments  on  MASCOTTE  by  Laboratoire  EM2C 


Figure  3:  Set-up  for  simultaneous  Planar  Laser-induced  Fluorescence  and 
Elastic  Scattering  imaging 
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Figure  4:  Schematic  diagram  of  the  flame  between  1  and  10  bar,  below  the 
critical  pressure  of  oxygen 

This  gives  the  position  of  the  jet  and  droplets. 

The  flame  front  ^  in  red,  is  not  seen  in  the  top  half  of  the  images  because 
very  little  of  the  laser  sheet  passes  undisturbed  through  the  jet.  The  bottom 
half  shows  that  the  flame  front  is  highly  wrinkled.  Scattering  from  the  lic[uid 
surface  is  shown  in  blue.  Further  from  the  injector,  the  intensity  is  lower, 
which  is  due  to  the  fact  that  the  droplets  are  smaller.  Also  the  thicker  droplet 
cloud  masks  the  sheet  of  droplets  which  are  illuminated  by  the  laser.  The 
scattered  signal  is  outside  the  dynamic  range  of  the  CCD  camera. 

Some  important  points  can  be  made  about  the  flame  shape.  The  flame 
starts  at  the  oxygen  injector  lip  and  then  burns  extremely  close  to  the  liquid 
jet.  The  liquid  jet  breaks  up  into  a  dense  cloud  of  droplets.  The  inter  droplet 
distance  divided  by  the  droplet  diameter  is  estimated  to  be  around  b/d  —  6. 
Because  the  cloud  is  so  dense,  external  group  combustion  takes  place  and 
the  flame  surrounds  but  does  not  penetrate  the  cloud.  A  schematic  diagram 
of  the  flame  shape  between  1  and  10  bar  is  shown  in  Figure  4. 


3  Application  of  The  Abel  Inversion 

Given  the  reduced  signal  level  of  PLIF  at  high  pressures,  other  techniques 
were  considered.  One  of  these  is  the  Abel  inversion  of  natural  emission 
images.  The  spectrum  of  the  flame  shows  that  is  it  possible  to  isolate  the 
natural  emission  of  OH*  and  H2O*  radicals,  as  shown  in  Figure  5. 

The  Abel  inversion  is  frequently  used  to  obtain  the  radial  intensity  dis¬ 
tribution  of  an  axisymmetric  object  from  a  projection  on  a  plane  placed  on 

^Colour  images  can  be  found  in  Candel  et  al:  '’Experimental  investigation  of  shear 
coaxial  cryogenic  jet  flames,”  Journal  of  Propulsion  and  Power,  Vol.  14  ,  1998,  pp.  826- 
834 
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Figure  5:  Spectrum  of  the  flame  at  10  bar  a)  300  <  A  <  400nm  ;  UV  region 
b)  500  <  A  <  700nm;  green-red  region 
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Figure  6:  Principle  of  the  Abel  inversion,  the  radial  structure  is  deduced 
from  the  projected  distribution 
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Figure  7:  Instantaneous  image  of  the  OH*  emission  of  a  turbulent  flame,  10 
bar,  J  =  14.5 


Figure  8:  Averaged  images  of  a  turbulent  flame.  This  gives  the  OH*  Intensity 
X  Residence  fraction  summed  along  the  line  of  sight.  10  bar,  J  =  14.5 

the  side.  The  principle  is  shown  in  Figure  6.  The  axisymmetric  distribu¬ 
tion  is  assumed  to  be  made  up  of  a  series  of  co-centric  rings.  To  calculate 
the  contribution  of  an  inside  ring,  the  contribution  of  the  rings  outside  it  is 
subtracted  from  the  cumulative  image  along  the  line  of  sight.  Working  from 
the  outside  to  the  inside,  this  enables  the  contribution  from  each  ring  to  be 
determined.  Around  the  axis,  the  process  is  particularly  sensitive  to  noise. 

Care  is  required  w'hen  applying  the  Abel  transform  to  a  non-axisymmetric 
turbulent  flame,  as  indicated  in  Figure  7.  Averaging  over  several  images  gives 
an  axisymmetric  field:  intensity  x  residence  fraction,  which  is  summed  over 
the  line  of  sight,  as  in  Figure  8.  The  residence  fraction  corresponds  to  the 
fraction  of  the  time  that  the  flame  spends  at  a  particular  point,  on  average. 
The  Abel  inversion  can  be  applied  to  this  field,  giving  a  slice  through  the 
injector  plane  of  the  intensity  of  OH*  radical  natural  emission  x  residence 
fraction.  Figure  9.  This  can  be  compared  with  the  average  of  the  PLIF 


Figure  9:  The  Abel  inversion  of  the  averaged  images.  This  gives  a  slice  of 
the  OH*  Intensity  x  the  Residence  fraction.  10  bar,  J  =  14.5 

images,  which  gives  a  slice  of  the  intensity  of  OH  radical  stimulated  emission 
X  residence  fraction,  as  in  Figure  10.  This  indicates  a  slightly  wider  high 
intensity  region  than  that  found  in  the  Abel-inverted  images.  However,  this 
is  to  be  expected  because  the  PLIF  images  show  all  OH  radicals,  wheras  the 
Abel-inverted  images  show  only  the  excited  OH*  radicals.  The  OH*  radical 
has  a  short  lifetime  compared  with  characteristic  flow  times  and  will  not  be 
seen  away  from  the  flame  front.  However,  it  is  likely  that  a  small  quantity  of 
OH  radicals  will  exist  in  the  high  temperature  regions  around  the  flame.  The 
corresponding  intensity  of  these  radicals  on  the  PLIF  images  will  be  small 
but  the  residence  fraction  will  be  high.  This  accounts  for  the  thick  flame 
zone  on  the  averaged  PLIF  images. 

To  compare  the  averaged  images  more  exactly  it  is  necessary  to  make 
the  reasonable  assumption  that  the  regions  where  OH  concentration  is  high 
correspond  to  the  regions  of  OH*.  Thus  when  the  low  intensity  regions  of 
the  PLIF  images  are  removed  and  the  images  averaged,  the  result  should  be 
similar  to  the  intensity  of  OH*  x  Residence  fraction.  The  comparison  turns 
out  to  be  excellent.  Figure  11  shows  the  positions  of  maximum  intensity  for 
the  Abel  inverted  images  and  for  the  average  of  the  PLIF  images. 

Whenever  the  Abel  inversion  is  used,  it  should  be  remembered  that  the 
slices  are  of  emission  intensity  x  residence  fraction,  rather  than  just  intensity. 
The  highest  values  are  usually  seen  near  the  injector,  where  the  residence 
fraction  is  high.  The  Abel-inverted  images  give  a  good  idea  of  the  change 
in  the  average  flame  shape  when  injection  conditions  are  altered.  It  should 
also  be  noted  that  the  Abel  inversion  assumes  that  no  light  is  absorbed  or 
deflected  by  the  flame. 

The  Abel-inverted  images  of  OH*  and  H2O*  I'adical  emission  are  shown 
in  Figures  12  and  13.  To  a  first  approximation,  the  intensity  of  OH*  emission 
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Figure  10:  Averaged  PLIF  images.  This  gives  a  slice  of  the  OH  Intensity  x 
the  Residence  fraction 


Figure  11:  Position  of  maximum  intensity  of  the  Abel-inverted  emission  im¬ 
ages  compared  with  that  of  the  average  PLIF  images 
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is  a  fiinctioTi  of  the  reaction  rate  and  the  intensity  of  the  thermally  excited 
reaction  product  Il-iO*  is  a  function  of  temperature.  From  these  im^rges, 
taken  with  an  oxygen  tube  recess  of  1  xdiox,  it  can  be  seen  that  tlu'  max¬ 
imum  ll-iO’'  emission  is  somewhat  further  downstream  than  the  maximum 
Oir  emissiom  as  would  be  expected. 


Figure  12:  OH*  Intensity  x  Residence  fraction  at  10  bar  in  a  slice  through 
the  axial  plane 


Figure  13:  H2O*  Intensity  x  Residence  fraction  at  10  bar  in  a  slice  through 
the  axial  plane 


4  Supercritical  effects 

The  critical  point  of  oxygen  is  50.8  bar,  155K.  The  Pressure- Temperature 
phase  equilibrium  diagram  is  shown  in  Figure  14.  Conditions  achieved  in 
Mascotte  version  VOl,  V03  and  in  the  Vulcain  motor  are  also  shown  in  this 
figure.  The  important  point  is  that  the  oxygen  will  have  no  definite  phase 
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Figure  14:  Pressure-Temperature  diagram  for  oxygen,  with  the  operating 
points  of  Mascotte  versions  VOl,  V03  and  the  Vulcain  rocket  motor 

boundary  in  V03  and  in  Vulcain,  unlike  the  low  pressure  VOl.  This  is  likely 
to  lead  to  differences  in  jet  and  droplet  break-up,  since  surface  tension  is  a 
factor  at  low'  pressure.  Above  the  critical  pressure,  two  gases  of  very  different 
density  mix. 

Studies  by  Delplanque  and  Sirignano  indicate  that  the  oxygen  jet  will  have 
a  dense,  licjuid-like  core,  where  the  temperature  is  low.  Tow'ards  the  edge  of 
the  jet,  where  the  temperature  is  higher,  the  oxygen  w411  be  supercritical. 
There  wdll  be  very  high  temperature  and  density  gradients  in  this  region, 
but  no  discontinuous  surface.  Around  the  critical  point,  small  temperature 
variations  lead  to  large  density  variations.  The  fluid  becomes  almost  opaque 
to  incident  light. 

5  Experiments  at  60  bar 

The  apparatus  for  flame  and  jet  visualisation  at  high  pressure  is  shown  in 
Figure  15.  The  technicpies  used  are  listed  below': 

•  Spectroscopy:  The  spectrum  of  the  flame  at  60  bar  was  obtained  over 
the  course  of  the  campaign.  It  showed  that  011"“  radical  emission  is 
well  isolated  by  the  filters. 

•  Images  of  OH*  radical  emission 
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Figure  15:  Experimental  set-up  used  in  the  high  pressure  tests 


Figure  16:  OH*  radical  emission  at  low  pressure  {P  ^  30bar) 

•  Shadowgraphy:  In  order  to  eliminate  as  much  of  the  flame’s  natural 
emission  as  possible,  the  ICCD  camera  had  a  very  short  exposure  time 
(300ns)  and  was  synchronised  with  the  spark  lamp.  In  addition,  after 
passing  through  the  chamber,  the  light  was  focussed  through  an  iris. 
This  cut  out  a  large  proportion  of  the  flame’s  natural  emission.  Finally, 
the  light  was  filtered  to  remove  as  much  natural  emission  as  possible, 
while  passing  the  wavelength  in  which  the  spark  was  most  intense, 
between  420  and  500  nm. 

Some  close-up  images  of  the  near-injector  region  were  also  taken. 

The  far-field  results  are  shown  in  Figures  16  to  19.  During  each  run 
the  pressure  took  30  seconds  to  build  up  to  60  bar  and  then  was  stable  for  7 
seconds.  Thus  images  at  low  pressure  were  obtained,  but  it  should  be  stressed 
that  the  conditions  were  not  steady  and  do  not  correspond  to  those  of  v’ersion 
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Figure  17:  OH*  radical  emission  at  60  bar 


Figure  18:  Shadowgraphy  at  low  pressure  (P  «  30bar) 


Figure  19:  Shadowgraphy  at  60  bar 
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Figure  20:  Emission  above  400nm  at  low  pressure,  near  injector 

\'01  and  V02  of  Mascotte.  It  is  clear  from  the  shadowgraphy  images  that 
the  shape  of  the  oxygen  jet  is  different  at  high  pressure.  Droplets  are  no 
longer  visible  and  it  is  more  clumpy  and  stringy.  Emission  is  more  intense  at 
high  pressure  and,  compared  with  the  stable  images  at  1,5  and  10  bar  from 
version  VOl  and  V02,  the  flame  expands  less  rapidly. 

The  near-field  emission  images  are  shown  in  Figures  20  and  21.  The 
flame  is  stabilised  just  behind  the  oxygen  injector  lip.  This  effect  has  also 
been  noted  by  Mayer  and  Tamura  on  the  cryogenic  combustion  facility  in 
Japan,  on  which  the  diameter  of  the  oxygen  jet  is  smaller  and  the  oxygen 
lip  proportionally  much  larger,  figure  (22).  The  near-field  shadowgraphy 
images  require  some  explanation.  Typical  samples  are  shown  in  Figures  23 
and  24.  Where  the  light  from  the  spark  lamp  has  been  deflected  only  slightly, 
a  shadowgraphy  effect  can  be  seen,  characterised  by  a  light  and  a  dark  band. 
Where  the  light  has  been  deflected  more,  it  is  not  captured  by  the  lens  and 
a  schlieren  effect  can  be  seen.  This  is  simply  a  dark  region.  A  shadowgraphy 
effect  can  be  seen  at  the  outside  edge  of  the  hydrogen  jet  and  in  some  parts 
of  the  gaseous  phase  in  the  low  pressure  image.  The  60  bar  image  is  almost 
entirely  determined  by  schlieren  effects.  Light  is  deflected  by  refractive  index 
gradients  perpendicular  to  the  ray  direction.  Refractive  index  gradients  will 
be  large  at  the  edge  of  the  oxygen  jet  and  at  the  edges  of  the  flame.  If  the 
flame  is  very  close  to  the  oxygen  jet,  the  two  will  be  indistinguishable. 

The  images  show  the  expansion  of  the  hydrogen  annulus,  the  wake  of  the 
oxygen  tube  lip,  the  position  of  the  oxygen  jet  (or  the  inner  edge  of  the  flame) 
and  the  outer  edge  of  the  flame.  An  image  taken  by  Mayer  and  Tamura  is 
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Figure  21:  Emission  above  400nm  at  60  bar,  near  injector 


Figure  22:  Oil*  radical  emission  at  45  bar  from  Mayer  and  Tamura,  Journal 
of  Propulsion  and  Po  wer  Vol  12,  1996 


Figure  24:  Shadowgraphy  at  60  bar,  near  injector 
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Figure  25:  Shadowgraphy  at  45  bar,  from  Mayer  and  Tamura,  Journal  of 
Propulsion  and  Power  Vol  12,  1996 


shown  in  Figure  25  for  comparison. 

An  Abel  inversion  of  the  near-injector  emission  images  at  60  bar  is  shown 
in  Figure  26.  This  will  not  be  valid  in  the  central  region  where  the  jet 
partially  shields  the  far  side  of  the  flame.  The  images  show  the  proximity  of 
the  flame  to  the  oxygen  jet.  The  bright  spot  behind  the  oxygen  injector  lip 
is  smaller  and  less  intense  than  that  observed  by  Mayer  and  Tamura  on  the 
1mm  diameter  injector. 

Important  points  can  be  made  about  the  flame  shape  at  60  bar.  The 
flame  and  the  oxygen  jet  are  extremely  turbulent.  The  jet  breaks  up  into 
packets,  rather  than  droplets.  The  flame  is  stabilized  on  the  oxygen  injector 
lip  in  a  recirculation  zone  and  then  grows  in  intensity.  This  is  schematized 
in  Figure  28. 


6  Conclusion 

Structure  of  subcritical  and  supercritical  cryogenic  flames  are  examined  using 
imaging  techniques.  The  flame  is  stabilized  in  the  near  vicinity  of  the  oxygen 
tube  and  initially  spreads  around  the  oxygen  core  as  a  corrugated  surface. 
The  momentum  ratio  has  a  decisive  effect  in  the  subcritical  cases.  It  has 
a  less  visible  effect  in  the  supercritical  case  but  this  aspect  needs  further 
study.  Subcritical  and  supercritical  flames  differ  notably  in  their  structure 
and  corresponding  rate  of  expansion. 


Figure  26:  Abel  inversion  of  the  near-injector  emission  images 


Figure  27:  Abel  transform  of  the  averaged  images  at  60  bar 


Figure  28:  Schematic  representation  of  the  flame  shape  at  60  bar 
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CNESOIa 


Coherent  anti-Stokes  Raman  Scattering  (CARS) 
Measurements  in  a  Shear-Coaxial  Cryogenic  Jet  Flame 
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Germany.: 


ObjeGtive 


•  improve  design  and  optimization  of  high  performance 
LOX-GH2  Rocket  Engines 


*  Development  of  optical  Diagnostics  for  combustion  studies 
at  high  pressure  (1 0  MPa) 


— ^  CARS  technique 

Instantaneous  measurements  of  temperature  -  concentration 
— ^  Heat  Release 


Turbulence  mixing  of  hot  reaction  products  and  Hydrogen 
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CNES13a 


M  ASCOTTE :  Comiustion;  ch 


Multiplex  CARS  technique  on  MASCOTTE 

-  Detection  of  Hj  (fuel  species) 

-  Single-shot  measurements  at  5  Hz  frequency 

O  150  measurements  during  a  run  {At=30s)  at  1  bar 
100  measurements  (20  s)  at  higher  pressure 

-  Temporal  resolution:  10  ns 

-  Planar  BOXCARS  arrangement  for  with  200  mm  focal  lens 

O  probe  volume:  1 .2  mm-long  x  1 00  pm  diameter 
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CNESOSa 


The  CARS  technique 

(coherent  anti-Stokes  Raman  spectroscopy) 


•  Three  laser  beams  are  applied  simultaneously 
to  the  gas  mixture 

-  two  pump  beams 

-  one  probe  beam 

-  A  new  laser  beam  is  created 

-  Its  amplitude  is  strong  when  the  difference 
(a)p-o)J  matches  a  Raman  frequency  a. 

•  Multiplex  CARS  spectroscopy 

-  Use  of  a  monochromatic  laser  for  top 
and  a  broadband  dye  laser  for  6)3 

-Instantaneous  Measurements 

-  Spatially  Resolved  Measurements 


1 

A  i  J 

[ 

“s  % 

frequency 

ON  ERA 


Multiplex  GARS  technique  on  MASCOTTE 


-  Detection  of  Hj  (fuel  species) 


-  Single-shot  measurements  at  5  Hz  frequency 

O  150  measurements  during  a  run  (At=30s)  at  1  bar 
100  measurements  (20  s)  at  higher  pressure 


-Temporal  resolution:  10  ns 

-  Planar  BOXCARS  arrangement  for  with  200  mm  focal  lens 
O  probe  volume:  1.2  mm-long  X  100  pm  diameter 
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Temperature  measurement  on  at  atmospheric  pressure 


Single-shot  spectrum  of  Hj 
Ydoppi.r»  rcoi.ision.rfor  all  rotational  lines 


Boltzmann  diagram 


ON  ER  A 


DLR 
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Probability  (%)  Probability  {%) 


P.D.F.  of  single-shot  temperature  at  1  bar 

Position:  Y=20  mm 


Temperature  (K)  Temperature  (K) 


mean  temperature  and  standard  deviation 


ONE  R  A 


P.D.F.  of  single-shot  temperature  at  1  bar 

X=1 80  mm  from  the  injector 
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Temperature  Profile  at  Y=20  mm 
P=1  bar 


H  ig  h  pressure  condition 


•  Dependance  on  the  rotational 
state  J 

•  influence  of  the  collisional 
partners  for  (Hj,  Oj,  H2O) 


•  CARS  line  depends  on  the 

linewidth  y  (Doppler  +  collisional)  ^  ^ 


»  Collisional  broadening 

4-  Dicke  effect  become  significant.  ^ 

N  0.4| 


E 


4173  4135  4097  4069  4021  3983 

Raman  shift  (cm‘^) 

Single-shot  spectrum  at  10  bar 


Impossibility  to  use  the  Boltzmann  Diagram  Technique 


^DLR 


ONERA 
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Collisional  broadening  (cm  '.amagat') 


Collisional  Linewidths  for  Hg-Hg  and  H2-H2O 


0.016 
0,014 
0,012 
0,010 
0,008 
0,006 
0.004 
0,002 
0,000 

0  500  1000  1500  2000  2500  5uuu  ^  ggg  ,jQ0Q  <1500  2OOO  2500  3000 

Temperature  (K)  Temperature  (K) 

Hi  i  Courtesy  of  X.  Michaut  and  H.  berger  M  M 

(University  of  Dijon) 


Temperature  measurement  at  high  pressure 

•  Temperatures  are  deduced  by  fitting  experimental  spectra  with 
the  theoretical  ones  calculated  for  different  temperatures  and 
for  different  molar  fractions  of  and  H2O 

•  Effects  of  O2  are  merged  with  those  of  H2 

(same  cpllisional  behavior  of  H2-H2  and  H2-O2 ) 


•  The  fitting  technique  used  is  a  simple  least  squares  where 
the  error  function  Is  given  by 
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CARS  Intensity 


High  Pressure  H2  CARS  spectra 


P=30bar  P=65bar 


ON  ERA 


Terriperatu re  evolution 


Condition  C 

Radial  position:  1 5  mm  from  the  flow  axis  onera 
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Principle  of  measurement 

•  Simultaneous  detection  of  two  species;  Hj  and  HjO 


CARS  spectrum  HjO  CARS  spectruni 


•  Data  processing 

-  H2O  concentration  estimated  from  the  intensity  of  the  HjO  CARS  spectra 

-  Species  composition  for  the  cpllisional  linewidth  can  be  estimated 

-  Temperature  deduced  from  the  Ha  and  Hp  spectra!  profiles 


•  Simultaneous  detection  of  Hj  and  HaO  using  CARS  have  been  performed 
for  pressures  up  to  65  bar. 
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Conclusion 


•  Detection  of  H2  (and  HjO)  using  CARS  have  been  successfully 
performed  for  pressures  up  to  65  bar  in  a  cryogenic  corribustor. 


•  Temperature  and  concentration  can  be  deduced  from 
the  experimental  CARS  signals 


•  Data  base  for  the  validation  of  combution  solvers 


•  Measurements  will  be  peiTormed  on  the  P8  combustor  (DLR) 
at  60  bar  next  year 
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Statistical  Investigation 
of  the  Turbulent  Flame  Geometrical  Structures 
in  a  Liquid  Oxygen  /  Gaseous  Hydrogen 
Shear-Coaxial  Jet 


Armelle  CESSOU,  Pascal  COLIN,  Denis  STEPOWSKI 

CORIA  UMR  6614 

CNRS  -  Universite  &  INSA  de  Rouen 
Place  Emile  Blondel 
76821  Mont  Saint  Aignan,  France 


Context 


•  rocket  engine  application 

•  cryogenic  propellant  combustion  facility  : 

MASCOTTE  test  bench  olOUERA 

•  developed  for  better  understanding  of  physical 
phenomona  encountered  in  rocket  engines : 

•  two-phase  combustion 

•  Investigation  by  planar  laser  indMced  f luorespence 

c64^ja 
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Investigated  Area 

investigated  area _ 


Flame  Description 


•  Flame  always  anchored  at  the  lip  of  the  injector 

•  Non-premixed  character 

•  conserved  scalar :  mixture  fraction 


Description :  constant-property  surface 
Hyp. :  infinitely  fast  chemistry  A  flame  at  stoichiometric 


constant-property  surface  Z  =  Zj 


Flame 


358 


Flame  Description 

•  chemistry  decoupled  from  turbulence 


mean  reaction  rate 
_ p.u.v. 


mean  reaction  rate 
p.u.a. 


flame  surface 


p.u.v. 


•  local  geometry  of  constant-property  surface 

•  presence  probability 

•  normal  to  the  surface 

•  curvature 

_ ■  flame  surface  production 


Operating  Conditions 


Point 

g/s 

g/s 

P 

bar 

TgH2 

K 

V,ox 

m/s 

VgH2 

m/s 

Roi 

xIO® 

J 

We 

x10® 

image 

number 

Cl 

1 

289 

2.18 

620 

67 

151 

C1p 

50 

9.5 

8 

289 

2.18 

163 

67 

M 

f482|,; 

:''fi 

1 

94 

2.18 

419 

16.9 

274 
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Experimental  Set-up 


Fluorescence  Scheme 


Effective  absorption  cross  section 
Fluorescence  efficiency 


>trongly  Increase  with  temperature 
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O _ _ 


Instantaneous  Image 


Flame  Identification 


•  reverse  configuration 

•  at  adiabatic  flame  temperature 

laser  irradiance  not  altered 

detection  limit :  Z  =  0.109 

before  reaching  the  flame 

while  Zs  =  0.11 

Hyp. :  infinitely  fast  chemistry  for  H2/O2 


bottom  outiine  of  fluorescence 
flame 


Contourline  Determination 


Binarization  S,  >  threshold  ■=>  S  =  1 00 
Sf  <  threshold  ‘=>8  =  0 
median  filter  (3  x  3) 
contourline  detection 


jjjllllllllllllHIlH 


High  frequency  filtering 


-  sampled  at  evenly  spaced 
intervals 

ds  =  +  d/ 

-  smoothed  through  a 
gaussian  filtering  (800  pm) 
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High  and  Low 
Frequency  Filterings 


high  freq.  low  freq. 

800  ii  m  profile  length 


UMR  6614 

j  CNRS  *UNIVERSITE  et  INSA  de  ROUEN  I 


^  -1,4 


2  /-■  ;:3 

x/d 


Presence 

Probability 

•  mesh  0.25  X  0.25  mm^ 

•  X  <  d,  sharp  distribution 

•  X  <  2dj 

fiamebows 
large  distribution 


0.0  0,2  0,4  0,6 


Length  Increase 


D) 

S  0.4 


a :  flapping 
b  :  flapping  +  wrinkling 

-  C1 


Curvature 


UMR  6614 
5  -UNIVERSITE  et  INSA  de  ROUEN  I 


6  -4  -2  0  2  4  6  -4  -2  0  2  4  6  -4  -2  0  2  4  6  8 


x(Si)y(Si)  -  y(S|)x(S|) 
(x2(Si)  -H  y^(Si)y^ 


UMR  6614 

I  CNRS  -UNIVERSITE  et  INSA  de  ROUEN  I 


Density 


DX 


•  2:  =  Px<i/dilZ3> 

•  mesh  0.25  x  0.25  mm^ 

• .  x<  2di 

S  underestimated 

froml 


UMR6614 

CNRS  -UNIVERSITE  et  INSA  de  ROUEN 


Comparison 


0  1  2  3  4  5 


x/d, 


0.02  0.32  0.62  0.92  1.22  1.52 


Conclusion 


-  reverse  configuration  N 

(oxygen  within  the  flame  envelope)  ^  stoichiometric  iSO-SUrface 


-  O2  fiuorescence 


statistical  description  of  flame  surface  properties 


validation  tools  for  numerical  modeling 


UMR  6614 
CNRS  -UNIVERSITE  et  INSA  de  ROUEN 
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I  Introduction 

II  Bifurcation  of  Dynamical  States  in  a  Dump 
Combustor 

III  Bifurcation  of  Steady  Combustion  with 
Recirculation  Zones 

IV  Incorporating  Bifurcation  of  Combustion  in 
Analysis  of  Combustor  Dynamics 

V  Some  Consequences  of  Subcritical  Bifurcations 
with  Noise 


VI  Concluding  Remarks 


HISTORICAL  BACKGROUND 


Bifurcation  of  Dynamical  States  of  Combustion 

In  the  mid- 1 980 ’s,  hysteresis  was  discovered  at  the 
boundary  between  stable  and  unstable  operation  of 
the  Caltech  dump  combustor  (Smith,  1985;  Sterling, 
1987) 


>  These  results  suggest  bifurcation  of  dynamical 
states  of  combustion  (limit  cycles?) 

>  The  behavior  has  since  been  found  by  others 
[Cohen  (UTRC);  Richards  (METC)] 


HISTORICAL  BACKGROUND  (conf  d) 


Bifiircation  of  Steady  States  of  Combustion 

•  In  the  1970’s,  Natanzon  et  al  proposed  and  investigated 
the  stability  of  bifurcated  steady  states  of  combustion  as 
a  possible  mechanism  for  random  occurrences  of 
combustion  instabilities. 


>  The  reasoning  rests  on  bifurcation  of  steady 
combustion  in  a  recirculation  zone  (Natanzon  et  al, 
1992) 


Inlet  Temperature  To 


HISTORICAL  BACKGROUND  (cont’d) 


>  One  state  (cold  recirculation  zone)  is  more  unstable  and 
prone  to  produce  combustion  instabilities. 

>  Which  state  is  reached  depends  on  history  (e.g.  during 
ignition  or  shut-down);  the  final  state  depends  on 
random  ‘accidents’  of  that  history. 

>  Hence  random  occurrences  of  combustion  instabilities 
may  be  observed. 


>  Analysis  of  the  recirculation  zone  is  an  extension  of 
that  for  a  perfectly  stirred  (homogeneous)  reactor 
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Radiation  Intensity 


HISTORICAL  BACKGROUND  (conf  d) 

7 

•  Doulatov  et  al  (1998)  (independently)  proposed 
existence  of  ‘multiple  modes’  of  combustion,  based  on 
statistical  analysis  of  pressure  records 

.  .statistical  error  equaled  almost  100%  while  the 
instrumentation  error  was  no  more  than  15%. . .” 

>  Observations  of  flame  structure  in  a  small  NH3/O2 
engine  (Beliy  et  al,  1983;  Abulkanov  et  al,  1974) 
showed  correlation  between  presence  of  instabilities 
and  spatial  characteristics  of  the  flame  stmcture. 


X ,  mm 

Mono-top 
Stable  Flame 


X ,  mm 

Migratory 

Flame 


X  ,  mm 

(1)  Mono-top  Stable  Flame 

(2)  Bi-top  flame  with 
Spontaneous  Instability 


Dependence  of  radiation  intensity  (ICuH)  on  the  chamber 
length;  a)  typical  mono-top  stable  flame;  b)  an  example 
of  migratory  flame;  c)  mono-top  stable  flame  and  bi-top 
flame  created  by  the  spontaneous  instability  in  the 
combustion  chamber  from  a  single  experiment. 
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HISTORICAL  BACKGROUND  (cont’d) 


•  Note  the  connection  in  the  Russian  reports  between 
statistical  analysis  of  pressure  records  (both  multiple  and 
single  firings)  and  the  conclusion  that  bifurcation  of  the 
combustion  process  exists  [Ref  1st  International 
Symposium  on  Liquid  Propellant  Rocket  Engines, 
Pennsylvania  State  University  (1993)] 

•  In  numerical  simulations  it  may  be  difficult  to  identify 
bifurcations  except  with  special  intentions  (?) 

•  Presence  of  hysteresis  constitutes  a  basis  for  control  by 
injecting  pulses  to  cause  transitions  between  two  stable 
states. 
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SOME  REMARKS  ON  AFFECTING 
COMBUSTION  INSTABILITIES 

The  problem  of  linear  stability  of  motions  in  a  combustion 
chamber  can  be  reduced  to  a  general  form  of  Rayleigh’s 
Criterion. 

Rayleigh’s  original  criterion  was  based  on  excitation  of 
acoustic  waves  in  a  tube  (e.g.  a  Rijke  tube)  by  energy 
addition  in  the  form  of  heat  transfer  or  energy  released  in 
combustion  porcessses. 

The  criterion  states  (roughly)  that  energy  addition  tends  to 
destabilize  (i.e.  encourage)  a  pressure  disturbance  if  the 
rate  of  energy  addition  is  increasing  locally  where  the 
pressure  distrubance  is  also  increasing. 

Q :  rate  of  energy  addition/volume 

— — ;  component  of  rate  of  change  of  Q  in  phase 
dt 

with  the  pressure  disturbance  hp 

f 

Change  of  energy  if  the  disturbance  ~  — — bp 


Formally: 


df  a'  dt 

n-Wp'=-f 


SOME  REMARKS  ON  AFFECTING 
COMBUSTION  DJSTABILITIES  (cont’d) 

Modal  expansion  and  spatial  averaging: 

p'(nt)=plp,(t)v,(0 : 

at 

F.=-^{\h'V.dV^^fiV.dS]  El  =  \^ldV 

Change  of  energy  of  an  oscillator  in  one  cycle 

^  iT'  r^n' 


‘Rayleigh’s  Criterion” 

AS.=^^Jn/.dF  |fi.^-dt’>0  =>  instability 

P^n  ‘ 
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SOME  REMARKS  ON  AFFECTING 
COMBUSTION  INSTABILITIES  (cont’d) 

•  The  destabilizing  tendency  of  energy  addition  depends  on  the 
local  phase  relation  and  spatial  disturbance. 

•  Stabilizing  or  controlling  combustion  instabilities  “according 
to  Rayleigh’s  criterion”  means  altering  the  spatial  distribution 
and  phasing  of  the  energy  released  by  combustion. 

•  e.g.  modulate  the  primary  or  a  secondary  supply  of  fuel. 

•  Neither  successes  nor  failures  in  treating  combustion 
instabilities  in  liquid  rockets  are  well-understood;  notably  the 
amplitudes  of  oscillations  before  and  after  application  of  a 
method  to  improve  stability  have  not  been  predicted. 

•  Why?  Because  the  mechanisms  are  not  well- 
understood:  relying  only  on  Rayleigh’s  criterion  for 
interpreting  problems  of  combustion  instabilities 
produces  no  information  about  the  mechanisms  unless 
a  theory  for  dQ/dt  is  given. 

•  The  situation  is  somewhat  better  for  solid  rockets 
because  the  overwhelming  dominant  mechanism  is 
associated  with  the  dynamics  of  combustion  processes 
at  a  burning  solid  surface.  A  secondary  cause  of 
instabilities  is  vortex  shedding,  strongly  dependent  on 
details  of  the  internal  geometry. 

•  Consideration  of  bifurcations  means  paying  attention 
to  a  special  class  of  mechanisms. 
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II  BIFURCATION  OF  DYNAMICAL 
STATES  IN  A  DUMP  COMBUSTOR 
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Combustion  Science  and  Technology  V.  126, 
pp.  381-388. 
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POTENTIAL  SIGNIFICANCE  OF 
HYSTERESIS 


In  combustors  for  which  hysteresis  exists,  the  range  of  stable 


If  a  disturbance  causes  the  system  to  be  unstable  in  state  A, 
active  control  of  pulsing  can  produce  transition  to  the  stable 
state  B. 


Understanding  hysteresis  means  understanding  true  nonlinear 
behavior  and  may  lead  to  a  method  of  nonlinear  control  in  the 
regions  of  vigorous  instabilities,  possibly  more  efficient  than 
methods  of  line  or  control. 

However,  when  it  exists,  hysteresis  occurs  in  some  lean  range  of 
mixture  ratio,  i.e.  in  the  region  including  the  lean  blowout  limit. 

Hence  a  practical  strategy  for  design  of  operational  combustors 
maybe 

•  design  the  combustor  for  stable  operation  for  moderate  and 
high  mixture  ratios,  by  using  'passive*"  control  (geomeby,  flow 
field  ...) 

•  use  active  control  based  on  the  presence  of  hysteresis  in  the 
range  of  low  mixture  ratios 


QUESTIONS  ON  THE  ROLE  OF 
HYSTERESIS  IN  ACTIVE  CONTROL  OF 

COMBUSTORS 


•  How  general  is  hysteresis  in  combustors? 

•  Under  what  conditions  does  hysteresis 
exist? 

•  What  distinguishes  the  conditions  under 
which: 

(i)  a  pulse  causes  no  transition?;  and 

(ii)  a  pulse  generates  transition  to  the 
lower  branch? 


MEASURED  STABILITY  BOUNDARY  AND 

HYSTERESIS  LOOP 
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EXAMPLE  OF  A  PULSED  TRANSITION 
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Ill  BIFURCATION  OF  STEADY  STATE 
COMBUSTION  WITH 
RECIRCULATION  ZONES 


References: 
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Natanzon,  M.S.  and  Men’shikova,  O.M.  (1992) 

“Bifurcation  of  Steady  Combustion  Regimes 
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Chambers,”  Physics  of  Combustion  and 
Explosion,  Vol.  28,  No.  4,  (pp.  10-18). 

Natanzon,  M.S.  (1999)  Unsteady  Combustion  in  Liquid 
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3.1  RTFTJRCATTON  OF  COMBUSTION  IN  A 
PERFECTLY-STIRRED  REACTOR 

References: 
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ZePdovich,  Y.  B.  Barenblatt,  G.  1.  And  Librovich,  V. 

B.  (1985)  Mathematical  Theory  of  Combustion 
and  Explosion. 

Bifurcations  in  Perfectly-Stirred  Reactors: 

Vulis,  L.  A.  (1960)  Thermal  Regimes  of  Combustion, 
McGraw-Hill  Book  Company,  New  York. 

Longwell,  J.  P.  and  Weiss,  M.  A.  (1955)  “High 

Temperature  Reaction  Rates  in  Hydrocarbon 
Combustion,”  Industrial  and  Engineering 
Chemistry,  Vol.  47,  pp.  1634ff-1643. 
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3.1  BIFURCATION  OF  COMBUSTION  IN  A 
PERFECTLY-STIRRED  REACTOR  (cont’d) 


References: 

Bifurcations  in  Perfectly-Stirred  Reactors  (cont’d): 

Hottel,  H.  C.,  Williams,  G.  C.  and  Bonnell,  A.  H. 

(1958)  “Application  of  Well-Stirred  Reactor 
Theory  to  the  Prediction  of  Combustor 
Performance.” 


•  Spalding  and  Hottel  et  al  treat  a  model  of  a  combustor 
comprising  a  perfectly  stirred  reactor  (mixing  of  reactants 
and  recirculated  products;  no  combustion)  followed  by  a 
plug  flow  reactor  (recirculation  zone  with  combustion). 


3.1  BIFURCATION  OF  COMBUSTION  IN  A 
PFBFF.CTT.Y-STIRRED  reactor  (cont’d) 


Reference:  Vulis,  L.  A.  Thermal  Regimes  of  Combustion 


„  m 
T  ,  c°  — 


T,c 


T,c 


Chemical  Reactor 


Assumption: 


Instantaneous  complete  mixing  of  inlet 
flow,  producing  uniform  conditions  in  the 
volume. 


mass  of  reactants  in  volume  _  p  /  >> 

T  :  residence  time,  r  =- - - - 

^  mass  flow  rate  m 

C.— C/  .,  \ 

w:  reaction  rate,  w= - ) 

Q ;  heat  of  combustion  per  unit  mass  (Energy/mass) 

Assumption:  Ignore  heat  losses  through  the  boundary 

condition  for  steady  reaction  (combustion). 

c:  mass  concentration  of  reactant  mixture 
V ;  volume  of  reactor 
p:  density  of  products 


3.1  BIFURCATION  OF  COMBUSTION  IN  A 
PERFECTLY-STIRRED  REACTOR  (cont’d) 


Rate  of  Energy  Release 
by  Combustion 


Net  Rate  of  Energy 
Loss  with  Flow 


0^  =  (per  unit  mass  of  reactants 

in  the  volume) 


Qi  -  ~ — ~Q  -  [^0^^  ]  C  for  a  first  order 

reaction 


^c(T-T^)  1 


//  rr 

pV 


(T-Tj 


Rewrite  the  condition  in  terms  of  four  dimensionless  parameters 


1 

CD 

(temperature) 

1 

O 

CD 

(inlet  temperature) 

^_Rc^Q 

^  c  E 

(heat  of  combustion) 

T  =T  A:  =— 

rc  r  0 

T. 


^residence  time  ^ 
V  reaction  time  J 


3.1  BIFURCATION  OF  COMBUSTION  IN  A 


PF.RFFCTLY-STIRRED  REACTOR  (cont’d) 


Define:  <p=l —  concentration  of  reacted  mixture  leaving  the  reactor, 
c. 


Complete 

combustion 


No  combustion 


Replot  to  show  hysteresis  loop: 


Temperature  in 
Reactor,  T* 


Extinction  ! 


^Stable  Branch 


Unstable  Branch 


Ignition 


Inlet  Temperature  T^ 


Subcritical  bifurcation  with  turning  points  occurs  because 

1)  Energy  loss  is  a  linear  function  of  temperature 

2)  Energy  gain  is  a  nonlinear  function  of  temperature 
having  a  ‘special’  form 
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3.2  BIFURCATION  OF  STEADY  COMBUSTION 
WITH  RECIRCULATION  ZONES; 
Connections  With  Loss  of  Stability 

Reference:  Natanzon  et  al. 

Motivation 

•  In  general  there  are  two  classes  of  mechanisms  for  loss  of 
stability: 

•  linear  instability:  ‘soft  excitation’,  supercritical 

bifurcation 

•  nonlinear  instability:  ‘hard  excitation’,  subcritical 

bifurcation 

•  Instabilities  arising  from  subcritical  bifurcations  may  be 
generated  in  two  ways: 

1)  sufficiently  large  disturbances  applied  to  the  system 
operating  in  the  lower  branch; 

2)  during  transient  changes  the  system  may  be  taken  to 
the  lower  turning  point,  then  making  the  transition  to 
the  upper  branch: 
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3.2  BIFURCATION  OF  STEADY  COMBUSTION 
WITH  RECIRCULATION  ZONES; 
Connections  With  Loss  of  Stability  (cont’d) 

•  If  instabilities  are  found  far  from  the  linear  stability 
boundary,  a  possible  cause  is  subcritical  bifurcation  of 
combustion. 

•  In  the  absence  of  large  disturbances,  path  (D  may  be 
the  cause,  identified  by  connection  with  the  history  of 
the  operating  state. 

Idea: 

-  •  Bifurcation  of  steady  combustion  is  associated  with 

behavior  of  recirculation  zones  treated  as  one¬ 
dimensional  perfectly  stirred  (homogeneous)  reactors. 

•  The  length  of  the  recirculation  zone  becomes  the 
bifurcation  parameter. 

•  An  essential  assumption  is  that  local  ignition  of  the 
injected  reactants  (gas/gas  or  liquid/gas)  may  occur  in 
two  ways. 

1)  presence  of  hot  combustion  products  from 
the  recirculation  zone; 

2)  ‘self-ignition’  of  the  injected  reactants  due  to 
high  initial  temperature  at  the  injector  outlet, 
or  low  activation  energy. 


3.2  BIFURCATION  OF  STEADY  COMBUSTION 
WITH  RECIRCULATION  ZONES: 

Connections  With  Loss  of  Stability  (cont’d) 

•  The  two  cases  of  ignition  produce  different  spatial 
structures  of  the  combustion  field  in  the  vicinity  of  the 
injector. 

•  If  one  of  the  fields  is  linearly  unstable,  then  the 
bifurcation  of  steady  combustion  leads  to  an  apparent 
bifurcation  of  stability. 


Ignition  by  Recirculation  Zone 


DISTANCE  FROM  INJECTOR  FACE  (mm) 


Self-Ignition  in  the  Bulk  Flow 
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3.3  NATANZON’S  ANALYSIS  OF  COMBUSTION 
WITH  RECIRCULATION  ZONES 

•  The  analysis  is  based  on  a  simplified  physical  model: 
the  results  should  be  regarded  as  suggestive  of  the 
behavior  described  above 

Physical  Model: 


-  C  - 

RECIRCULATION 

ZONE 


•  The  field  is  approximated  as  a  one-dimensional  steady 
flow: 


•  Recirculation  zone  is  followed  by  the  combustion  region 
downstream. 
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3.3  NATANZON^S  ANALYSIS  OF  COMBUSTION 
WITH  RECIRCULATION  ZONES  (cont’d) 

•  A  condition  for  steady  combustion  in  the  recirculation 
zone  is  set  in  analogy  with  the  condition  Qj  =  Qjj 
constructed  for  a  perfectly  stirred  reactor,  but  flow  enters 
due  to  entrainement  and  leaves  at  the  upstream  face, 
having  the  temperature  7^  of  the  zone. 

•  For  analysis  of  the  recirculation  zone,  the  external  flow 
from  the  injection  port  is  implicit,  being  the  source  of  the 
entrained  flow. 

•  As  for  the  perfectly  stirred  reactor,  there  exist  conditions 
under  which  steady  states  exist:  one  is  stable,  and  two  are 
unstable. 

i)  Low  temperature  stable  state  with  ‘no’ 
combustion  in  the  recirculation  zone.  The 
main  flow  must  ‘self-ignite’  at  the  end  or 
downstream  of  the  recirculation  zone. 

ii)  High  temperature  stable  state  representing 
combustion  in  the  recirculation  zone.  The 
flow  leaving  the  zone  at  the  injector  face 
ignites  the  main  flow. 
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3.3  NATANZON’S  ANALYSIS  OF  COMBUSTION 
WITH  RFCTRCULATION  ZONES  (cont’d) 


Condition  for  steady  combustion: 


(pi,Cj)'  -  ipuCl) 

Net  Energy  Flow 
out  of 

Recirculation 

Zone 


=  Q„iT’) 

=  wH/  - 

Energy 

Released 

by 

Combustion 


Energy 

Absorbed 

by 

Liquid  Fuel 
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3.3  NATANZON’S  ANALYSIS  OF  COMBUSTION 
WITH  RECIRCULATION  ZONES  (cont’d) 


Equations  Solved  for  Recirculation  Zone 


Gas  Flow: 


d_ 

dx 


:  vaporization  rate 

:  mass  concentrations 


A{pu)=GJ 


v,+v„+v  =1 

f  0  p 


Combustion  Products: 


dx 


(p^,)= 


w 


V  A(pii)=wi 


A(pw)=(pw)*  -(pi/)° 


Fuel  Vapor: 


—{piiv.)=G^  — —w 
dx^^  \+t 


:  vaporization  rate 
of  fuel  drops 
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3.3  NATANZON’S  ANALYSIS  OF  COMBUSTION 
WITH  RECIRCULATION  ZONES  (cont’d) 

Energy:  ^ 

Upuh) = wH,  -  G,.  [c,  (r  -  r;  )+§(r )]+ g,c/, 

ax 

r  :  surface  temperature  of  fiiel  drops 
:  heat  capacity  of  fuel 
heat  of  vaporization 


(puC  t) -{puC  l)  =  wH/-B^{pu) 


Motions  of  Fuel  Drops: 


b~c,{t-t:)v4t,)-ct, 

l-.=l+hjc. 


dw 

dx 


v:  velocity 
m:  mass 
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3.3  NATANZON’S  ANALYSIS  OF  COMBUSTION 
WITH  RECIRCULATION  ZONES  (cont’d) 

•  The  true  meaning  of  the  model  and  the  calculations  is 
NOT  clear  in  the  two  references  but  it  appears  that  the 
model  is  virtually  identical  to  that  for  a  PSR. 

•  Solution  is  obtained  by  iteration  giving  the  results 
sketched  above. 

•  The  results  of  the  calculations  for  the  recirculation  zone 
can  be  displayed  in  the  bifurcation  diagram,  a  replot  of  the 
results  given  above: 

2000 

* 

t- 

.  1500 

H 

N 

g; 

I  1000 


I  500 

V 


0  5  10  15  20  25  30 

Length  of  Recirculation  Zone,  I  (mm) 

•  This  figure  shows  that  for  initial  temperatures  T^  higher 
than  about  830°K  only  the  high  temperature  state  of 
combustion  exists  in  the  recirculation  zone. 
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3.3  NATANZON’S  ANALYSIS  OF  COMBUSTION 
WITH  RECIRCULATION  ZONES  (cont’d) 

•  For  a  range  of  T  °  <  830°K,  say  T  ^  =  775°K,  two  stable  states 
are  possible  for,  in  this  case,  1.2cm  <  ^  <  2.0cm. 

-  For  r  °  =  775°K,  as  increases,  the  temperature  remains 
low  in  the  recirculation  zone  until  t  «  2cm  when  ignition 
occurs. 

-  Then  t  is  reduced,  combustion  is  sustained  until  t  »  1.2cm 
when  extinction  occurs. 

•  Multiple  states  can  also  be  obtained  if  for  a  fixed  length  of 
recirculation  zone  the  equivalence  ratio  a  (i.e.  the  F/0  mixture 
ratio)  is  changed,  consequently  changed  the  heat  of 
combustion. 

The  next  figure  shows  the  bifurcation  diagram  in  the 
coordinates:  temperature  7^  in  the  recirculation  zone 
versus  (|),  the  equivalence  ratio. 
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Equivalence  Ratio,  ^ 

•  Note  that  in  this  setting  two  hysteresis  loops  exist. 
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3.3  NATANZON’S  ANALYSIS  OF  COMBUSTION 
WITH  RECIRCULATION  ZONES  (cont’d) 

Combustion  and  The  Temperature  Field  Outside 
the  Recirculation  Zone 

•  Natanzon  and  Men’shikova  (1992)  have  reported  the 
following  results  for  the  temperature  field  outside  the 
recirculation  zone,  obtained  by  solving  one-dimensional 
equations  with  T^  =  775°K  and  (|)  =  ? 
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3.3  NATANZON’S  ANALYSIS  OF  COMBUSTION 


WITH  RECIRCULATION  ZONES  (cont’d) 


•  The  sequence  of  these  cases  falls  as  shown  in  the  bifurcation 
diagram: 


500  h 


Length  of  Recirculation  Zone,  t  (mm) 

Some  of  the  detailed  behavior  of  these  curves  depends  on  the 
model  of  droplet  combustion,  fragmentation,... 

Note  that  the  overall  steady  combustion  process  in  the 
chamber  is  hardly  affected  by  the  particular  condition  in  the 
recirculation  zone: 

~  if  the  chamber  is  sufficiently  long,  the  combustion 
efficiency  and  exit  temperature  are  essentially  identical 
in  all  cases 

In  respect  to  the  excitation  of  combustion  instabilities, 
when  there  are  multiple  steady  states  ,  the  question  is: 

Which  of  the  steady  states  is  less  stable? 
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3.3  NATANZQN’S  ANALYSIS  OF  COMBUSTION 
WITH  RECIRCULATION  ZONES  (cont’d) 

•  Because  the  current  existence  of  a  particular  steady  state 
depends  on  the  history  of  the  operating  point  of  the 
chamber,  the  possibility  for  ‘randomly’  occurring 
instabilities  is  exposed. 

A  Two-Dimensional  Model  for  Combustion  with  Bifurcations 

•  A  model  is  constructed  for  gas/gas  reactants. 

•  The  recirculation  zone  is  treated  as  a  perfectly  stirred  reactor 
but  flows  in  and  out  are  included. 

•  A  condition  for  stationary  combustion  enforced. 

•  Details  of  the  flow  outside  the  recirculation  zones  are 
incomplete. 

•  Results  show  the  existence  of  two  states  of  steady  combustion: 

16 
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Ignition  by  Recirculation  Zone 
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3.3  NATANZON’S  ANALYSIS  OF  COMBUSTION 
WITH  RF.CTRCULATION  ZONES  (cont’d) 


Self-Ignition  in  the  Bulk  Flow 


Stability  of  the  Two  States  of  Steady  Combustion 

•  Natanzon  establishes  the  stability  of  the  two  steady  states  by 
applying  his  method  based  on  the  amplitude-phase-frequency 
characteristic  (APFC)  interpreted  as  a  “hodograph.” 

•  A  recirculation  zone  is  a  passiye  capacitance,  proyiding 
distributed  sources  along  its  boundary. 
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3.3  NATANZON’S  ANALYSIS  OF  COMBUSTION 


WITH  RECIRCULATION  ZONES  (cont’d) 

•  The  method  comes  down  to  examining  the  behavior  of  the 
transfer  function  R{m)  for  the  combustion  zone: 

R(i(i))=— 

8p 

bu:  velocity  fluctuation  at  the 
bp:  pressure  fluctuation 


6  —  r -  Cold 

\  Recirculation 
\  Zone 

\  0o)  =  2:fl5OO 

V 


=2n.500  ^ 


\(o-2;r-l000 


271-1000  2;r '450^^271 -375 


o)  =  2n-l500 


-101 _ ^ ^ _ L 

-10  -8  -6  -4 


4  6  8 


•  IRe  —  >0  =>  tendency  to  instability.  (Natanzon) 

\^P) 

•  Consistent  with  Rayleigh:  move  Q  where  p  is  high 


HENCE  the  state  with  a  hot  (combusting) 
recirculation  zone  is  more  unstable. 


3.3  NATANZON\S  ANALYSIS  OF  COMBUSTION 
WTTH  RECIRCULATION  ZONES  (cont’d) 


Experimental  Results  Confirming  the  Possible 
Influences  of  Recirculation  Zones 

•  The  figure  shows  laboratory  apparatus  specially  constructed 
to  demonstrate  qualitative  changes  of  stability  when 
recirculation  zones  are  present. 

•  Tests  were  run  with  and  without  recirculation  zones; 
recirculation  zones  were  created  by  mounting  a  cascade 
of  rods  immediately  downstream  of  the  pylons  through 
which  fuel  was  injected  in  an  air  stream. 


ZONES  ZONES 


In  all  tests,  the  chief  experimental  variables  were  the 
temperature  7^  of  the  inlet  air  and  the  equivalent  ratio. 


The  reactants  were  liquid  fuel  and  heated  air. 
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3.3  NATANZON’S  ANALYSIS  OF  COMBUSTION 


WITH  RECIRCULATION  ZONES  (cont’d) 


combustion  and  no  ignition  or  extinction. 


- boundary  of  extinction 

A  steady  combustion 
A  no  combustion 


•  The  extinction  boundary  is  essentially  coincident  with  the  boundary 
of  self-ignition  (apparently  no  ignition  was  used)  and  is  minimum 
when  (])  =  1. 

•  Tests  with  turbulence-inducing  cascade  (and  recirculation  zones) 

-  boundaries  of  extinction 

O  steady  combustion 
#  no  combustion 

B  region  in  which  hysteresis  is  observed 


IV  INCORPORATING  BIFURCATION  OF 
COMBUSTION  IN  ANALYSIS  OF 
COMBUSTOR  DYNAMICS 
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AGARDCP512. 

Burnley,  V.  (1990)  “Nonlinear  Combustion  Instabilities 
and  Stochastic  Sources,”  Ph.D.  Dissertation, 
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INCORPORATING  A  GLOBAL  BIFURCATION  OF 
COMBUSTION  IN  ANALYSIS  OF  COMBUSTOR 

DYNAMICS 


•  Two  apparently  different  types  of  subcritical  bifurcations 
have  been  described,  both  involving  recirculation  zones  in 
essential  ways. 

1)  Bifurcation  of  dynamical  states,  showing  the  existence 
of  two  steady  oscillatory  states  in  a  dump  combustor: 


A:  loss  of  linear  stability  and  development 
of  steady  limit  cycles  for  all  ^ 

B:  absolute  linear  stability  for  ([)  >  (l)^ 

Lower  branch:  weak  oscillations:  stable  linear  modes 

driven  by  noise  sources  independent  of 
the  acoustic  field. 

Upper  branch:  strong  oscillations  in  stable  limit 

cycles  characterized  by  the  balance 
of  energy  (power)  gains  and  losses; 
the  behavior  is  intrinsically 
nonlinear  (gains  of  losses  or  both). 
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TNCORPORATING  A  GLOBAL  BIFURCATION  OF 
COMBUSTION  IN  ANALYSIS  OF  COMBUSTOR 

DYNAMICS  (cont’d) 


2)  Bifurcation  of  steady  states,  one  of  which  is  linearly 
stable  and  one  linearly  unstable,  developing  stable 
limit  cycles. 


•  The  two  types  of  bifurcation  involve  physical  behavior 
differing  in  detail  but  which  can  be  formally  accommodated 
in  the  same  theoretical  framework. 


•  Note  that  the  scales  of  the  recirculation  zones  differ  in 
the  two  cases  because  the  scales  of  the  combustor  differ. 


•  In  liquid  rockets,  many  ‘small’  recirculation  zones  exist 
at  the  injector  face.  The  problem  of  cooperative 
behavior  to  produce  a  bifurcation  of  dynamical  behavior 
in  the  combustor  remains  unsolved:  interactions  among 
the  zones  is  sufficient  and  probably  necessary. 
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INCORPORATING  A  GLOBAL  BIFURCATION  OF 
COMBUSTION  IN  ANALYSIS  OF  COMBUSTOR 

DYNAMICS  (cont’d) 


•  Analysis  based  on  modal  expansion  and  spatial  averaging 
suggests  the  strategy  to  follow. 

Modal  Expansion  p'(f^t)=?Zihi(0¥i(^) 

Spatial  Averaging  =F, 

at 

F.=~{\h'VdV^^h,dy} 

P^n 

where  h,f  represent  volumetric  and  surface 
sources  of  pressure  distrubances: 

;  «-Vp'=-/ 

a  dv 


•  The  functions  h  and/ contain  purely  gasdynamical  parts  and 
sources  of  mass,  momentum,  and  energy  due,  primarily,  to 
combustion  processes  and  interaction  between  gaseous  and 
condensed  phases. 

•  Some  of  the  sources  may  be  due  to  actuation  associated 
with  some  means  of  active  control. 
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INCORPORATING  A  GLOBAL  BIFURCATION  OF 


COMBUSTION  IN  ANALYSIS  OF  COMBUSTOR 


DYNAMICS  (cont’d) 

•  Gasdynamics  alone  cannot  produce  subcritical  bifurcations. 
Hence  the  behavior  sought  must  be  associated  with  other 
sources,  in  particular  the  dynamics  of  the  energy  release 
accompanying  combustion  processes. 

•  As  a  result  of  the  spatial  averaging,  the  dynamical 
consequences  of  volumetric  and  surface  sources  (h,  f) 
are  equivalent. 


a  dt 


n-Vp'  =  -f 


Modal  Expansion  and  Spatial  Averaging: 

^+co;n.=0;  F,=^-j^{\h^v.dV+<^fy,.dS} 

With  volumetric  energy  release  2^  and  mass  injection  at 
the  lateral  boundary: 

Accounting  for  volumetric  energy  addition  only: 
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TNCORPORATING  A  GLOBAL  BIFURCATION  OF 
COMBUSTION  IN  ANALYSIS  OF  COMBUSTOR 

DYNAMICS  (cont’d) 


•  Two  models  of  Q[  proposed: 

1)  Periodic  Delayed  Combustion  Response 

Physical  Basis:  (observations) 

(i)  For  small  amplitudes,  a  small  shed  vortex  is  surroimded 
by  hot  combustion  products  from  the  previous  cycle: 
shorter  time  delay  and  less  energy  released. 

(ii)  For  large  amplitudes,  a  large  shed  vortex  interacts  with 
the  wall,  creating  a  region  of  cold  reactants:  a  longer 
time  delay  is  associated  with  larger  energy  release. 

Model:  Energy  released  in  the  vortex  shedding  is 
proportional  to  the  acoustic  velocity,  with  a 
time  delay  (relative  to  the  pressure  oscillation) 
linear  in  the  acoustic  oscillator  energy 


TiyrORPORATING  A  GLOBAL  BIFURCATION  OF 


COMBUSTION  IN  ANALYSIS  OF  COMBUSTOR 


DYNAMICS  (cont’d) 


Realization: 


Q' 

ri  =  V^  sinco/ 


(Single  Mode) 


Substitute  in  Rayleigh’s  Criterion 

=6il+‘^=sinS,] 

A  one-dimensional  map  with  bifurcation  parameter  C2. 
2)  “Twice-Kicked”  Oscillator 
Physical  Basis:  (observations) 

(i)  A  vortex  is  shed  when  the  velocity  at  the  dump  plane  is 
rising  during  each  cycle. 

(ii)  For  ’small’  vortices,  combustion  occurs  later  in  the 
cycle,  nearly  to  completion,  producing  a  broad  ‘pulse’ 
of  pressure. 


(iii)  For  ‘large’  vortices,  combustion  is  not 
complete  in  the  first  pulse  and  occasionally 
smaller  secondary  pulses  of  pressure  occur. 
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INCORPORATING  A  GLOBAL  BIFURCATION  OF 


COMBUSTION  IN  ANALYSIS  OF  COMBUSTOR 


DYNAMICS  (cont’d) 

Model:  The  energy  released  is  modeled  as  two  pulses 
(5-functions)  in  time.  A  large  amplitude  pulse 
always  occurs  when  the  pressure  is  near 
maximum,  and  rising,  causing  growth.  For 
Targe’  amplitudes  a  pulse  also  occurs  when  the 
pressure  is  minimum,  tending  to  alternate  the 
oscillations. 

Realization: 

Use  Rayleigh’s  Criterion  to  construct  the  relations 
for  energies  at  maxima  and  minima  of  pressure: 


Free  Oscillation: 

®  =$  . 


ma.x+  mm- 


Pulse  near  » 


&  +C.& 

max+  max-  1  ma.x- 


Pulse  near  p-: 


S  .  =  S'  .  “  c,  S  . 

mm+  min-  2  mm- 


Combine  and  renormalize: 


Bifurcation  parameter 
C3  for  C4  fixed 


INCORPORATING  A  GLOBAL  BIFURCATION  OF 
COMBUSTION  IN  ANALYSIS  OF  COMBUSTOR 

DYNAMICS  (cont’d) 

The  bifurcation  diagrams  for  the  two  cases  are  shown  in  the 
following  two  figures.  In  each  case,  period-doubling 
bifurcations  occur,  leading  to  chaotic  behavior  not  observed 
in  the  experimental  results. 

DELAYED  COMBUSTION  RESPONSE 
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INCORPORATING  A  GLOBAL  BIFURCATION  OF 
COMBUSTION  IN  ANALYSIS  OF  COMBUSTOR 

DYNAMICS  (cont’d) 


•  The  question  remains:  can  these  ideas  be  adapted  in  such  a 
fashion  as  to  provide  an  explanation  for  the  subcritical 
bifurcation  and  hysteresis  shown  by  the  experimental  results? 

V.  Burnley  (1996)  “Nonlinear  Combustion  Instabilities 
and  Stochastic  Sources,”  Ph.D.  Thesis,  California  Institute 
of  Technology, 

•  Modeling  the  time-dependent  heat  release  for  distributed 
combustion  following  the  models  of  Levine  and  Baum(1983) 
and  Greene  (1990). 

m' -ii^  or  Q' =Qq^F^u-uj) 


•  The  threshold  velocity  is  necessary  to  give  a 
subcritical  bifurcation. 
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INCORPORATING  A  GLOBAL  BIFURCATION  OF 
COMBUSTION  IN  ANALYSIS  OF  COMBUSTOR 

DYNAMICS  (cont’d) 


•  With  the  above  representation  of  m  or  Q ,  application  of  a 
continuation  method  gives  the  following  bifurcation  diagram: 


BIFURCATION  DIAGRAM  FOR  THE  DETERMINISTIC  SYSTEM; 
THRESHOLD  VELOCITY  MODEL  OF  NONLINEAR  COMBUSTION 


•  For  reasons  not  known,  the  usual  supercritical  bifurcation 
occurs  at  a,  =  0,  followed  by  two  turning  points,,  giving  the 
possibility  for  hysteresis  and  nonlinear  (pulsed)  instabilities. 


Note:  existence  of  a  threshold  velocity  0 
is  essential  for  this  result:  the  reasons 
are  not  known. 
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V  SOME  CONSEQUENCES  OF 
SUBCRITICAL  BIFURCATIONS  WITH  NOISE 


5.1  REPRESENTING  NOISE  IN  ANALYSIS  OF 
THE  GLOBAL  DYNAMICS 

•  There  are  three  sorts  of  problems: 

•  Formal  construction  of  noise  (stochastic)  sources  in 
the  framework  of  spatial  averaging  equations  for  the 
fluctuating  motions  (oscillator  equations). 

•  Modeling  the  noise  sources 

•  Solving  the  stochastic  differential  equations 


A  Construction  of  Noise  Sources  in  the  Oscillator  Equations 

•  The  first  step  is  application  of  the  principle  of  splitting 
small  disturbances. 

Chu,  B.T.  and  Kovazonay,  L.S.G.  (1958)  “Nonlinear 
Interactions  in  a  Viscous  Heat-Conducting  Compressible 
Gas,”  Journal  of  Fluid  Mechanics,  V.  3,  No.  5  (pp.  494- 
514). 
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5.1  REPRESENTING  NOISE  IN  ANALYSIS  OF 
THE  GLOBAL  DYNAMICS  (cont’d) 

SPLITTING  THE  UNSTEADY  FLOW  FIELD 


•  Chu  and  Kovazsnay  (1958) 

•  In  the  limit  of  small  amplitudes,  disturbances  in  a 
compressible  fluid  may  be  synthesized  of  three  types  of 
waves:  acoustic  waves,  vorticity  waves,  and  entropy 
waves. 

•  acoustic  waves  carry  pressure  and  velocity 

•  vorticity  waves  carry  vorticity  and  velocity 

•  entropy  waves  carry  entropy  and  velocity 

P^=P[ 

n'=Q;, 

a  5 

•  To  second  order,  acoustic  waves  remain  free  of  entropy  and 
vorticity. 

•  Apply  the  principle  of  additive  splitting  to  problems  of 
second-order  acoustics  and  nonuniform  mean  flow. 
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5.1  REPRESENTING  NOISE  IN  ANALYSIS  OF 
THE  GLOBAL  DYNAMICS  (cont’d) 

•  Apply  the  principle  of  additive  splitting  to  problems  of 
second-order  acoustics  and  nonimiform  mean  flow.  Write  the 
flow  variables  as  sums  of  the  three  contributions: 

P  =  P>  p's 

To  first  order; 

Acoustic  Waves :  p'  =  p[  p’^=p[=0 
Vorticity Waves:  Q'=Q[=0 

Entropy  W  aves :  s' 5'  =  5^^  =  0 


The  velocity  field  contains  all  three: 

u' —U 

Q  S 

Note:  The  vorticity  field  may  contain,  to  first  order, 
random  contributions  (noise)  and  terms  arising, 
for  example,  from  generation  of  vorticity  at  a 
lateral  surface  with  mass  injection. 
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5.1  REPRESENTING  NOISE  IN  ANALYSIS  OF 
THE  GLOBAL  DYNAMICS  (cont’d) 


A.  Procedure: 

(i)  Write  p=p+p' ,  u=u+u'  etc.  in  the  general 
conservation  equations  expanded  to  3rd  order 
in  O' 

(ii)  Form  the  nonlinear  wave  equation  for  the 
pressure 

(hi)  Assume  and  spatial  average  to 

give  the  oscillator  equations  h,  =F„ : 

a'  a'  a 

+<ffp— dS-  fS-Vv  dr+\  f^\|/  dV 
"  dt  •'  d’’dt 

/,  =  ^•Wu+u'-Vu]-W\\f  ^dV 


etc 


(iv)  Apply  the  splitting  principle  and  rewrite 
i=i:+i?+i: 


•1  "I 


=  j[v(w  •t7')-w'xVxw}V\}/^JF 


427 


5.1  REPRESENTING  NOISE  IN  ANALYSIS  OF 


THE  GLOBAL  DYNAMICS  (cont’d) 


(v)  The  fluctuations  appearing  in  /j,  ...  must  be 
modeled: 

Assume  to  first  order, 

yk- 

and  set  p'„=p',=0 

To  second  order  the  oscillator  equations  have  the  form 


h.  +0);r|.  =2a,f|,  +2(O.0.ti. 

+Sfch,+LT]+H.+(2=:)1 


t].,  are  stochastic  sources  defined  as 

integrals  but  they  must  be 
modeled. 


B.  Modeling  the  Noise  Sources 


•  No  modeling  based  on  experimental,  theoretical,  or 
phenomenological  results  has  been  accomplished. 

•  Approximate  the  sources  and  as  white 
noise  processes. 


5.1  REPRESENTING  NOISE  IN  ANALYSIS  OF 
THE  GLOBAL  DYNAMICS  (cont’d) 

C.  Solving  the  Stochastic  Differential  Equations 

•  The  random  character  of  the  stochastic  sources  will 
appear  as  randomness  imposed  in  the  amplitudes 
r|„(0  of  the  modes. 

•  Thus  the  Fourier  synthesis  of  the  pressure  field 
serves  in  the  first  instance  as  an  approximate 
representation  of  the  deterministic  (acoustics)  part  of 
the  unsteady  motions  and  good  results  can  be 
obtained  by  retaining  only  a  small  number  of  terms 
in  the  modal  expansion. 

•  Results  have  been  obtained  only  for  the  simplest  case 
of  two  modes,  with  noise  sources  only  in  the 
fundamental  mode.  Nonlinear  gasdynamic  coupling 
transfers  stochastic  behavior  to  the  second  mode. 
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5.1  REPRESENTING  NOISE  IN  ANALYSIS  OF 
THE  GLOBAL  DYNAMICS  (cont’d) 

C.  Solving  the  Stochastic  Differential  Equations  (cont’d) 

•  There  are  two  types  of  stochastic  influences 
included  here: 

(i)  represent  stochastic  influences  on  the 
‘spring’  of  natural  frequency  and  in  the 
damping  of  growth  rate  of  the  n^^  mode; 

(ii)  represents  a  stochastic  driving  source 
causing  excitation  of  the  n^^  mode  even  in  the 
absence  of  driving  by  combustion  processes 
[cf.  Russian  interpretation  of  test  results  and 
analysis  by  R.  Messier  for  solid  rockets  in  the 
1970’s-  1980’s. 


The  equations  to  be  solved  are  (Burnley,  1996) 

Ti,  +coJri,  =2(  +2co/ 0,  +(T;  )'^  +E:, 

V  2  J  y  Z(D^  J 

fi,  +C0^T1,  =  2aX  +2co,e,ri3  +{F^  f" 
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Normalized  Spectrum  Pressure  at  the  head  end 


5.1  REPRESENTING  NOISE  IN  ANALYSIS  OF 
THE  GLOBAL  DYNAMICS  (cont’d) 

•  Results  for  probability  density  functions  have  been 
constructed  using  the  Monte-Carlo  method 
(Paparrizos,  1991;  Burnley,  1996)  applied  to  a 
Stratonovich  form  of  the  stochastic  differential 
equations  (Burnley,  1996). 

•  The  following  figures  show  a  pressure  trace  and 
spectrum  for  one  simulation. 


0.2r - 1 - ' - r 


Nondimensional  Tim©,  ta/L 


Nondimensiona!  Frequency,  fUa 


•  Many  simulations  with  initial  conditions  varied, 
and  constructing  probability  distributions  after 
the  initial  transients  have  decayed. 


5.1  REPRESENTING  NOISE  IN  ANALYSIS  OF 
THE  GLOBAL  DYNAMICS  (contM) 


C.  Solving  the  Stochastic  Differential  Equations  (cont’d) 

•  A  series  of  simulations  have  been  carried  out  with  the  nonlinear 
combustion  model  described  above,  and  the  growth  (decay) 
constant  a  ,  of  the  fundamental  mode  varied.  The  probability 
density  functions  are  shown  in  the  following  figures. 
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VI  CONCLUDING  REMARKS 


•  Indisputable  theoretical  and  experimental  evidence 
exists  for  the  existence  of  various  sorts  of  subcritical  as 
well  as  supercritical  bifurcations  in  combustors. 

•  A  framework  is  available  for  investigating  the 
consequences  of  global  bifurcations  in  combustor 
dynamics. 

•  The  consequences  of  many  local  subcritical  bifurcations 
have  not  been  examined. 

•  The  proposal  that  hysteresis  associated  with 

recirculation  zones  and  different  stability  of  the 
associated  steady  combustion  fields  merits  further 
investigation  as  a  possible  explanation  for  random 
occurrences  of  combustion  instabilities. 

•  Numerical  simulations  of  these  problems  remain  in  a 
primitive  state. 
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Some  Aspects  of  Nonreproducibility  of  Self-Oscillation 
Excitation 

in  Combustion  Chamber  of  Liquid  Rocket  Engine 

K.P.  Denisov*,  V.N.  Ivanov**,  V.P.  Pikalov*** 

Research  Institute  of  Chemical  Machine  Building  (NIICHIMMASH), 

Sergiev-Posad,  Russia 

Abstract 

Some  aspects  of  low  reproducibility  of  high-frequency  self-oscillation 
excitation  in  combustion  chamber  of  liquid  rocket  engine. 

In  the  investigation  of  model  combustion  chambers  it  is  shown  that  the 
observed  effects  of  unreproducibility  arise  from: 

•  possibility  of  realization  of  combustion  zone  ambiguous  (bistable) 
stationary  states  with  invariant  parameters  of  the  propellant  feeding  system.; 

•  combustion  zone  nonline*ar  sensitivity  to  chamber  natural  acoustic 
field  that  results  in  the  occurrence  of  self-oscillation  unstable  limit  cycle  of  one  or 
other  mode. 


*  Vice-General  Director. 

**  Senior  Reseach  Engineer. 

***  Head,  Reseach  Division. 
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Nomenclature 


A  ^=root  mean  square  of  total  pulsating  signal 

F=frequency 

H=potential  barrier 

=static  combustion  chamber  pressure 
F  =pulsating  pressure  component 

=  root  mean  square  of  pulsating  pressure  component 
=stationary  probability  density  function 
Q'  -  pulsating  heat  release  component 
P'  =  F  complex  amplitude 
Q'  =  Q'  complex  amplitude 

R(T)=autocorrelation  function 
SLC=stable  limit  cycle 
ULC=  unstable  limit  cycle 
T=oscillation  period 
Vo=Oxygen  velocity  stream 
Vp=Propan  velocity  stream 

a=equivalense  ratio  (oxygen  excess  relative  to  the  stoichiometric  ratio) 

P=decay  coefficient 

pT  =effective  oscillation  decrement 

5(T)=Dirac  delta  function 

X,=amplitude  (envelope)  of  narrow-band  signal 

A-slc“SLC  amplitude 

A,ulc=ULC  amplitude 

^(t)=normal  delta-correlated  (white)  noise 

T=temporal  interval 

(p=oscillation  phase 

=potencial  function 

co=  oscillation  angular  frequency 

(jt)o=  natural  angular  frequency 
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Introduction 


Dynamic  ambiguity  of  combustion  regime  in  liquid  rocket  engine  (LRE) 
chambers  manifests  itself  in  spontaneous  transitions  from  a  state  of  small  random 
oscillation  ("noises')  to  a  state  of  regular  self-oscillation  of  high  amplitude  with 
propellant  supply  condition  and  other  monitored  "external"  factors  being  constant. 
A  distinctive  feature  of  those  transitions  is  the  fact  that  they  occur  occasionally, 
for  no  apparent  reason.  No  observable  disturbance  (of  pulse  or  explosive  nature, 
for  example)  generally  precedes  the  self-excitation  event.  Self-oscillation 
excitation  is  of  "soft",  smooth  character.  The  self-oscillation  termination  occurs 
in  a  similar  way.  The  random  character  of  such  spontaneous  transitions  causes  us 
to  speak  of  self-oscillation  excitation  nonreproducibility  (such  terms  as  "limited 
reproducibility",  "random  reproducibility'  are  also  used). 

At  least  two  basic  concepts  exist  of  the  nature  of  low  reproducibility  of 
acoustic  self-oscillation  excitation  in  combustion  chambers. 

One  of  the  concepts  is  based  on  postulating  a  possibility  of  combustion 
zone  space  pattern  bifurcations  that  are  significant  with  respect  to  stability  to 
small  acoustic  disturbance.  This  is  equivalent  to  the  stating  that  a  stationary 
regime  of  propellant  feeding  to  the  chamber  is  not  a  factor  that  defines  the  flame 
pattern  unambiguously. 

Thus  the  given  concept  is  based  on  the  assumption  that,  in  principle, 
generation  of  multiple  combustion  zone  space  patterns  that  differ  in  regimes  of 
flame  stabilization  is  possible  at  identical  values  of  parameters  defining  the 
combustion  chamber  stationary  mode.  These  combustion  zone  patterns  may  have 
different  sensitivity  to  small  acoustic  disturbances. 

In  particular,  the  realization  of  one  or  other  combustion  zone  space 
pattern  may  depend  on  a  combination  of  random  conditions  that  take  place  in  the 
course  of  engine  starting  and  operating  in  the  specified  regime.  As  a  consequence 
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of  the  random  character  of  one  or  other  combustion  zone  space  pattern  realization, 
the  self-oscillation  excitation  process  acquires  a  statistic  character  also. 

The  second  concept  of  nonreproducibility  is  based  on  the  idea  that  a 
nonlinear  ("hard')  excitation  of  a  potentially  self-oscillating  system  occurs  by  the 
action  of  a  stationary  wide-band  acoustic  noise  of  turbulent  combustion,  or  by  the 
action  of  single  random  pressure  disturbances  of  explosive  nature  that  occur  in 
combustion  chamber  sometimes.  This  model  postulates  the  presence  of  unstable 
limit-cycle  (ULC)  self-oscillation  that  results  from  the  fact  that  n,  the  effective 
index  of  interaction  between  oscillatory  components  Q'  of  heat  release  (or  mass 
release)  rate  and  oscillatory  components  P'  of  pressure  is  a  nonlinear  function  of 
pressure  oscillation  amplitude 


where  Q  and  P  are  complex  amplitude  of  Q'  and  P' ,  respectively. 

The  considered  concepts  of  nonreproducibility  are  alternate  in  a  sense 
because  they  lead  to  different  scenarios  of  random  beginning  of  self-oscillations. 
The  adequacy  of  used  diagnostic  model  is  very  important  for  successful 
application  of  methods  and  tools  of  diagnostics  of  combustion  stability  in  LRE. 
In  this  coimection  the  study  of  the  nature  of  self-oscillation  excitation  low 
reproducibility,  as  well  as  the  search  of  diagnostics  methods  and  ways  of 
predicting  this  phenomenon  are  of  critical  importance.  Some  results  of  the  study 
of  self-oscillation  excitation  low  reproducibility  in  combustion  chamber  are 
presented.  The  data  are  based  on  experimental  research  work  carried  out  at 
NnCHIMMASH's  set-ups  modeling  oscillatory  combustion. 

1.  External  Symptoms  of  Spontaneous  Combustion 

instability 

A  wide  variety  of  external  symptoms  of  spontaneous  combustion 
instability  are  observed  in  engine  tests  of  which  the  following  two  are  most 
common: 
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a)  instability  is  observed  in  a  series  of  m  identical  tests  of  engines  of  the 
same  type  in  n>m  cases  immediately  or  in  some  random  time  t  after  the  engine 
reaches  the  specified  stationary  operation  regime  (Fig. la) 

b)  instability  occurs  in  some  time  after  a  certain  engine  reaches  the 
specified  operation  regime  and  manifests  itself  as  a  single  "spike"  of  self¬ 
oscillation  (Fig.  lb)  or  as  alternating  'noise'  and  self-oscillatory  combustion 
regimes.  (Fig.lc) 

2.  Model  Combustion  Chamber  with  a  Single 
injection  Element 

Investigation  was  carried  out  in  a  model  combustion  chamber  with  a 
single  coaxial  injection  element  (Fg.2).  Gaseous  propane  and  gaseous  oxygen 
were  used  as  propellants.  The  combustion  chamber  is  equipped  with  instruments 
for  flame  photographing  and  reaction  zone  oscillatory  pressure  measurement. 

The  following  result  was  obtained.  In  the  course  of  tests  a  distinct 
ambiguity  of  flame  space  pattern  was  recorded  in  a  certain  range  of  operating 
governing  parameters:  fuel  flow  rate  and  oxidizer  flow  rate.  Depending  on  the 
ignition  condition,  flame  stabilization  regime  (and  flame  space  pattern, 
accordingly)  could  be  in  one  of  two  stable  stationary  state  (Fig.3).  First  of  them 
corresponded  to  the  flame  stabilization  inside  the  injection  element  in  the  vicinity 
of  the  fuel  post  edge  (Scheme  1).  The  flame  had  an  elongated  form  and  radiated  a 
relatively  low  level  of  noise  (1...30  Pa  in  the  0...10  kHz  range).  The  second 
stationary  state  corresponded  to  flame  stabilization  outside  of  the  injection 
element  at  some  distance  fi'om  the  oxidizer  tube  exit  (Scheme  2).  This  state  is 
characterized  by  a  shorter  smeared  flame  and  a  relatively  high  integral  noise  level 
(26...  90  Pa).  The  noise  frequency  contents  were  different  in  both  cases.  The 
mentioned  ambiguity  of  stabilization  mode  (bistability)  was  observed  in  a  wide 
range  of  the  governing  parameters  (fuel  flow  velocity  and  oxidizer  flow  velocity). 

A  principally  important  effect  was  found  in  the  tests:  an  acoustically 
stable  'noise'  combustion  corresponded  to  the  flame  stabilization  Scheme  1 
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throughout  the  investigated  range  of  the  governing  parameters  (  Fig  4a).  In  the 
case  of  flame  stabilization  according  to  Scheme  2,  in  the  same  range  of  the 
governing  parameters  there  was  a  boundary  beyond  which  strong  self-oscillation 
of  the  first  radial  mode  (of  2620  Hz  ,  320  Pa)  excited.  (Fig.4b). 

This  series  of  experiments  showed  that  the  bistability  of  combustion  zone 
pattern  could  really  give  rise  to  dynamic  ambiguity  that  manifests  itself  in  the 
presence  or  absence  of  self-oscillation.  With  unmonitored  ignition  process  (if  it  is 
unknown  what  particular  flame  stabilization  regime  is  realized),  the  excitation  of 
self-oscillation  in  an  ensemble  of  similar  tests  becomes  of  a  random,  statistical 
nature.  Thus  it  can  be  stated  that  the  ambiguity  of  dynamic  state  of  combustion 
may  a  consequence  of  a  stochastic  process  of  initial  ignition  and  subsequent 
generation  of  one  or  another  flame  self-consistent  spatial  pattern. 

It  should  be  noted  that  in  the  given  experiments,  once  the  mixture  was 
ignited  and  flame  stabilization  regime  was  established  as  per  Scheme  1,  an 
artificial  transition  from  Scheme  1  to  Scheme  2  was  achieved  by  disturbing 
oxygen  concentration  in  the  reaction  zone.  Some  air  was  injected  in  the  oxygen 
feeding  line  for  a  short  time  for  this  purpose.  (Fig.5).  After  the  air  injection  was 
terminated,  the  system  retained  stabilization  Scheme  2.  In  this  case,  after  the 
system  transition  to  the  region  of  potential  instability,  self-oscillation  excitation 
became  a  completely  deterministic,  reproducible  phenomenon. 

Thus  in  this  example  the  self-oscillation  excitation  nonreproducibility  is 
only  apparent  and  is  a  result  of  the  absence  of  information  on  flame  stabilization 
regime  on  ignition.  However  the  mere  fact  of  recognizing  a  possibility  of 
bifurcation  of  flame  spatial  pattern  with  subsequent  self-oscillation  excitation 
deserves  attention. 

3.  Model  Combustion  Chamber  with  Gas-Liquid 
injection  Eiements 

A  combustion  chamber  of  104  mm  in  diameter  was  investigated  (Fig.6). 
Propellants,  liquid  kerosene  and  gaseous  (cold)  oxygen  were  supplied  from 
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cylinders.  The  chamber  was  equipped  with  an  independent  igniter,  using  methane 
and  oxygen.  A  device  for  applying  artificial  pulse  disturbances  to  the  reaction 
zone  and  fuel  feed  cavity  is  moimted  on  the  chamber. 

The  injection  head  includes  four  injection  elements  with  radial  input  of 
liquid  fuel  into  the  flow  of  gaseous  oxidizer. 

The  range  of  test  regimes  on  coordinates:  chamber  static  pressure,  - 
equivalence  ratio  (oxygen  excess  relative  to  the  stoichiometric  ratio),  a,  is  shown 
in  Fig.  7. 

In  the  course  of  testing  this  injection  head  type  the  following  dynamic 
state  of  combustion  process  were  recorded; 

•  noise  stable  state,  (Pe<  0.4  MPa),  Fig.8; 

•  chaotic  alternation  of  noise  and  self-oscillatory  states.  Each 
dynamic  state  occurrence  instants  and  duration  are  random  in  the  ensemble  of 
realizations  similar  by  P^  and  a,  (Pc=  0.4  MPa),  Fig.9; 

•  self-oscillatory  combustion  (f,o,  =  6300  Hz,  P^  =  0.6... 0.7  MPa), 

Fig.l0. 

Thus  it  can  be  stated  that  in  the  P  =  0.4. .  .0.6  MPa,  a  =  0.65 . .  .0.75  range 
the  dynamic  state  of  the  combustion  process  is  characterized  by  stochastic 
ambiguity  and  the  ambiguity  symptoms  differ  fi-om  those  described  in  the 
previous  section. 

To  understand  the  observed  transition  nature,  let  us  consider  a  model  of 
potentially  self-oscillatory  system  being  in  the  mode  of  hard  excitation. 

3.1.  Model  of  potentially  self-oscillatory  system  in  the  mode 

of  hard  excitation. 

It  should  be  noted  that  in  the  case  under  consideration  the  self-oscillation 
is  excited  at  a  fi'equency  of  combined  first  longitudinal  mode  -first  tangential 
mode  (f,o]  =  6300  Hz,).  Therefore  a  simple  one-mode  approximation  will  be  used 
in  devising  a  model  of  dynamic  transitions.  In  such  approximation  the  combustion 
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chamber  stochastic  dynamics  can  be  described  by  the  following  model  of 
potentially  self-oscillatory  system 

ipi 

^  +  2y9(l)—  -f  (1) 

at  at 

where:  P'  is  oscillatory  pressure; 

t  is  time; 

(Oq  is  a  circular  frequency  ; 

|3(X)  is  an  effective  coefficient  of  decay  or  generation,  a  nonlinear 
fimction  of  amplitude  X‘, 

^(t)  is  a  stationary  normal  delta-correlated  (white)  noise  of  turbulent 
combustion. 

With  the  proviso  that  ’the  following  evolutionary  equation 

for  pressure  oscillation  amplitude  X  corresponds  to  equation  (1)*^^: 

^  =  -/3(A)  •  A  +  a),N, /SA  +  ^(0-  (2). 

at 

Here:  Nq  is  the  spectral  intensity  of  random  noise  ^(t); 

(t)  is  normal  delta-correlated  random  function: 

R(t)=<  at)  ■  at  +  r)  >=  •  N,  ■  S(r)IA. 

Within  the  framework  of  the  self-oscillatory  system  with  hard  excitation 
a  particular  P(X,)  relationship  is  postulated  which  provides  the  existence  of  self¬ 
oscillation  unstable  limit  cycle  and  stable  limit  cycle. 

Evolutionary  equation  (2)  describes  one-dimension  Markov  process  with 
coefficients  of  "drift"  and  "diffusion",  equations  (3)  and  (4),  respectively: 


a(A)  =  —P(A)  •  A  -t-  q}q  NqI^A, 

(3) 

b(A)  =  co,^NJA. 

(4) 

Equation  (2)  has  the  potential  function  in  the  following  form: 
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(5) 


where: 


A  ^ ) 


=0)q  -N^/ 8 j3o. 

(6) 

Pq  =  4t  A  — ^  0 , 

(7) 

and  a  stationary  state  probability  density 


P„(A)  =  C,  exp[-20(A)/or']  =  QAexp 


-^]Aj3(A)dZ 


(8). 


The  C,  constant  is  determined  from  the  norming  condition 


jpjA)dA  =  L 


0 


The  calculated  relations  0(x)  and  Pjt  are  shown  in  Fig  1 1  (x=  Mr). 

As  can  be  seen,  0(x)  potential  has  two  minima  that  correspond  to  locally 
stable  state  of  dynamic  equilibrium:  the  left  minimum  corresponds  to  the  noise 
state  while  the  right  one  corresponds  to  the  self-oscillatory  state.  Unstable 
dynamic  equilibrium  corresponds  to  the  unstable  limit  cycle. 

The  distribution  of  state  probability  density  P(x)  for  a  specific  form  of 
function  P  =  p(x)  is  bimodal  and  the  maxima  of  probability  density  correspond  to 
the  minima  of  potential  0(x).  The  probability  density  minimum  corresponds  to 
the  maximum  of  potential  0(x),  i.e.  to  the  unstable  limit  cycle. 

The  self-excitation  of  the  dynamic  system  described  by  Eq.  (1)  is 
stochastic  in  principle,  in  the  sense  that  it  is  impossible  to  know  a  priori  the  exact 
time  instant  when  excitation  will  occur  i.  e.  the  time  of  system  transition  from 
noise  to  self-oscillatory  state. 

The  mean  expectation  time  of  self-excitation  event  depends  on  the  noise 
intensity  No  and  the  potential  barrier  height  h,.  If  the  potential  barrier  is 
sufficiently  high  (which  corresponds  to  a  high  amplitude  of  the  imstable  limit 
cycle  compared  to  ct,  Eq.  (6)),  and  the  noise  intensity  Nq  is  small,  the  expectation 
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time  can  be  very  large,  and  self-oscillation  excitation  will  be  a  very  rare  event, 
accordingly.  And  vice  versa,  with  low  potential  barrier  and  high  noise  intensity 
random  transitions  of  "noise  < — >self-oscillation"  will  occur.  In  this  case 
pressure  oscillogram  for  i-realization  will  look  like  alternating  noise  and  self¬ 
oscillation  regions,  similar  to  those  shown  in  Fig.9. 

Thus  within  the  limits  of  Model  1,  the  nonreproducibility  of  self- 
oscillation  excitation,  with  'external'  conditions  being  constant,  is  attributed  to  the 
presence  of  unstable  limit  cycle  and  random  wide-band  noise  ^(t). 

3.2  System  Identification 

The  task  of  identification  of  a  potentially  self-oscillatory  system  is  stated 
in  this  case  as  follows:  to  what  extent  the  initial  uncertainty  of  the  statistic  model 
of  the  observed  dynamic  transitions  can  be  reduced. 

The  model  of  potentially  self-oscillatory  system  with  hard  excitation 
(Model  1)  is  taken  as  a  base  model.  The  following  relationship  connecting  P(A,) 
with  the  narrow-band  noise  amplitude  probability  density  is  obtained: 

where 

'y(ji)  =  ,  (10) 

dX 

and  Pst  (A.)  is  a  function  of  amplitude  probability  density,  estimated  experimentally 
based  on  the  realization  of  the  random  process  A,(t). 

The  realization  of  oscillatory  process  P'  =P'(t)  filtered  in  the  band 
/= 1000...  10000  Hz  and  the  corresponding  noise  spectrogram  for  the  region  that 
precedes  oscillation  self-excitation  (t  =4... 5  s)  are  shown  in  Fig.  12  (a  compressed 

oscillogram  is  presented).  The  narrow-band  peak  with  decrement  fiT  »  0.04 
observed  in  the  vicinity  of/,o,  «  6150  Hz  corresponds  to  combined  first  tangential 
mode-first  longitudinal  mode  of  the  chamber  volume  natural  acoustical 
oscillations. 
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The  estimate  of  probability  density  of  envelope  (amplitude)  X,(t)  of 
narrow-band  noise  P'(t)  =X(t)-cos  {27tfji-^(p(t)')  in  the  region  that  precedes 
oscillation  self-excitation  is  shown  in  Fig.  13.  The  corresponding  estimate  of  the 
decay  coefficient  [P]  -  amplitude  [  A,]  relation.  As  can  be  seen  in  Fig.  13,  the 
reconstructed  relation  P=P(X)  is  of  a  falling  character  and  predicts  the  existence  of 
unstable  limit  cycle  with  amplitude  X  ulc  *  8ct.  The  found  fact  of  unstable  limit 
cycle  existence  allows  identifying  the  system  being  studied  within  the  base  model 
of  a  potentially  self-oscillatory  system  with  hard  excitation  of  self-oscillation 
(Model  1). 

It  should  be  noted  that  for  the  case  of  dissipative  system  of  the  second 
order,  P(X)=P„=const,  the  distribution  (8)  changes  to  the  known  Rayleigh 
distribution  for  the  amplitude  (envelope)  of  linear  narrow-band  noise: 


In  the  fi'amework  of  Model  1,  actually  deterministic  transition  fi'om  a 
locally  stable  noise  state  into  a  self-  oscillatory  one  can  be  performed  by  applying 
a  pulse  to  the  system  (combustion  zone).  The  pulse  should  be  high  sufficiently  to 
overcome  the  potential  barrier,  (hard  pulse  excitation  of  the  system).  In  this  case 
oscillation  with  an  amplitude  of  X  >X  ^lc  should  develop  in  the  system  after  the 
disturbance  application. 

In  the  same  time.  Model  1  assumes  a  possibility  of  a  reverse  transition 
from  a  self-oscillatory  state  to  noise  in  case  the  pulse  is  applied  as  antiphase  to  the 
existent  self-oscillation. 

Such  situation  is  shown  in  Fig  14  for  the  chamber  under  investigation. 
The  pulse  disturbance  was  applied  to  the  chamber  reaction  zone  by  exploding  an 
explosive  charge. 

It  should  be  mentioned  however  that  treatment  of  effects  of  self¬ 
oscillation  excitation  and  termination  under  the  action  of  artificial  disturbance  in 
the  framework  of  Model  1  is  not  quite  correct.  There  is  reason  to  believe  that, 
under  the  conditions  of  the  described  experiments,  the  artificial  pulse  disturbance. 
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apart  from  acoustic  effect,  produce  a  parametric  effect  on  combustion  (it  changes 
the  combustion  zone  pattern,  possibly  including  change  in  the  flame  stabilization 
regime).  This  is  suggested,  in  particular,  by  a  specific  character  of  the  dynamic 
transition  of  "noise  to  self-oscillation":  the  decay  of  oscillatory  response  actually 
to  the  initial  noise  level,  with  a  subsequent  smooth  transition  to  the  stable  limit 
cycle.  Fig.  15.  A  possibility  of  self-oscillation  termination  on  application  of  a  pulse 
to  the  fuel  feed  cavity  is  another  evidence  of  parametric  effects  (Fig.  16).  The 
establishment  of  a  new  'noise'  stationary  state  on  pulse  application  is  also  shown 
in  Fig.  16.  This  state  differs  in  noise  spectrum  and  level  from  the  state  that  preceds 
oscillation  self-excitation. 

Thus  the  given  series  of  experiments  reveals  a  possibility  of  realization  of 
combustion  zone  multiple  stationary  patterns,  which  is  heyond  the  frames  of  the 
base  Model  1. 

The  treatment  presented  shows  that  while  devising  diagnostic  models  of 
oscillatory  combustion  in  LRE  chambers  it  is  necessary  to  consider,  along  with 
'hard  excitation'  effects,  a  possibility  of  bifurcation  of  combustion  stationary 
modes  and  dynamic  state  ambiguity  resulted  from  the  bifurcation.. 

Conclusion 

From  the  conducted  studies  it  is  concluded  that  at  least  two  aspects 
should  be  taken  into  consideration  for  devising  adequate  models  of  unstable 
combustion  in  LRE: 

•  possibility  of  realization  of  combustion  zone  ambiguous  (bistable) 
stationary  states  with  invariant  parameters  of  the  propellant  feeding  system.; 

•  combustion  zone  nonlinear  sensitivity  to  chamber  natural  acoustic 
field  that  results  in  the  occurrence  of  self-oscillation  unstable  limit  cycle  of  one  or 
other  mode. 

Both  aspects  of  realization  of  flame  self-consistent  pattern  may  lead  to 
the  nonreproducibility  of  self-oscillation  excitation  in  LRE  combustion  chamber. 
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Pressure 

oscillations 


Fig.  3.  Photograph  of  the  diffusion  flame. 

1 .  Flame  stabilization  inside  the  injection  element 
in  the  vicinity  of  the  fuel  post  edge. 

2.  Flame  stabilization  outside  of  the  injection  element 
above  the  oxidizer  tube. 


450 


V„  (m/s) 


V„(m/s) 

30 


20 


10 


Fig.  4.  Dynamic  .states  and  path  of  the  system  movement  in  the  domain  ot  govetTiing  parameters: 

a)  flame  stabilization  as  per  Scheme  1; 

b)  flame  stabilization  as  per  Scheme  2; 

(propellant  feeding  parameters,  Vo  and  Vp ,  maiiced  with  similar  figures  arc  completely  identical). 
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Flame  stabilization  as  Flame  stabilization  as 

per  Scheme  I  \  per  Scheme  2 
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Fig.  6.  Schematic  of  combustion  instability  study 
(all  dimensions  are  in  millimeters) 


Static  pressure  (MPa) 


Equivalence  ratio 


Fig.  7.  Combustion  chamber  operation  regimes. 
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Noise  state 


Amplitude  spectram 
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Fig.  8.  Noise  state  of  combustion  process. 
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Fig.  10.  Self-oscillatory  combustion. 
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Fig.  1 1 .  Deca\  coefficienl  (?.  probability  density  Pst  and  potential  (P 
versus  oscillation  relative  amplitude. 


h  1-2  -  Potential  barrier  of  the  transition  from  state  1  to  state  2 
h  2-1  -  Potential  barrier  of  the  transition  from  state  2  to  state  1 
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Fig.  12.  Spontaneous  dynamic  transition  “noise  self-oscillation”; 

a)  compressed  oscillogram; 

b)  noise  .spectrum  in  the  region  preceding  oscillation  sell  e.vcitation. 
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Disturbance 


Sell-oscillation 
spontaneous  teimi  nation 
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Time 


Disturbance 

Self-oscillation 

spontaneous  excitation  ^ 
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Time 


Fig.  14.  Self-oscillation  excitation  (a)  and  termination  (b)  artificial 
disturbance  (disturbance  applied  to  the  chamber  reaction  zone). 
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l  ig.  16.  Self-oscillalion  termination  (disturbance  applied  to  lied  cavity). 
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Fundamental  Combustion  Instability  Research 
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40K  L0X/CH4  Stability  testing 
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Unique  Combustion  Stability  Validation  Database 
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"Coaxial  Injector  Spray  Characterization  Using  Water/Air  as  Simulants" 
Michelle  Zaller,  Mark  klem,  28Th  JANNAF  Combustion  Meeting, 
Volume  2,  pp151-160,  1993. 
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Gaseous  and  Liquid  Nitrogen  Gaseous  and  Liquid  Nitrogen 

No  swirl  No 

Gaseous  Nitrogen;  T  =  SST,  .006  tb/s  Gaseous  Nitrogen:  T  =  55°F,  .006  tb/s 


Vaporizing  Liquid  Nitrogen  Spray  Data  for  Varying  Chamber  Pressures 
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Radial  Location  (in)  ttadial  Location  (iti) 

Liquid  nitrogen  flow  rale  .024  Ib/s;  gaseous  nitrogen  flow  rate  .012-.015  Ib/s. 


LO 


>.  -L. 


—  n 
0) 
CO 

o 

/  rs 


> 


o; 

N 


5^. 

X 


-D 

O  .  ‘d 
rr 


iT: 


!t  h 
cr  cr* 


i 

■H 

;  o 


(S^ 


_L 


o 

1 

ID 

O 

i 

d 

p 

1 

in 

1; 

— i 

Ll, 

1 

d 

c 

r> 

i 

> — ^ 

i 

X 

1 

o 

• 

d 

Q 

i/) 

1 

j 

u 

V-‘ 

N^. 

— - 

— 

H 

O 

<X 

1 

.9 

•D 

q 

D 

q 

a: 

• 

d 

] 

d 

1 

o 

d 

o  00  ^  ^ 

(s/uj>  A:|IDCj0A  dOJQ  IDIXV  UOB]/^ 


i 

cn  in 

d  d 

• 

li  li 

O'  O' 

O  O 

— 

d  d 

d  d 

1)  II 

tj-  u_  “ 

C  O 

X  X 

- 

O  O  -J 

-J  -J  } 

*» 

i 

V)  kfi 

W 

» 

h 

j 

.•-N 

\j 

n  :C 

-  CD 

q  q 

i 

1 

U 

1 

•/♦ 

: 

/ 

C  • 

/ 

'4 

/ 

/ 

-  •r' 

i 

o'o 

o 

O 

o 

o 

lO 

O 

m 

c 

CM 

CM 

6 


tD 

6 


0.5 

c 

c 

o 

d 

“c 

c 

0 

tn 

_J 

d 

c 

*0 

CM 

cr 

d 

o 

d 


fLun')  j?^?luD!G  UD?1^ 


s.  n! 

d 


< 

■  r . 

C£i 

Lu 

Th 

(/) 

LJ 

cr 

□_ 

q: 

LU 

CD 

2 

< 

( 

'd 


c 

'X' 


<!'• 

fv 


475 


TEST  PARAMETERS 

FOR  PHASE/DOPPLER  DROP  SIZE  MEASUREMENTS 
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Rgure  8. — Case 
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Figure  9. — Power 
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Rgure  1 2.— Variation  of  osciHation  amplitude  with  mixture  ratio 
(configuration  3,  Df  •  .235'*). 
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High  Energy  Laser  as  a  Stability  Rating  Device 
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3.  45  -  degree  mirror 

4.  Focusing  lens 

5.  Test  combustor 


C02  Laser 

Synrad  57-1  sealed  gas,  RF  excitation 


O 

> 

O 

> 

D 

■u 

o 

LO 


■(3  ”5  N 


o 


£ 

E 


0) 

c 

o 


(D 

m 

o 


N 

X 

o 

I 

^  I 


£  X 


ur> 


c  N 

E  X 

CM 


O 


489 


SUMMARY 

Despite  far  from  optimal  conditions  (unsteady  combustion,  inability 
to  move  the  focal  point)  we  were  consistently  able  to  drive  very  small 
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Comparison  of  multiple-exposure  (time-averaged)  with  instantaneous  photographs. 


Ultrafast  Imaging  of  Sprays 
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Probe  beam  :  1054  nm;  Pump  beam  :  527  nm 
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Intensity  (arb.  units) 
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cooperation  in  the  field  of 
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ANALYSIS  OF  THE  FULL  3D  FLOW  IN  CENTRIFUGAL 

IMPELLER 


o 

(/) 

o 

& 

*S 

o 


O  (/)  o 
o  c  «*=- 
c  o 

W  ‘^3 

i  §i 

O  0*  0 

/» 

^■rp 


CL 

Q. 

(0 


0) 


a=  o 
§.« 


.ii  P  .  > 

m*  •■■  ■—■ 


(0 

3 

(0 

> 

0) 

o 


(/) 

S 

o 

75 

(Q 


0) 

> 

O 

0) 

o 


o  ci>  !5 

£  >  w 

C  Z  (D 

°Q  I 

W  CO  c 

5=0 

©  3 

E“- 

2  -d  75 
IQ  o  -- 
Q.  Q.  o 

r-  J2  S 

E  o>  t: 

.2^  >  CO 

(/)  o 

O  *0 

■o 

^  =  O) 

s  ®  •- 

©  JQ  -D 

©  =3 

©  73 


©  .£  ^ 


<  ■D 


:=  © 


E't^  ;©'/|l 

E  o 

s  ©  ©  ^ 
^  Tr  *=  s 

g  ©  ©  **- 

©  **  O  #A 
©  T3  2 
C  o  “O  © 
O  £  ©  = 
o  ^  ©  © 
.±  ©  3  a 

« 1-2.1 

3  ©  I. 

cr  E  © 

3  > 

=  1  ° 

© 

^  © 

'+■<  1(M  _ 

‘C  '?  H  ^ 

>  .E  .  c 

>  O)  © 

£  ©  c  S 

©  3>‘q.%: 

^  Qm  "Q 

S  »  i5 

o  ^  ® 

^  ^  ©  "O 

©  E  © 

o:  ^  c 

^  ©  .— 

w  >  S 

58  o  'S  -o 

^  P  A  *5: 

El.  © 

©  o 

ilil 

.E  n  Q£  c 


^  © 

*_  O 
o  « 

s 

c  o 
o  s: 

-fS  o 


©  I4S 


© 

© 

o 

c 

© 


© 

© 


Q. 

E 

3 

Q. 


E  n 
8- 
©  ©  L. 

£  ©  O 
K  >  t 
£  © 

©  ©  ^ 
■Q  3  ® 

-5  ©  © 
©  £  Q. 

E  3) 
.J2  3  c 
© 

m  © 

«  C  "5 

5^2 

-O  Q. 

©  ©  *0 

S’®® 

©  ©  rr 

©  © 
C  ■=  CL 

2  o  ® 

Q.^  O 

^  ©  © 
©  o  .2 


3 

73 

© 

3 


S 

(§1 

© 

u 

© 

•  • 

s 


529 


FLOW  INSIDE  IMPELLER  BLADES 
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Grid  Structure  for  Centrifugal  Impeller  Blades 


Comparison  of  Numerical  and  Experimental 
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EXPERIMENTAL  AND  NUMERICAL  ANALYSIS  OF  FALLING 

FILM  BEHAVIOR 
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NTAL  ANALYSIS  OF  LOUVERED  FIN  HEAT 
EXCHANGERS  PERFORMANCE 
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EFFICIENCY  ENHANCEMENT  OF  THE  COAL  FIRED  DOMESTIC  HEATERS 

AND  STOVES 


Domestic  Downdraft  Coal  Burner 
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INFLUENCE  OF  HORIZONTAL  FINS  WITH  LOW 
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FORCED  CONVECTIVE  HEAT  TRANSFER  IN  POROUS  MEDIUM  OF  WIRE 
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16) Copper  plate, 

1 7) heaters 
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second  medium. 


LOW-SPEED  WIND  TUNNEL 
800x800  mm,  10  m/s  <  Velocity  <30  m/s 


HIGH-SPEED  WIND  TUNNEL 
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150x150  mm 

0.6  <  Mach  <  2.2  and  2.3  <  Mach  <  4.0 


HIGH-SPEED  WIND  TUNNEL 
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Institute  (SAGE,  TOrkiye) 
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UPERSONIC  FLOWS 
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Pitot-tube  and  Pressure-based 
Instrumentations 


<  short  view  of  experiments 
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Photographs  of  supersonic  free  jets 
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Shadowgraphs  of  a  supersonic 
jet  impingement  with 
1 0  degrees  inclination 


Rayleigh  measurement  of  temperature  in  the  reaction 
zone  of  a  combustion  flpwityt^ithout  flame 
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AUTOMOTIVE  LABORATORIES 

E-mail ;  sorusbav  @  inkn.itu.edu.tr 
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investigated 
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AUTOMATION  TECHNOLOGIES 
DEVELOPMENT  AND  EDUCATION  CENTER 
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Tachometers,  Power  Supplies,  etc.) 


